aboutsummaryrefslogtreecommitdiff
path: root/sys/mwcs/mwvmul.gx
diff options
context:
space:
mode:
authorJoseph Hunkeler <jhunkeler@gmail.com>2015-07-08 20:46:52 -0400
committerJoseph Hunkeler <jhunkeler@gmail.com>2015-07-08 20:46:52 -0400
commitfa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 (patch)
treebdda434976bc09c864f2e4fa6f16ba1952b1e555 /sys/mwcs/mwvmul.gx
downloadiraf-linux-fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4.tar.gz
Initial commit
Diffstat (limited to 'sys/mwcs/mwvmul.gx')
-rw-r--r--sys/mwcs/mwvmul.gx22
1 files changed, 22 insertions, 0 deletions
diff --git a/sys/mwcs/mwvmul.gx b/sys/mwcs/mwvmul.gx
new file mode 100644
index 00000000..1c5f4867
--- /dev/null
+++ b/sys/mwcs/mwvmul.gx
@@ -0,0 +1,22 @@
+# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
+
+# MW_VMUL -- Vector multiply.
+
+procedure mw_vmul$t (a, b, c, ndim)
+
+PIXEL a[ndim,ndim] #I input matrix
+PIXEL b[ndim] #I input vector
+PIXEL c[ndim] #O output vector
+int ndim #I system dimension
+
+int i, j
+PIXEL v
+
+begin
+ do j = 1, ndim {
+ v = 0
+ do i = 1, ndim
+ v = v + a[i,j] * b[i]
+ c[j] = v
+ }
+end