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

---
 pkg/xtools/ranges/rgpack.gx | 37 +++++++++++++++++++++++++++++++++++++
 1 file changed, 37 insertions(+)
 create mode 100644 pkg/xtools/ranges/rgpack.gx

(limited to 'pkg/xtools/ranges/rgpack.gx')

diff --git a/pkg/xtools/ranges/rgpack.gx b/pkg/xtools/ranges/rgpack.gx
new file mode 100644
index 00000000..d77a3a09
--- /dev/null
+++ b/pkg/xtools/ranges/rgpack.gx
@@ -0,0 +1,37 @@
+# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
+
+include	<pkg/rg.h>
+
+# RG_PACK -- Pack input data to include only points in the ranges.
+#
+# The output array must be large enough to contain the desired points.
+# If the ranges are merged then the input and output arrays may be the same.
+
+procedure rg_pack$t (rg, a, b)
+
+pointer	rg					# Ranges
+PIXEL	a[ARB]					# Input array
+PIXEL	b[ARB]					# Output array
+
+int	i, j, k, n
+
+begin
+	# Error check the range pointer.
+
+	if (rg == NULL)
+	    call error (0, "Range pointer undefined")
+
+	j = 0
+	do i = 1, RG_NRGS(rg) {
+	    if (RG_X1(rg, i) > RG_X2(rg, i)) {
+		do k = RG_X1(rg, i), RG_X2(rg, i), -1 {
+		    j = j + 1
+		    b[j] = a[k]
+		}
+	    } else {
+		n = RG_X2(rg, i) - RG_X1(rg, i) + 1
+	        call amov$t (a[RG_X1(rg, i)], b[j + 1], n)
+	        j = j + n
+	    }
+	}
+end
-- 
cgit