Repatch (from linux) of OSX IRAF

1 files changed, 83 insertions, 0 deletions

diff --git a/math/gsurfit/gsstatd.x b/math/gsurfit/gsstatd.x
new file mode 100644
index 00000000..b8c551f1
--- /dev/null
+++ b/math/gsurfit/gsstatd.x
@@ -0,0 +1,83 @@
+# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
+
+include <math/gsurfit.h>
+include "dgsurfitdef.h"
+
+# GSGET -- Procedure to fetch a gsurfit parameter
+double procedure dgsgetd (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 dgsset (sf, parameter, val)
+
+pointer	sf		# pointer to the surface fit
+int	parameter	# parameter to be fetched
+double	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 dgsgeti (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