diff options
Diffstat (limited to 'math/curfit/cvevald.x')
| -rw-r--r-- | math/curfit/cvevald.x | 51 |
1 files changed, 51 insertions, 0 deletions
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 |