diff options
author | Joe Hunkeler <jhunkeler@gmail.com> | 2015-08-11 16:51:37 -0400 |
---|---|---|
committer | Joe Hunkeler <jhunkeler@gmail.com> | 2015-08-11 16:51:37 -0400 |
commit | 40e5a5811c6ffce9b0974e93cdd927cbcf60c157 (patch) | |
tree | 4464880c571602d54f6ae114729bf62a89518057 /sys/mwcs/gen/mwltranr.x | |
download | iraf-osx-40e5a5811c6ffce9b0974e93cdd927cbcf60c157.tar.gz |
Repatch (from linux) of OSX IRAF
Diffstat (limited to 'sys/mwcs/gen/mwltranr.x')
-rw-r--r-- | sys/mwcs/gen/mwltranr.x | 26 |
1 files changed, 26 insertions, 0 deletions
diff --git a/sys/mwcs/gen/mwltranr.x b/sys/mwcs/gen/mwltranr.x new file mode 100644 index 00000000..9cafe4d2 --- /dev/null +++ b/sys/mwcs/gen/mwltranr.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_ltranr (p1, p2, ltm, ltv, ndim) + +real p1[ndim] #I input point +real p2[ndim] #O transformed output point +real ltm[ndim,ndim] #I linear transformation matrix +real ltv[ndim] #I linear translation vector +int ndim #I dimension of system + +int i, j +real p3[MAX_DIM] + +begin + call amovr (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 |