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/amaxk.gx | 29 +++++++++++++++++++++++++++++
 1 file changed, 29 insertions(+)
 create mode 100644 sys/vops/amaxk.gx

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

diff --git a/sys/vops/amaxk.gx b/sys/vops/amaxk.gx
new file mode 100644
index 00000000..f45bca09
--- /dev/null
+++ b/sys/vops/amaxk.gx
@@ -0,0 +1,29 @@
+# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
+
+# AMAXK -- Compute the maximum of a constant and a vector (generic).
+
+procedure amaxk$t (a, b, c, npix)
+
+PIXEL	a[ARB]
+PIXEL	b
+PIXEL	c[ARB]
+int	npix, i
+$if (datatype == x)
+real	abs_b
+$endif
+
+begin
+	$if (datatype == x)
+	abs_b = abs (b)
+	$endif
+
+	do i = 1, npix
+	    $if (datatype == x)
+		if (abs(a[i]) >= abs_b)
+		    c[i] = a[i]
+		else
+		    c[i] = b
+	    $else
+		c[i] = max (a[i], b)
+	    $endif
+end
-- 
cgit