From 40e5a5811c6ffce9b0974e93cdd927cbcf60c157 Mon Sep 17 00:00:00 2001
From: Joe Hunkeler <jhunkeler@gmail.com>
Date: Tue, 11 Aug 2015 16:51:37 -0400
Subject: Repatch (from linux) of OSX IRAF

---
 math/curfit/cvevald.x | 51 +++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 51 insertions(+)
 create mode 100644 math/curfit/cvevald.x

(limited to 'math/curfit/cvevald.x')

diff --git a/math/curfit/cvevald.x b/math/curfit/cvevald.x
new file mode 100644
index 00000000..c1c1f052
--- /dev/null
+++ b/math/curfit/cvevald.x
@@ -0,0 +1,51 @@
+# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
+
+include <math/curfit.h>
+
+include "dcurfitdef.h"
+
+# CVEVAL -- Procedure to evaluate curve at a given x. The CV_NCOEFF(cv)
+# coefficients are assumed to be in COEFF.
+
+double procedure dcveval (cv, x)
+
+pointer	cv		# curve descriptor
+double	x		# x value
+
+int	left
+pointer	cptr, xptr
+double	yfit
+
+double	adotd()
+
+begin
+
+	# calculate the non-zero basis functions
+	switch (CV_TYPE(cv)) {
+	case CHEBYSHEV:
+	    left = 0
+	    call dcv_b1cheb (x, CV_ORDER(cv), CV_MAXMIN(cv), CV_RANGE(cv),
+			   XBASIS(CV_XBASIS(cv)))
+	case LEGENDRE:
+	    left = 0
+	    call dcv_b1leg (x, CV_ORDER(cv), CV_MAXMIN(cv), CV_RANGE(cv),
+			   XBASIS(CV_XBASIS(cv)))
+	case SPLINE3:
+	    call dcv_b1spline3 (x, CV_NPIECES(cv), -CV_XMIN(cv),
+	    		      CV_SPACING(cv), XBASIS(CV_XBASIS(cv)), left)
+	case SPLINE1:
+	    call dcv_b1spline1 (x, CV_NPIECES(cv), -CV_XMIN(cv),
+	    		      CV_SPACING(cv), XBASIS(CV_XBASIS(cv)), left)
+	case USERFNC:
+	    left = 0
+	    call dcv_b1user (cv, x)
+	}
+
+
+	# accumulate the fitted value
+	cptr = CV_COEFF(cv) + left
+	xptr = CV_XBASIS(cv)
+	yfit = adotd (XBASIS(xptr), COEFF(cptr), CV_ORDER(cv))
+
+	return (yfit)
+end
-- 
cgit