1 files changed, 61 insertions, 0 deletions
diff --git a/math/gsurfit/gsurfitdef.h b/math/gsurfit/gsurfitdef.h
new file mode 100644
index 00000000..7ee6cc1d
--- /dev/null
+++ b/math/gsurfit/gsurfitdef.h
@@ -0,0 +1,61 @@
+# Header file for the surface fitting package
+
+# set up the curve descriptor structure
+
+define	LEN_GSSTRUCT	64
+
+define	GS_TYPE		Memi[$1]	# Type of curve to be fitted
+define	GS_XORDER	Memi[$1+1]	# Order of the fit in x
+define	GS_YORDER	Memi[$1+2]	# Order of the fit in y
+define	GS_XTERMS	Memi[$1+3]	# Cross terms for polynomials
+define	GS_NXCOEFF	Memi[$1+4]	# Number of x coefficients
+define	GS_NYCOEFF	Memi[$1+5]	# Number of y coefficients
+define	GS_NCOEFF	Memi[$1+6]	# Total number of coefficients
+define	GS_XREF		Memr[P2R($1+7)]	# x reference value
+define	GS_YREF		Memr[P2R($1+8)]	# y reference value
+define	GS_ZREF		Memr[P2R($1+9)]	# z reference value
+define	GS_XMAX		Memr[P2R($1+10)]# Maximum x value
+define	GS_XMIN		Memr[P2R($1+11)]# Minimum x value
+define	GS_YMAX		Memr[P2R($1+12)]# Maximum y value
+define	GS_YMIN		Memr[P2R($1+13)]# Minimum y value		
+define	GS_XRANGE	Memr[P2R($1+14)]# 2. / (xmax - xmin), polynomials
+define	GS_XMAXMIN	Memr[P2R($1+15)]# - (xmax + xmin) / 2., polynomials
+define	GS_YRANGE	Memr[P2R($1+16)]# 2. / (ymax - ymin), polynomials
+define	GS_YMAXMIN	Memr[P2R($1+17)]# - (ymax + ymin) / 2., polynomials
+define	GS_NPTS		Memi[$1+18]	# Number of data points
+
+define	GS_MATRIX	Memi[$1+19]	# Pointer to original matrix
+define	GS_CHOFAC	Memi[$1+20]	# Pointer to Cholesky factorization
+define	GS_VECTOR	Memi[$1+21]	# Pointer to  vector
+define	GS_COEFF	Memi[$1+22]	# Pointer to coefficient vector
+define	GS_XBASIS	Memi[$1+23]	# Pointer to basis functions (all x)
+define	GS_YBASIS	Memi[$1+24]	# Pointer to basis functions (all y)
+define	GS_WZ		Memi[$1+25]	# Pointer to w * z (gsrefit)
+
+# matrix and vector element definitions
+
+define	XBASIS		Memr[P2P($1)]	# Non zero basis for all x
+define	YBASIS		Memr[P2P($1)]	# Non zero basis for all y
+define	XBS		Memr[P2P($1)]	# Non zero basis for single x
+define	YBS		Memr[P2P($1)]	# Non zero basis for single y
+define	MATRIX		Memr[P2P($1)]	# Element of MATRIX
+define	CHOFAC		Memr[P2P($1)]	# Element of CHOFAC
+define	VECTOR		Memr[P2P($1)]	# Element of VECTOR
+define	COEFF		Memr[P2P($1)]	# Element of COEFF
+
+# structure definitions for restore
+
+define	GS_SAVETYPE	$1[1]
+define	GS_SAVEXORDER	$1[2]
+define	GS_SAVEYORDER	$1[3]
+define	GS_SAVEXTERMS	$1[4]
+define	GS_SAVEXMIN	$1[5]
+define	GS_SAVEXMAX	$1[6]
+define	GS_SAVEYMIN	$1[7]
+define	GS_SAVEYMAX	$1[8]
+
+# data type
+
+define	DELTA		EPSILON
+
+# miscellaneous