diff options
Diffstat (limited to 'sys/vops/adot.gx')
| -rw-r--r-- | sys/vops/adot.gx | 28 |
1 files changed, 28 insertions, 0 deletions
diff --git a/sys/vops/adot.gx b/sys/vops/adot.gx new file mode 100644 index 00000000..baadd952 --- /dev/null +++ b/sys/vops/adot.gx @@ -0,0 +1,28 @@ +# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc. + +# ADOT -- Vector inner or dot product. The function value is the sum of the +# products of each pair of elements of the input vectors. + +$if (datatype == ld) +double procedure adot$t (a, b, npix) +$else +real procedure adot$t (a, b, npix) +$endif + +PIXEL a[ARB], b[ARB] + +$if (datatype == ld) +double sum +$else +real sum +$endif + +int npix, i + +begin + sum = 0$f + do i = 1, npix + sum = sum + a[i] * b[i] + + return (sum) +end |