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/lz/arcpr.x | 24 ++++++++++++++++++++++++ 1 file changed, 24 insertions(+) create mode 100644 sys/vops/lz/arcpr.x (limited to 'sys/vops/lz/arcpr.x') diff --git a/sys/vops/lz/arcpr.x b/sys/vops/lz/arcpr.x new file mode 100644 index 00000000..f52a1651 --- /dev/null +++ b/sys/vops/lz/arcpr.x @@ -0,0 +1,24 @@ +# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc. + +# ARCP -- Reciprocal of a constant divided by a vector. No divide by zero +# checking is performed. + +procedure arcpr (a, b, c, npix) + +real a # constant numerator +real b[ARB] # vector denominator +real c[ARB] # output vector +int npix +int i + +begin + if (a == 0.0) { + call aclrr (c, npix) + } else if (a == 1.0) { + do i = 1, npix + c[i] = 1.0 / b[i] + } else { + do i = 1, npix + c[i] = a / b[i] + } +end -- cgit