aboutsummaryrefslogtreecommitdiff
path: root/sys/vops/lz/asqrd.x
diff options
context:
space:
mode:
authorJoe Hunkeler <jhunkeler@gmail.com>2015-08-11 16:51:37 -0400
committerJoe Hunkeler <jhunkeler@gmail.com>2015-08-11 16:51:37 -0400
commit40e5a5811c6ffce9b0974e93cdd927cbcf60c157 (patch)
tree4464880c571602d54f6ae114729bf62a89518057 /sys/vops/lz/asqrd.x
downloadiraf-osx-40e5a5811c6ffce9b0974e93cdd927cbcf60c157.tar.gz
Repatch (from linux) of OSX IRAF
Diffstat (limited to 'sys/vops/lz/asqrd.x')
-rw-r--r--sys/vops/lz/asqrd.x23
1 files changed, 23 insertions, 0 deletions
diff --git a/sys/vops/lz/asqrd.x b/sys/vops/lz/asqrd.x
new file mode 100644
index 00000000..e6cf3f70
--- /dev/null
+++ b/sys/vops/lz/asqrd.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 asqrd (a, b, npix, errfcn)
+
+double a[ARB], b[ARB]
+int npix, i
+extern errfcn()
+double errfcn()
+errchk errfcn
+
+begin
+ do i = 1, npix {
+ if (a[i] < 0)
+ b[i] = errfcn (a[i])
+ else
+ {
+ b[i] = sqrt (a[i])
+ }
+ }
+end