diff options
| author | Joseph Hunkeler <jhunkeler@gmail.com> | 2015-07-08 20:46:52 -0400 |
|---|---|---|
| committer | Joseph Hunkeler <jhunkeler@gmail.com> | 2015-07-08 20:46:52 -0400 |
| commit | fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 (patch) | |
| tree | bdda434976bc09c864f2e4fa6f16ba1952b1e555 /sys/mwcs/mwmmul.gx | |
| download | iraf-linux-fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4.tar.gz | |
Initial commit
Diffstat (limited to 'sys/mwcs/mwmmul.gx')
| -rw-r--r-- | sys/mwcs/mwmmul.gx | 23 |
1 files changed, 23 insertions, 0 deletions
diff --git a/sys/mwcs/mwmmul.gx b/sys/mwcs/mwmmul.gx new file mode 100644 index 00000000..b86449c2 --- /dev/null +++ b/sys/mwcs/mwmmul.gx @@ -0,0 +1,23 @@ +# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc. + +# MW_MMUL -- Matrix multiply. + +procedure mw_mmul$t (a, b, c, ndim) + +PIXEL a[ndim,ndim] #I left input matrix +PIXEL b[ndim,ndim] #I right input matrix +PIXEL c[ndim,ndim] #O output matrix +int ndim #I dimensionality of system + +int i, j, k +PIXEL v + +begin + do j = 1, ndim + do i = 1, ndim { + v = 0 + do k = 1, ndim + v = v + a[k,j] * b[i,k] + c[i,j] = v + } +end |