diff options
| author | Joseph Hunkeler <jhunkeler@gmail.com> | 2015-07-08 20:46:52 -0400 |
|---|---|---|
| committer | Joseph Hunkeler <jhunkeler@gmail.com> | 2015-07-08 20:46:52 -0400 |
| commit | fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 (patch) | |
| tree | bdda434976bc09c864f2e4fa6f16ba1952b1e555 /sys/vops/ak/abeqkx.x | |
| download | iraf-linux-fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4.tar.gz | |
Initial commit
Diffstat (limited to 'sys/vops/ak/abeqkx.x')
| -rw-r--r-- | sys/vops/ak/abeqkx.x | 31 |
1 files changed, 31 insertions, 0 deletions
diff --git a/sys/vops/ak/abeqkx.x b/sys/vops/ak/abeqkx.x new file mode 100644 index 00000000..c2767408 --- /dev/null +++ b/sys/vops/ak/abeqkx.x @@ -0,0 +1,31 @@ +# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc. + +# ABEQK -- Vector boolean equals constant. C[i], type INT, is set to 1 if +# A[i] equals B, else C[i] is set to zero. + +procedure abeqkx (a, b, c, npix) + +complex a[ARB] +complex 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.0,0.0)) { + do i = 1, npix + if (a[i] == (0.0,0.0)) + c[i] = 1 + else + c[i] = 0 + } else { + do i = 1, npix + if (a[i] == b) + c[i] = 1 + else + c[i] = 0 + } +end |