From 40e5a5811c6ffce9b0974e93cdd927cbcf60c157 Mon Sep 17 00:00:00 2001 From: Joe Hunkeler Date: Tue, 11 Aug 2015 16:51:37 -0400 Subject: Repatch (from linux) of OSX IRAF --- math/gsurfit/gsvector.gx | 65 ++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 65 insertions(+) create mode 100644 math/gsurfit/gsvector.gx (limited to 'math/gsurfit/gsvector.gx') diff --git a/math/gsurfit/gsvector.gx b/math/gsurfit/gsvector.gx new file mode 100644 index 00000000..60044dd6 --- /dev/null +++ b/math/gsurfit/gsvector.gx @@ -0,0 +1,65 @@ +# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc. + +include +$if (datatype == r) +include "gsurfitdef.h" +$else +include "dgsurfitdef.h" +$endif + +# GSVECTOR -- Procedure to evaluate the fitted surface at an array of points. +# The GS_NCOEFF(sf) coefficients are stored in the +# vector COEFF. + +$if (datatype == r) +procedure gsvector (sf, x, y, zfit, npts) +$else +procedure dgsvector (sf, x, y, zfit, npts) +$endif + +pointer sf # pointer to surface descriptor structure +PIXEL x[ARB] # x value +PIXEL y[ARB] # y value +PIXEL zfit[ARB] # fits surface values +int npts # number of data points + +begin + # evaluate the surface along the vector + switch (GS_TYPE(sf)) { + case GS_POLYNOMIAL: + if (GS_XORDER(sf) == 1) { + call $tgs_1devpoly (COEFF(GS_COEFF(sf)), y, zfit, npts, + GS_YORDER(sf), GS_YMAXMIN(sf), GS_YRANGE(sf)) + } else if (GS_YORDER(sf) == 1) { + call $tgs_1devpoly (COEFF(GS_COEFF(sf)), x, zfit, npts, + GS_XORDER(sf), GS_XMAXMIN(sf), GS_XRANGE(sf)) + } else + call $tgs_evpoly (COEFF(GS_COEFF(sf)), x, y, zfit, npts, + GS_XTERMS(sf), GS_XORDER(sf), GS_YORDER(sf), GS_XMAXMIN(sf), + GS_XRANGE(sf), GS_YMAXMIN(sf), GS_YRANGE(sf)) + case GS_CHEBYSHEV: + if (GS_XORDER(sf) == 1) { + call $tgs_1devcheb (COEFF(GS_COEFF(sf)), y, zfit, npts, + GS_YORDER(sf), GS_YMAXMIN(sf), GS_YRANGE(sf)) + } else if (GS_YORDER(sf) == 1) { + call $tgs_1devcheb (COEFF(GS_COEFF(sf)), x, zfit, npts, + GS_XORDER(sf), GS_XMAXMIN(sf), GS_XRANGE(sf)) + } else + call $tgs_evcheb (COEFF(GS_COEFF(sf)), x, y, zfit, npts, + GS_XTERMS(sf), GS_XORDER(sf), GS_YORDER(sf), GS_XMAXMIN(sf), + GS_XRANGE(sf), GS_YMAXMIN(sf), GS_YRANGE(sf)) + case GS_LEGENDRE: + if (GS_XORDER(sf) == 1) { + call $tgs_1devleg (COEFF(GS_COEFF(sf)), y, zfit, npts, + GS_YORDER(sf), GS_YMAXMIN(sf), GS_YRANGE(sf)) + } else if (GS_YORDER(sf) == 1) { + call $tgs_1devleg (COEFF(GS_COEFF(sf)), x, zfit, npts, + GS_XORDER(sf), GS_XMAXMIN(sf), GS_XRANGE(sf)) + } else + call $tgs_evleg (COEFF(GS_COEFF(sf)), x, y, zfit, npts, + GS_XTERMS(sf), GS_XORDER(sf), GS_YORDER(sf), GS_XMAXMIN(sf), + GS_XRANGE(sf), GS_YMAXMIN(sf), GS_YRANGE(sf)) + default: + call error (0, "GSVECTOR: Unknown surface type.") + } +end -- cgit