From fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 Mon Sep 17 00:00:00 2001
From: Joseph Hunkeler <jhunkeler@gmail.com>
Date: Wed, 8 Jul 2015 20:46:52 -0400
Subject: Initial commit

---
 sys/vops/asum.gx | 32 ++++++++++++++++++++++++++++++++
 1 file changed, 32 insertions(+)
 create mode 100644 sys/vops/asum.gx

(limited to 'sys/vops/asum.gx')

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
-- 
cgit