aboutsummaryrefslogtreecommitdiff
path: root/sys/vops/abnek.gx
diff options
context:
space:
mode:
Diffstat (limited to 'sys/vops/abnek.gx')
-rw-r--r--sys/vops/abnek.gx31
1 files changed, 31 insertions, 0 deletions
diff --git a/sys/vops/abnek.gx b/sys/vops/abnek.gx
new file mode 100644
index 00000000..4643cd89
--- /dev/null
+++ b/sys/vops/abnek.gx
@@ -0,0 +1,31 @@
+# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
+
+# ABNEK -- Vector boolean not equals constant. C[i], type INT, is set to 1 if
+# A[i] is not equal to B, else C[i] is set to zero.
+
+procedure abnek$t (a, b, c, npix)
+
+PIXEL a[ARB]
+PIXEL b
+int c[ARB]
+int npix
+int i
+
+begin
+ # The case b==0 is perhaps worth optimizing. On many machines this
+ # will save a memory fetch.
+
+ if (b == 0$f) {
+ do i = 1, npix
+ if (a[i] != 0$f)
+ c[i] = 1
+ else
+ c[i] = 0
+ } else {
+ do i = 1, npix
+ if (a[i] != b)
+ c[i] = 1
+ else
+ c[i] = 0
+ }
+end