diff options
| author | Joe Hunkeler <jhunkeler@gmail.com> | 2015-08-11 16:51:37 -0400 |
|---|---|---|
| committer | Joe Hunkeler <jhunkeler@gmail.com> | 2015-08-11 16:51:37 -0400 |
| commit | 40e5a5811c6ffce9b0974e93cdd927cbcf60c157 (patch) | |
| tree | 4464880c571602d54f6ae114729bf62a89518057 /sys/vops/asum.gx | |
| download | iraf-osx-40e5a5811c6ffce9b0974e93cdd927cbcf60c157.tar.gz | |
Repatch (from linux) of OSX IRAF
Diffstat (limited to 'sys/vops/asum.gx')
| -rw-r--r-- | sys/vops/asum.gx | 32 |
1 files changed, 32 insertions, 0 deletions
diff --git a/sys/vops/asum.gx b/sys/vops/asum.gx new file mode 100644 index 00000000..716d2b53 --- /dev/null +++ b/sys/vops/asum.gx @@ -0,0 +1,32 @@ +# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc. + +# ASUM -- Vector sum. Take care to prevent integer overflow by returning +# a floating point sum. + +$if (datatype == csir) +real procedure asum$t (a, npix) +$else $if (datatype == ld) +double procedure asum$t (a, npix) +$else +PIXEL procedure asum$t (a, npix) +$endif $endif + +PIXEL a[ARB] +int npix +int i + +$if (datatype == csir) +real sum +$else $if (datatype == ld) +double sum +$else +PIXEL sum +$endif $endif + +begin + sum = 0$f + do i = 1, npix + sum = sum + a[i] + + return (sum) +end |