1 files changed, 20 insertions, 0 deletions
diff --git a/sys/mwcs/gen/mwvmuld.x b/sys/mwcs/gen/mwvmuld.x
new file mode 100644
index 00000000..0af8dfa7
--- /dev/null
+++ b/sys/mwcs/gen/mwvmuld.x
@@ -0,0 +1,20 @@
+# MW_VMUL -- Vector multiply.
+
+procedure mw_vmuld (a, b, c, ndim)
+
+double	a[ndim,ndim]		#I input matrix
+double	b[ndim]			#I input vector
+double	c[ndim]			#O output vector
+int	ndim			#I system dimension
+
+int	i, j
+double	v
+
+begin
+	do j = 1, ndim {
+	    v = 0
+	    do i = 1, ndim
+		v = v + a[i,j] * b[i]
+	    c[j] = v
+	}
+end