From fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 Mon Sep 17 00:00:00 2001
From: Joseph Hunkeler <jhunkeler@gmail.com>
Date: Wed, 8 Jul 2015 20:46:52 -0400
Subject: Initial commit

---
 sys/imio/tf/impl3r.x | 51 +++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 51 insertions(+)
 create mode 100644 sys/imio/tf/impl3r.x

(limited to 'sys/imio/tf/impl3r.x')

diff --git a/sys/imio/tf/impl3r.x b/sys/imio/tf/impl3r.x
new file mode 100644
index 00000000..675fd8d0
--- /dev/null
+++ b/sys/imio/tf/impl3r.x
@@ -0,0 +1,51 @@
+# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
+
+include	<syserr.h>
+include	<imhdr.h>
+include	<imio.h>
+
+# IMPL3? -- Put a line to an apparently three dimensional image.  If there
+# is only one input buffer, no image section, we are not referencing out of
+# bounds, and no datatype conversion needs to be performed, directly access
+# the pixels to reduce the overhead per line.
+
+pointer procedure impl3r (im, line, band)
+
+pointer	im		# image header pointer
+int	line		# line number within band
+int	band		# band number
+
+int	fd, nchars
+long	vs[3], ve[3], offset
+pointer	bp, impgsr(), fwritep()
+errchk	imopsf, imerr
+
+begin
+	repeat {
+	    if (IM_FAST(im) == YES && IM_PIXTYPE(im) == TY_REAL) {
+		fd = IM_PFD(im)
+		if (fd == NULL) {
+		    call imopsf (im)
+		    next
+		}
+		if (line < 1 || line > IM_LEN(im,2) ||
+		    band < 1 || band > IM_LEN(im,3))
+		    call imerr (IM_NAME(im), SYS_IMREFOOB)
+
+		offset = (((band - 1) * IM_PHYSLEN(im,2) + line - 1) *
+		    IM_PHYSLEN(im,1)) * SZ_REAL + IM_PIXOFF(im)
+		nchars = IM_PHYSLEN(im,1) * SZ_REAL
+		ifnoerr (bp = (fwritep (fd, offset, nchars) - 1) / SZ_REAL + 1)
+		    return (bp)
+	    }
+
+	    vs[1] = 1
+	    ve[1] = IM_LEN(im,1)
+	    vs[2] = line
+	    ve[2] = line
+	    vs[3] = band
+	    ve[3] = band
+
+	    return (impgsr (im, vs, ve, 3))
+	}
+end
-- 
cgit