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/plio/plregrop.x | 76 +++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 76 insertions(+)
 create mode 100644 sys/plio/plregrop.x

(limited to 'sys/plio/plregrop.x')

diff --git a/sys/plio/plregrop.x b/sys/plio/plregrop.x
new file mode 100644
index 00000000..6e37183a
--- /dev/null
+++ b/sys/plio/plregrop.x
@@ -0,0 +1,76 @@
+# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
+
+include <plset.h>
+include <plio.h>
+
+# PL_REGIONROP -- Perform a rasterop operation on an irregular region in a
+# 2-dimensional plane of a mask.  The boundary of the region is defined by
+# a user supplied function which, given the line number (Y) returns a range
+# list defining the region for that line.
+#
+#	rl_new = ufcn (ufd, y, rl_out, xs, npix)
+#
+# where
+#	rl_new	= true if range list for line Y is different than for Y-1
+#	ufd	= user region descriptor (parameters defining region)
+#	y	= input y value
+#	rl_out 	= output range list
+#	xs	= first pixel to be edited in dst mask
+#	npix	= number of pixels in edit region
+#
+# If the dimensionality of the mask exceeds 2, the pl_setplane() procedure
+# should be called first to define the plane of the mask to be modified.
+
+procedure pl_regionrop (pl, ufcn, ufd, y1, y2, rop)
+
+pointer	pl			#I mask descriptor
+extern	ufcn()			#I user supplied region tracing procedure
+pointer	ufd			#I user region descriptor
+int	y1			#I first mask line to be modified
+int	y2			#I last mask line to be modified
+int	rop			#I rasterop defining operation
+
+bool	rl_new
+long	v[PL_MAXDIM]
+int	ll_len, xs, npix
+pointer	sp, ll_out, ll_reg, ll_dst, ol_dst, rl_out
+pointer	pl_access()
+int	pl_r2li()
+bool	ufcn()
+errchk	plvalid
+
+begin
+	call plvalid (pl)
+
+	call smark (sp)
+	call salloc (ll_out, LL_MAXLEN(pl), TY_SHORT)
+	call salloc (ll_reg, LL_MAXLEN(pl), TY_SHORT)
+	call salloc (rl_out, RL_MAXLEN(pl), TY_INT)
+
+	call amovl (PL_PLANE(pl,1), v, PL_MAXDIM)
+	ol_dst = 1
+
+	for (v[2]=y1;  v[2] <= y2;  v[2]=v[2]+1) {
+	    ll_dst = pl_access (pl, v)
+	    rl_new = ufcn (ufd, v[2], Memi[rl_out], xs, npix)
+	    if (rl_new)
+		ll_len = pl_r2li (Memi[rl_out], 1, Mems[ll_reg], npix)
+
+	    if (ll_dst != ol_dst || rl_new) {
+	    call pl_linestencil (Mems[ll_reg],  1, 1,
+				 Mems[ll_dst], xs, PL_MAXVAL(pl),
+				 Mems[ll_reg],  1,
+				 Mems[ll_out], npix, rop)
+		ol_dst = ll_dst
+	    }
+
+	    # Update the affected line of the destination mask.
+	    call pl_update (pl, v, Mems[ll_out])
+	}
+
+	# Compress the mask if excessive free space has accumulated.
+	if (PL_NEEDCOMPRESS(pl))
+	    call pl_compress (pl)
+
+	call sfree (sp)
+end
-- 
cgit