From 40e5a5811c6ffce9b0974e93cdd927cbcf60c157 Mon Sep 17 00:00:00 2001 From: Joe Hunkeler Date: Tue, 11 Aug 2015 16:51:37 -0400 Subject: Repatch (from linux) of OSX IRAF --- sys/vops/ak/abltkx.x | 29 +++++++++++++++++++++++++++++ 1 file changed, 29 insertions(+) create mode 100644 sys/vops/ak/abltkx.x (limited to 'sys/vops/ak/abltkx.x') diff --git a/sys/vops/ak/abltkx.x b/sys/vops/ak/abltkx.x new file mode 100644 index 00000000..04527b7f --- /dev/null +++ b/sys/vops/ak/abltkx.x @@ -0,0 +1,29 @@ +# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc. + +# ABLTK -- Vector boolean less than constant. C[i], type INT, is set to 1 if +# A[i] is less than B, else C[i] is set to zero. + +procedure abltkx (a, b, c, npix) + +complex a[ARB] +complex b +int c[ARB] +int npix +int i +real abs_b + +begin + # The case b==0 is perhaps worth optimizing. On many machines this + # will save a memory fetch. + + if (b == (0.0,0.0)) { + call aclri (c, npix) + } else { + abs_b = abs (b) + do i = 1, npix + if (abs (a[i]) < abs_b) + c[i] = 1 + else + c[i] = 0 + } +end -- cgit