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/lz/asqrl.x | |
| download | iraf-linux-fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4.tar.gz | |
Initial commit
Diffstat (limited to 'sys/vops/lz/asqrl.x')
| -rw-r--r-- | sys/vops/lz/asqrl.x | 23 |
1 files changed, 23 insertions, 0 deletions
diff --git a/sys/vops/lz/asqrl.x b/sys/vops/lz/asqrl.x new file mode 100644 index 00000000..3b0d23f0 --- /dev/null +++ b/sys/vops/lz/asqrl.x @@ -0,0 +1,23 @@ +# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc. + +# ASQR -- Compute the square root of a vector (generic). If the square root +# is undefined (x < 0) a user supplied function is called to compute the value. + +procedure asqrl (a, b, npix, errfcn) + +long a[ARB], b[ARB] +int npix, i +extern errfcn() +long errfcn() +errchk errfcn + +begin + do i = 1, npix { + if (a[i] < 0) + b[i] = errfcn (a[i]) + else + { + b[i] = sqrt (double (a[i])) + } + } +end |