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/imio/tf/impakl.x | 46 ++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 46 insertions(+) create mode 100644 sys/imio/tf/impakl.x (limited to 'sys/imio/tf/impakl.x') diff --git a/sys/imio/tf/impakl.x b/sys/imio/tf/impakl.x new file mode 100644 index 00000000..884f931b --- /dev/null +++ b/sys/imio/tf/impakl.x @@ -0,0 +1,46 @@ +# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc. + +# IMPAK? -- Convert an array of pixels of a specific datatype to the +# datatype given as the final argument. + +procedure impakl (a, b, npix, dtype) + +long a[npix] +int b[npix], npix, dtype + +pointer bp + +begin + switch (dtype) { + case TY_USHORT: + call achtlu (a, b, npix) + case TY_SHORT: + call achtls (a, b, npix) + case TY_INT: + if (SZ_INT == SZ_INT32) + call achtli (a, b, npix) + else { + call malloc (bp, npix, TY_INT) + call achtli (a, Memi[bp], npix) + call ipak32 (Memi[bp], b, npix) + call mfree (bp, TY_INT) + } + case TY_LONG: + if (SZ_INT == SZ_INT32) + call achtll (a, b, npix) + else { + call malloc (bp, npix, TY_LONG) + call achtll (a, Meml[bp], npix) + call ipak32 (Meml[bp], b, npix) + call mfree (bp, TY_LONG) + } + case TY_REAL: + call achtlr (a, b, npix) + case TY_DOUBLE: + call achtld (a, b, npix) + case TY_COMPLEX: + call achtlx (a, b, npix) + default: + call error (1, "Unknown datatype in imagefile") + } +end -- cgit