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/gsgcoeff.gx | |
| download | iraf-osx-40e5a5811c6ffce9b0974e93cdd927cbcf60c157.tar.gz | |
Repatch (from linux) of OSX IRAF
Diffstat (limited to 'math/gsurfit/gsgcoeff.gx')
| -rw-r--r-- | math/gsurfit/gsgcoeff.gx | 53 |
1 files changed, 53 insertions, 0 deletions
diff --git a/math/gsurfit/gsgcoeff.gx b/math/gsurfit/gsgcoeff.gx new file mode 100644 index 00000000..3d8a294d --- /dev/null +++ b/math/gsurfit/gsgcoeff.gx @@ -0,0 +1,53 @@ +# 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 + +# GSGCOEFF -- Procedure to fetch a particular coefficient. +# If the requested coefficient is undefined then INDEF is returned. + +$if (datatype == r) +real procedure gsgcoeff (sf, xorder, yorder) +$else +double procedure dgsgcoeff (sf, xorder, yorder) +$endif + +pointer sf # pointer to the surface fitting descriptor +int xorder # X order of desired coefficent +int yorder # Y order of desired coefficent + +int i, n, maxorder, xincr + +begin + if ((xorder > GS_XORDER(sf)) || (yorder > GS_YORDER(sf))) + return (INDEF) + + switch (GS_XTERMS(sf)) { + case GS_XNONE: + if (yorder == 1) + n = xorder + else if (xorder == 1) + n = GS_NXCOEFF(sf) + yorder - 1 + else + return (INDEF) + case GS_XHALF: + maxorder = max (GS_XORDER(sf) + 1, GS_YORDER(sf) + 1) + if ((xorder + yorder) > maxorder) + return (INDEF) + 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) + } + + return (COEFF(GS_COEFF(sf) + n - 1)) +end |