diff options
| author | Joseph Hunkeler <jhunkeler@gmail.com> | 2015-07-08 20:46:52 -0400 |
|---|---|---|
| committer | Joseph Hunkeler <jhunkeler@gmail.com> | 2015-07-08 20:46:52 -0400 |
| commit | fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 (patch) | |
| tree | bdda434976bc09c864f2e4fa6f16ba1952b1e555 /math/gsurfit/gsvectorr.x | |
| download | iraf-linux-fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4.tar.gz | |
Initial commit
Diffstat (limited to 'math/gsurfit/gsvectorr.x')
| -rw-r--r-- | math/gsurfit/gsvectorr.x | 57 |
1 files changed, 57 insertions, 0 deletions
diff --git a/math/gsurfit/gsvectorr.x b/math/gsurfit/gsvectorr.x new file mode 100644 index 00000000..38213ecc --- /dev/null +++ b/math/gsurfit/gsvectorr.x @@ -0,0 +1,57 @@ +# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc. + +include <math/gsurfit.h> +include "gsurfitdef.h" + +# GSVECTOR -- Procedure to evaluate the fitted surface at an array of points. +# The GS_NCOEFF(sf) coefficients are stored in the +# vector COEFF. + +procedure gsvector (sf, x, y, zfit, npts) + +pointer sf # pointer to surface descriptor structure +real x[ARB] # x value +real y[ARB] # y value +real 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 rgs_1devpoly (COEFF(GS_COEFF(sf)), y, zfit, npts, + GS_YORDER(sf), GS_YMAXMIN(sf), GS_YRANGE(sf)) + } else if (GS_YORDER(sf) == 1) { + call rgs_1devpoly (COEFF(GS_COEFF(sf)), x, zfit, npts, + GS_XORDER(sf), GS_XMAXMIN(sf), GS_XRANGE(sf)) + } else + call rgs_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 rgs_1devcheb (COEFF(GS_COEFF(sf)), y, zfit, npts, + GS_YORDER(sf), GS_YMAXMIN(sf), GS_YRANGE(sf)) + } else if (GS_YORDER(sf) == 1) { + call rgs_1devcheb (COEFF(GS_COEFF(sf)), x, zfit, npts, + GS_XORDER(sf), GS_XMAXMIN(sf), GS_XRANGE(sf)) + } else + call rgs_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 rgs_1devleg (COEFF(GS_COEFF(sf)), y, zfit, npts, + GS_YORDER(sf), GS_YMAXMIN(sf), GS_YRANGE(sf)) + } else if (GS_YORDER(sf) == 1) { + call rgs_1devleg (COEFF(GS_COEFF(sf)), x, zfit, npts, + GS_XORDER(sf), GS_XMAXMIN(sf), GS_XRANGE(sf)) + } else + call rgs_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 |