1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <math/curfit.h>
include "curfitdef.h"
# CVSET -- Procedure to store the fit parameters derived from outside
# the CURFIT package inside the curve descriptor structure for use
# by the CVEVAL and CVVECTOR proocedures. The curve_type is one of
# LEGENDRE, CHEBYSHEV or SPLINE3. For the polynomials the number of
# coefficients is equal to one plus the order of the polynomial. In the
# case of the cubic spline the number of coefficients equals three plus
# the number of polynomial pieces. The polynomials are normalized over
# from xmin to xmax.
procedure cvset (cv, curve_type, xmin, xmax, coeff, ncoeff)
pointer cv # curve descriptor
int curve_type # the functional form of the curve
real xmin # the minimum x value
real xmax # the maximum x value
real coeff[ncoeff] # the coefficient array
int ncoeff # the number of coefficients
errchk malloc
begin
# allocate space for curve descriptor
call malloc (cv, LEN_CVSTRUCT, TY_STRUCT)
if (ncoeff < 1)
call error (0, "CVSET: Illegal number of coefficients.")
if (xmin >= xmax)
call error (0, "CVSET: xmax <= xmin.")
# set curve_type dependent curve descriptor parameters
switch (curve_type) {
case CHEBYSHEV, LEGENDRE:
CV_ORDER(cv) = ncoeff
CV_NCOEFF(cv) = ncoeff
CV_RANGE(cv) = 2. / (xmax - xmin)
CV_MAXMIN(cv) = - (xmax + xmin) / 2.
case SPLINE3:
CV_ORDER(cv) = SPLINE3_ORDER
CV_NCOEFF(cv) = ncoeff
CV_NPIECES(cv) = ncoeff - SPLINE3_ORDER
CV_SPACING(cv) = (CV_NPIECES(cv) + 1) / (xmax - xmin)
case SPLINE1:
CV_ORDER(cv) = SPLINE1_ORDER
CV_NCOEFF(cv) = ncoeff
CV_NPIECES(cv) = ncoeff - SPLINE1_ORDER
CV_SPACING(cv) = (CV_NPIECES(cv) + 1) / (xmax - xmin)
case USERFNC:
CV_ORDER(cv) = ncoeff
CV_NCOEFF(cv) = ncoeff
CV_RANGE(cv) = 2. / (xmax - xmin)
CV_MAXMIN(cv) = - (xmax + xmin) / 2.
default:
call error (0, "CVSET: Unknown curve type.")
}
# set remaining curve parameters
CV_TYPE(cv) = curve_type
CV_XMIN(cv) = xmin
CV_XMAX(cv) = xmax
# allocate space for xbasis and coefficient arrays, set remaining
# pointers to NULL
call malloc (CV_XBASIS(cv), CV_ORDER(cv), TY_REAL)
call malloc (CV_COEFF(cv), CV_NCOEFF(cv), TY_REAL)
CV_MATRIX(cv) = NULL
CV_CHOFAC(cv) = NULL
CV_VECTOR(cv) = NULL
CV_BASIS(cv) = NULL
CV_WY(cv) = NULL
CV_LEFT(cv) = NULL
CV_USERFNC(cv) = NULL
# restore coefficients
call amovr (coeff, COEFF(CV_COEFF(cv)), CV_NCOEFF(cv))
end
|
