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