From fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 Mon Sep 17 00:00:00 2001 From: Joseph Hunkeler Date: Wed, 8 Jul 2015 20:46:52 -0400 Subject: Initial commit --- pkg/xtools/ranges/rgexcluder.x | 56 ++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 56 insertions(+) create mode 100644 pkg/xtools/ranges/rgexcluder.x (limited to 'pkg/xtools/ranges/rgexcluder.x') diff --git a/pkg/xtools/ranges/rgexcluder.x b/pkg/xtools/ranges/rgexcluder.x new file mode 100644 index 00000000..44cb90fe --- /dev/null +++ b/pkg/xtools/ranges/rgexcluder.x @@ -0,0 +1,56 @@ +# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc. + +include + +# RG_EXCLUDE -- Exclude points given by 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_excluder (rg, a, nin, b, nout) + +pointer rg # Ranges +real a[nin] # Input array +int nin # Number of input points +real b[ARB] # Output array +int nout # Number of output points + +int i, j, k, n, ntemp + +begin + # Error check the range pointer. + + if (rg == NULL) + call error (0, "Range descriptor undefined") + + if (RG_NRGS(rg) == 0) { + call amovr (a[1], b[1], nin) + nout = nin + } else { + ntemp = 0 + + i = 1 + j = 1 + k = min (nin, min (RG_X1(rg, i), RG_X2(rg, i)) - 1) + n = max (0, k - j + 1) + call amovr (a[j], b[ntemp+1], n) + ntemp = ntemp + n + + do i = 2, RG_NRGS(rg) { + j = max (1, max (RG_X1(rg, i-1), RG_X2(rg, i-1)) + 1) + k = min (nin, min (RG_X1(rg, i), RG_X2(rg, i)) - 1) + n = max (0, k - j + 1) + call amovr (a[j], b[ntemp+1], n) + ntemp = ntemp + n + } + + i = RG_NRGS (rg) + j = max (1, max (RG_X1(rg, i), RG_X2(rg, i)) + 1) + k = nin + n = max (0, k - j + 1) + call amovr (a[j], b[ntemp+1], n) + ntemp = ntemp + n + } + + nout = ntemp +end -- cgit