diff options
| author | Joe Hunkeler <jhunkeler@gmail.com> | 2015-08-11 16:51:37 -0400 |
|---|---|---|
| committer | Joe Hunkeler <jhunkeler@gmail.com> | 2015-08-11 16:51:37 -0400 |
| commit | 40e5a5811c6ffce9b0974e93cdd927cbcf60c157 (patch) | |
| tree | 4464880c571602d54f6ae114729bf62a89518057 /sys/vops/lz/allnx.x | |
| download | iraf-osx-40e5a5811c6ffce9b0974e93cdd927cbcf60c157.tar.gz | |
Repatch (from linux) of OSX IRAF
Diffstat (limited to 'sys/vops/lz/allnx.x')
| -rw-r--r-- | sys/vops/lz/allnx.x | 23 |
1 files changed, 23 insertions, 0 deletions
diff --git a/sys/vops/lz/allnx.x b/sys/vops/lz/allnx.x new file mode 100644 index 00000000..b4527117 --- /dev/null +++ b/sys/vops/lz/allnx.x @@ -0,0 +1,23 @@ +# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc. + +# ALLN -- Compute the natural logarithm of a vector (generic). If the natural +# logarithm is undefined (x <= 0) a user supplied function is called to get +# the pixel value to be returned. + +procedure allnx (a, b, npix, errfcn) + +complex a[ARB], b[ARB] +int npix, i +extern errfcn() +complex errfcn() +errchk errfcn + +begin + do i = 1, npix { + if (a[i] == (0.0,0.0)) + b[i] = errfcn (a[i]) + else { + b[i] = log (a[i]) + } + } +end |