From fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 Mon Sep 17 00:00:00 2001 From: Joseph Hunkeler Date: Wed, 8 Jul 2015 20:46:52 -0400 Subject: Initial commit --- sys/vops/lz/arcpi.x | 24 ++++++++++++++++++++++++ 1 file changed, 24 insertions(+) create mode 100644 sys/vops/lz/arcpi.x (limited to 'sys/vops/lz/arcpi.x') diff --git a/sys/vops/lz/arcpi.x b/sys/vops/lz/arcpi.x new file mode 100644 index 00000000..193f35e1 --- /dev/null +++ b/sys/vops/lz/arcpi.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 arcpi (a, b, c, npix) + +int a # constant numerator +int b[ARB] # vector denominator +int c[ARB] # output vector +int npix +int i + +begin + if (a == 0) { + call aclri (c, npix) + } else if (a == 1) { + do i = 1, npix + c[i] = 1 / b[i] + } else { + do i = 1, npix + c[i] = a / b[i] + } +end -- cgit