From 40e5a5811c6ffce9b0974e93cdd927cbcf60c157 Mon Sep 17 00:00:00 2001 From: Joe Hunkeler Date: Tue, 11 Aug 2015 16:51:37 -0400 Subject: Repatch (from linux) of OSX IRAF --- sys/mwcs/gen/mwmmuld.x | 21 +++++++++++++++++++++ 1 file changed, 21 insertions(+) create mode 100644 sys/mwcs/gen/mwmmuld.x (limited to 'sys/mwcs/gen/mwmmuld.x') diff --git a/sys/mwcs/gen/mwmmuld.x b/sys/mwcs/gen/mwmmuld.x new file mode 100644 index 00000000..ae35f082 --- /dev/null +++ b/sys/mwcs/gen/mwmmuld.x @@ -0,0 +1,21 @@ +# MW_MMUL -- Matrix multiply. + +procedure mw_mmuld (a, b, c, ndim) + +double a[ndim,ndim] #I left input matrix +double b[ndim,ndim] #I right input matrix +double c[ndim,ndim] #O output matrix +int ndim #I dimensionality of system + +int i, j, k +double 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 -- cgit