diff options
| author | Joseph Hunkeler <jhunkeler@gmail.com> | 2015-07-08 20:46:52 -0400 |
|---|---|---|
| committer | Joseph Hunkeler <jhunkeler@gmail.com> | 2015-07-08 20:46:52 -0400 |
| commit | fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 (patch) | |
| tree | bdda434976bc09c864f2e4fa6f16ba1952b1e555 /sys/mwcs/gen/mwvmuld.x | |
| download | iraf-linux-fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4.tar.gz | |
Initial commit
Diffstat (limited to 'sys/mwcs/gen/mwvmuld.x')
| -rw-r--r-- | sys/mwcs/gen/mwvmuld.x | 20 |
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 |