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/gsscoeff.gx | |
| download | iraf-osx-40e5a5811c6ffce9b0974e93cdd927cbcf60c157.tar.gz | |
Repatch (from linux) of OSX IRAF
Diffstat (limited to 'math/gsurfit/gsscoeff.gx')
| -rw-r--r-- | math/gsurfit/gsscoeff.gx | 54 |
1 files changed, 54 insertions, 0 deletions
diff --git a/math/gsurfit/gsscoeff.gx b/math/gsurfit/gsscoeff.gx new file mode 100644 index 00000000..09b894e3 --- /dev/null +++ b/math/gsurfit/gsscoeff.gx @@ -0,0 +1,54 @@ +# 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 + +# GSSCOEFF -- Procedure to set a particular coefficient. +# If the requested coefficient is undefined then the coefficient is not set. + +$if (datatype == r) +procedure gsscoeff (sf, xorder, yorder, coeff) +$else +procedure dgsscoeff (sf, xorder, yorder, coeff) +$endif + +pointer sf # pointer to the surface fitting descriptor +int xorder # X order of desired coefficent +int yorder # Y order of desired coefficent +PIXEL coeff # Coefficient value + +int i, n, maxorder, xincr + +begin + if ((xorder > GS_XORDER(sf)) || (yorder > GS_YORDER(sf))) + return + + switch (GS_XTERMS(sf)) { + case GS_XNONE: + if (yorder == 1) + n = xorder + else if (xorder == 1) + n = GS_NXCOEFF(sf) + yorder - 1 + else + return + case GS_XHALF: + maxorder = max (GS_XORDER(sf) + 1, GS_YORDER(sf) + 1) + if ((xorder + yorder) > maxorder) + return + n = xorder + xincr = GS_XORDER(sf) + do i = 2, yorder { + n = n + xincr + if ((i + GS_XORDER(sf) + 1) > maxorder) + xincr = xincr - 1 + } + case GS_XFULL: + n = xorder + (yorder - 1) * GS_NXCOEFF(sf) + } + + COEFF(GS_COEFF(sf) + n - 1) = coeff +end |