1 files changed, 83 insertions, 0 deletions
diff --git a/math/gsurfit/gsstatr.x b/math/gsurfit/gsstatr.x
new file mode 100644
index 00000000..826bcafa
--- /dev/null
+++ b/math/gsurfit/gsstatr.x
@@ -0,0 +1,83 @@
+# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
+
+include <math/gsurfit.h>
+include "gsurfitdef.h"
+
+# GSGET -- Procedure to fetch a gsurfit parameter
+real procedure gsgetr (sf, parameter)
+
+pointer	sf		# pointer to the surface fit
+int	parameter	# parameter to be fetched
+
+begin
+	switch (parameter) {
+	case GSXMAX:
+	    return (GS_XMAX(sf))
+	case GSXMIN:
+	    return (GS_XMIN(sf))
+	case GSYMAX:
+	    return (GS_YMAX(sf))
+	case GSYMIN:
+	    return (GS_YMIN(sf))
+	case GSXREF:
+	    return (GS_XREF(sf))
+	case GSYREF:
+	    return (GS_YREF(sf))
+	case GSZREF:
+	    return (GS_ZREF(sf))
+	}
+end
+
+
+# GSSET -- Procedure to set a gsurfit parameter
+procedure gsset (sf, parameter, val)
+
+pointer	sf		# pointer to the surface fit
+int	parameter	# parameter to be fetched
+real	val		# value to set
+
+begin
+	switch (parameter) {
+	case GSXREF:
+	    GS_XREF(sf) = val
+	case GSYREF:
+	    GS_YREF(sf) = val
+	case GSZREF:
+	    GS_ZREF(sf) = val
+	}
+end
+
+
+# GSGETI -- Procedure to fetch an integer parameter
+
+int procedure gsgeti (sf, parameter)
+
+pointer sf		# pointer to the surface fit
+int	parameter	# integer parameter
+
+begin
+	switch (parameter) {
+	case GSTYPE:
+	    return (GS_TYPE(sf))
+	case GSXORDER:
+	    switch (GS_TYPE(sf)) {
+	    case GS_LEGENDRE, GS_CHEBYSHEV, GS_POLYNOMIAL:
+		return (GS_XORDER(sf))
+	    }
+	case GSYORDER:
+	    switch (GS_TYPE(sf)) {
+	    case GS_LEGENDRE, GS_CHEBYSHEV, GS_POLYNOMIAL:
+	        return (GS_YORDER(sf))
+	    }
+	case GSXTERMS:
+	    return (GS_XTERMS(sf))
+	case GSNXCOEFF:
+	    return (GS_NXCOEFF(sf))
+	case GSNYCOEFF:
+	    return (GS_NYCOEFF(sf))
+	case GSNCOEFF:
+	    return (GS_NCOEFF(sf))
+	case GSNSAVE:
+	    return (GS_SAVECOEFF+GS_NCOEFF(sf))
+	}
+end