diff options
Diffstat (limited to 'math/gsurfit/gsscoeffr.x')
| -rw-r--r-- | math/gsurfit/gsscoeffr.x | 46 |
1 files changed, 46 insertions, 0 deletions
diff --git a/math/gsurfit/gsscoeffr.x b/math/gsurfit/gsscoeffr.x new file mode 100644 index 00000000..dacb4bd0 --- /dev/null +++ b/math/gsurfit/gsscoeffr.x @@ -0,0 +1,46 @@ +# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc. + +include <math/gsurfit.h> +include "gsurfitdef.h" + +# GSSCOEFF -- Procedure to set a particular coefficient. +# If the requested coefficient is undefined then the coefficient is not set. + +procedure gsscoeff (sf, xorder, yorder, coeff) + +pointer sf # pointer to the surface fitting descriptor +int xorder # X order of desired coefficent +int yorder # Y order of desired coefficent +real 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 |