Initial commit

1 files changed, 40 insertions, 0 deletions

diff --git a/math/gsurfit/gszeror.x b/math/gsurfit/gszeror.x
new file mode 100644
index 00000000..f7c4e5ed
--- /dev/null
+++ b/math/gsurfit/gszeror.x
@@ -0,0 +1,40 @@
+# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
+
+include <math/gsurfit.h>
+include "gsurfitdef.h"
+
+# GSZERO -- Procedure to zero the accumulators before doing
+# a new fit in accumulate mode. The inner products of the basis functions
+# are accumulated in the GS_NCOEFF(sf) ** 2
+# array MATRIX, while
+# the inner products of the basis functions and the data ordinates are
+# accumulated in the  NCOEFF(sf)-vector VECTOR.
+
+procedure gszero (sf)
+
+pointer	sf	# pointer to surface descriptor
+errchk	mfree
+
+begin
+	# zero the accumulators
+	switch (GS_TYPE(sf)) {
+	case GS_LEGENDRE, GS_CHEBYSHEV, GS_POLYNOMIAL:
+
+	    GS_NPTS(sf) = 0
+	    call aclrr (VECTOR(GS_VECTOR(sf)), GS_NCOEFF(sf))
+	    call aclrr (MATRIX(GS_MATRIX(sf)), GS_NCOEFF(sf) ** 2)
+
+	    # free the basis functions defined from previous calls to sfrefit
+	    if (GS_XBASIS(sf) != NULL)
+	        call mfree (GS_XBASIS(sf), TY_REAL)
+	    GS_XBASIS(sf) = NULL
+	    if (GS_YBASIS(sf) != NULL)
+	        call mfree (GS_YBASIS(sf), TY_REAL)
+	    GS_YBASIS(sf) = NULL
+	    if (GS_WZ(sf) != NULL)
+	        call mfree (GS_WZ(sf), TY_REAL)
+	    GS_WZ(sf) = NULL
+	default:
+	    call error (0, "GSZERO: Unknown surface type.")
+	}
+end