diff options
| author | Joe Hunkeler <jhunkeler@gmail.com> | 2015-08-11 16:51:37 -0400 |
|---|---|---|
| committer | Joe Hunkeler <jhunkeler@gmail.com> | 2015-08-11 16:51:37 -0400 |
| commit | 40e5a5811c6ffce9b0974e93cdd927cbcf60c157 (patch) | |
| tree | 4464880c571602d54f6ae114729bf62a89518057 /math/gsurfit/gszero.gx | |
| download | iraf-osx-40e5a5811c6ffce9b0974e93cdd927cbcf60c157.tar.gz | |
Repatch (from linux) of OSX IRAF
Diffstat (limited to 'math/gsurfit/gszero.gx')
| -rw-r--r-- | math/gsurfit/gszero.gx | 60 |
1 files changed, 60 insertions, 0 deletions
diff --git a/math/gsurfit/gszero.gx b/math/gsurfit/gszero.gx new file mode 100644 index 00000000..e99cbe4d --- /dev/null +++ b/math/gsurfit/gszero.gx @@ -0,0 +1,60 @@ +# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc. + +include <math/gsurfit.h> +$if (datatype == r) +include "gsurfitdef.h" +$else +include "dgsurfitdef.h" +$endif + +# GSZERO -- Procedure to zero the accumulators before doing +# a new fit in accumulate mode. The inner products of the basis functions +# are accumulated in the GS_NCOEFF(sf) ** 2 +# array MATRIX, while +# the inner products of the basis functions and the data ordinates are +# accumulated in the NCOEFF(sf)-vector VECTOR. + +$if (datatype == r) +procedure gszero (sf) +$else +procedure dgszero (sf) +$endif + +pointer sf # pointer to surface descriptor +errchk mfree + +begin + # zero the accumulators + switch (GS_TYPE(sf)) { + case GS_LEGENDRE, GS_CHEBYSHEV, GS_POLYNOMIAL: + + GS_NPTS(sf) = 0 + call aclr$t (VECTOR(GS_VECTOR(sf)), GS_NCOEFF(sf)) + call aclr$t (MATRIX(GS_MATRIX(sf)), GS_NCOEFF(sf) ** 2) + + # free the basis functions defined from previous calls to sfrefit + $if (datatype == r) + if (GS_XBASIS(sf) != NULL) + call mfree (GS_XBASIS(sf), TY_REAL) + GS_XBASIS(sf) = NULL + if (GS_YBASIS(sf) != NULL) + call mfree (GS_YBASIS(sf), TY_REAL) + GS_YBASIS(sf) = NULL + if (GS_WZ(sf) != NULL) + call mfree (GS_WZ(sf), TY_REAL) + GS_WZ(sf) = NULL + $else + if (GS_XBASIS(sf) != NULL) + call mfree (GS_XBASIS(sf), TY_DOUBLE) + GS_XBASIS(sf) = NULL + if (GS_YBASIS(sf) != NULL) + call mfree (GS_YBASIS(sf), TY_DOUBLE) + GS_YBASIS(sf) = NULL + if (GS_WZ(sf) != NULL) + call mfree (GS_WZ(sf), TY_DOUBLE) + GS_WZ(sf) = NULL + $endif + default: + call error (0, "GSZERO: Unknown surface type.") + } +end |