From 40e5a5811c6ffce9b0974e93cdd927cbcf60c157 Mon Sep 17 00:00:00 2001
From: Joe Hunkeler <jhunkeler@gmail.com>
Date: Tue, 11 Aug 2015 16:51:37 -0400
Subject: Repatch (from linux) of OSX IRAF

---
 sys/pmio/plprop.gx | 177 +++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 177 insertions(+)
 create mode 100644 sys/pmio/plprop.gx

(limited to 'sys/pmio/plprop.gx')

diff --git a/sys/pmio/plprop.gx b/sys/pmio/plprop.gx
new file mode 100644
index 00000000..60d1c603
--- /dev/null
+++ b/sys/pmio/plprop.gx
@@ -0,0 +1,177 @@
+# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
+
+include	<plset.h>
+include	<plio.h>
+
+# PL_PIXROP -- Rasterop between source and destination pixel arrays.
+
+procedure pl_pixrop$t (px_src,xs,src_maxval, px_dst,ds,dst_maxval, npix, rop)
+
+PIXEL	px_src[ARB]		#I source pixel array
+int	xs			#I starting pixel index in src
+int	src_maxval		#I max pixel value in src mask
+PIXEL	px_dst[ARB]		#O destination pixel array
+int	ds			#I starting pixel index in dst
+int	dst_maxval		#I max pixel value in dst mask
+int	npix			#I number of pixels to convert
+int	rop			#I rasterop
+
+pointer	sp, src
+int	opcode, i
+PIXEL	data, ceil, src_value
+int	and(), or(), xor(), not()
+define	out_ 91
+
+begin
+	opcode = R_OPCODE(rop)
+	data   = R_DATA(rop)
+	ceil   = 0
+
+	# Pixel value to be used if input mask is boolean.
+	if (src_maxval == 1) {
+	    src_value = data
+	    if (src_value <= 0)
+		src_value = dst_maxval
+	}
+
+	# Handle the easy cases first.
+	switch (opcode) {
+	case PIX_CLR:
+	    call aclr$t (px_dst[ds], npix)
+	    return
+	case PIX_SET:
+	    call amovk$t (data, px_dst[ds], npix)
+	    goto out_
+	case PIX_SRC:
+	    if (src_maxval != 1)
+		call amov$t (px_src[xs], px_dst[ds], npix)
+	    else {
+		do i = 1, npix
+		    if (px_src[xs+i-1] > 0)
+			px_dst[ds+i-1] = src_value
+		    else
+			px_dst[ds+i-1] = 0
+	    }
+
+	    goto out_
+	case PIX_DST:
+	    return	# no-op
+	}
+
+	# Integer or boolean source mask?
+	if (src_maxval != 1) {
+	    # Integer source mask; operate directly on source mask.
+
+	    switch (opcode) {
+	    case PIX_NOTSRC:
+		do i = 1, npix
+		    px_dst[ds+i-1] = not (px_src[xs+i-1])
+	    case PIX_NOTDST:
+		do i = 1, npix
+		    px_dst[ds+i-1] = not (px_dst[xs+i-1])
+
+	    case PIX_SRC_AND_DST:
+		do i = 1, npix
+		    px_dst[ds+i-1] = and (px_src[xs+i-1], px_dst[ds+i-1])
+	    case PIX_SRC_OR_DST:
+		do i = 1, npix
+		    px_dst[ds+i-1] =  or (px_src[xs+i-1], px_dst[ds+i-1])
+	    case PIX_SRC_XOR_DST:
+		do i = 1, npix
+		    px_dst[ds+i-1] = xor (px_src[xs+i-1], px_dst[ds+i-1])
+
+	    case PIX_SRC_AND_NOTDST:
+		do i = 1, npix
+		    px_dst[ds+i-1] = and (px_src[xs+i-1], not(px_dst[ds+i-1]))
+	    case PIX_SRC_OR_NOTDST:
+		do i = 1, npix
+		    px_dst[ds+i-1] =  or (px_src[xs+i-1], not(px_dst[ds+i-1]))
+	    case PIX_NOTSRC_AND_DST:
+		do i = 1, npix
+		    px_dst[ds+i-1] = and (not(px_src[xs+i-1]), px_dst[ds+i-1])
+	    case PIX_NOTSRC_OR_DST:
+		do i = 1, npix
+		    px_dst[ds+i-1] =  or (not(px_src[xs+i-1]), px_dst[ds+i-1])
+
+	    case PIX_NOT_SRC_AND_DST:
+		do i = 1, npix
+		    px_dst[ds+i-1] = not (and (px_src[xs+i-1], px_dst[ds+i-1]))
+	    case PIX_NOT_SRC_OR_DST:
+		do i = 1, npix
+		    px_dst[ds+i-1] = not ( or (px_src[xs+i-1], px_dst[ds+i-1]))
+	    case PIX_NOT_SRC_XOR_DST:
+		do i = 1, npix
+		    px_dst[ds+i-1] = not (xor (px_src[xs+i-1], px_dst[ds+i-1]))
+	    }
+
+	} else {
+	    # Boolean source mask; use integer DATA value from ROP if source
+	    # mask pixel is set.
+
+	    call smark (sp)
+	    call salloc (src, npix, TY_PIXEL)
+
+	    do i = 1, npix
+		if (px_src[xs+i-1] > 0)
+		    Mem$t[src+i-1] = src_value
+		else
+		    Mem$t[src+i-1] = 0
+
+	    switch (opcode) {
+	    case PIX_NOTSRC:
+		do i = 1, npix
+		    px_dst[ds+i-1] = not (Mem$t[src+i-1])
+	    case PIX_NOTDST:
+		do i = 1, npix
+		    px_dst[ds+i-1] = not (px_dst[xs+i-1])
+
+	    case PIX_SRC_AND_DST:
+		do i = 1, npix
+		    px_dst[ds+i-1] = and (Mem$t[src+i-1], px_dst[ds+i-1])
+	    case PIX_SRC_OR_DST:
+		do i = 1, npix
+		    px_dst[ds+i-1] =  or (Mem$t[src+i-1], px_dst[ds+i-1])
+	    case PIX_SRC_XOR_DST:
+		do i = 1, npix
+		    px_dst[ds+i-1] = xor (Mem$t[src+i-1], px_dst[ds+i-1])
+
+	    case PIX_SRC_AND_NOTDST:
+		do i = 1, npix
+		    px_dst[ds+i-1] = and (Mem$t[src+i-1], not(px_dst[ds+i-1]))
+	    case PIX_SRC_OR_NOTDST:
+		do i = 1, npix
+		    px_dst[ds+i-1] =  or (Mem$t[src+i-1], not(px_dst[ds+i-1]))
+	    case PIX_NOTSRC_AND_DST:
+		do i = 1, npix
+		    px_dst[ds+i-1] = and (not(Mem$t[src+i-1]), px_dst[ds+i-1])
+	    case PIX_NOTSRC_OR_DST:
+		do i = 1, npix
+		    px_dst[ds+i-1] =  or (not(Mem$t[src+i-1]), px_dst[ds+i-1])
+
+	    case PIX_NOT_SRC_AND_DST:
+		do i = 1, npix
+		    px_dst[ds+i-1] = not (and (Mem$t[src+i-1], px_dst[ds+i-1]))
+	    case PIX_NOT_SRC_OR_DST:
+		do i = 1, npix
+		    px_dst[ds+i-1] = not ( or (Mem$t[src+i-1], px_dst[ds+i-1]))
+	    case PIX_NOT_SRC_XOR_DST:
+		do i = 1, npix
+		    px_dst[ds+i-1] = not (xor (Mem$t[src+i-1], px_dst[ds+i-1]))
+	    }
+
+	    call sfree (sp)
+	}
+out_
+	# If writing to an integer mask, mask the data to the indicated max
+	# value (necessary to avoid very large values if any NOT operations
+	# occurred).  If writing to a boolean mask, map positive integer mask
+	# values to 1.
+
+	if (dst_maxval == 1) {
+	    data = 1
+	    call argt$t (px_dst[ds], npix, ceil, data)
+	} else if (dst_maxval > 1) {
+	    data = dst_maxval
+	    call aandk$t (px_dst[ds], data, px_dst[ds], npix)
+	}
+end
-- 
cgit