From fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 Mon Sep 17 00:00:00 2001 From: Joseph Hunkeler Date: Wed, 8 Jul 2015 20:46:52 -0400 Subject: Initial commit --- sys/mwcs/gen/mwltrand.x | 26 ++++++++++++++++++++++++++ 1 file changed, 26 insertions(+) create mode 100644 sys/mwcs/gen/mwltrand.x (limited to 'sys/mwcs/gen/mwltrand.x') diff --git a/sys/mwcs/gen/mwltrand.x b/sys/mwcs/gen/mwltrand.x new file mode 100644 index 00000000..d35670c7 --- /dev/null +++ b/sys/mwcs/gen/mwltrand.x @@ -0,0 +1,26 @@ +# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc. + +include "../mwcs.h" + +# MW_LTRAN -- Perform a general N-dimensional linear transformation, i.e., +# matrix multiply and translation. + +procedure mw_ltrand (p1, p2, ltm, ltv, ndim) + +double p1[ndim] #I input point +double p2[ndim] #O transformed output point +double ltm[ndim,ndim] #I linear transformation matrix +double ltv[ndim] #I linear translation vector +int ndim #I dimension of system + +int i, j +double p3[MAX_DIM] + +begin + call amovd (p1, p3, ndim) + do j = 1, ndim { + p2[j] = ltv[j] + do i = 1, ndim + p2[j] = p2[j] + ltm[i,j] * p3[i] + } +end -- cgit