Initial commit

1 files changed, 48 insertions, 0 deletions

diff --git a/noao/digiphot/daophot/daolib/mvmul.x b/noao/digiphot/daophot/daolib/mvmul.x
new file mode 100644
index 00000000..352def99
--- /dev/null
+++ b/noao/digiphot/daophot/daolib/mvmul.x
@@ -0,0 +1,48 @@
+# MVMUL -- Multply a matrix (left-hand side) by a one dimensional vector
+# (right-hand side) and return the resultant vector.
+
+procedure mvmul (matrix, maxdim, dim, vector, result)
+
+real	matrix [maxdim, maxdim]	# input matrix
+int	maxdim			# maximum size of input matrix
+int	dim			# dimension of matrix and vectors
+real	vector[maxdim]		# input vector
+real	result[maxdim]		# iutput vector
+
+double	sum
+int	i, j
+
+begin
+	do i = 1, dim {
+	    sum = 0.0
+	    do j = 1, dim
+		sum = sum + double (matrix[j,i]) * double(vector[j])
+	    result[i] = sum
+	}
+
+end
+
+
+# DMVMUL -- Multply a matrix (left-hand side) by a one dimensional vector
+# (right-hand side) and return the resultant vector.
+
+procedure dmvmul (matrix, maxdim, dim, vector, result)
+
+double	matrix [maxdim, maxdim]	# input matrix
+int	maxdim			# maximum size of input matrix
+int	dim			# dimension of matrix and vectors
+double	vector[maxdim]		# input vector
+double	result[maxdim]		# iutput vector
+
+double	sum
+int	i, j
+
+begin
+	do i = 1, dim {
+	    sum = 0.0d0
+	    do j = 1, dim
+		sum = sum + (matrix[j,i] * vector[j])
+	    result[i] = sum
+	}
+
+end