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
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <math/gsurfit.h>
include "dgsurfitdef.h"
# GSVECTOR -- Procedure to evaluate the fitted surface at an array of points.
# The GS_NCOEFF(sf) coefficients are stored in the
# vector COEFF.
procedure dgsvector (sf, x, y, zfit, npts)
pointer sf # pointer to surface descriptor structure
double x[ARB] # x value
double y[ARB] # y value
double zfit[ARB] # fits surface values
int npts # number of data points
begin
# evaluate the surface along the vector
switch (GS_TYPE(sf)) {
case GS_POLYNOMIAL:
if (GS_XORDER(sf) == 1) {
call dgs_1devpoly (COEFF(GS_COEFF(sf)), y, zfit, npts,
GS_YORDER(sf), GS_YMAXMIN(sf), GS_YRANGE(sf))
} else if (GS_YORDER(sf) == 1) {
call dgs_1devpoly (COEFF(GS_COEFF(sf)), x, zfit, npts,
GS_XORDER(sf), GS_XMAXMIN(sf), GS_XRANGE(sf))
} else
call dgs_evpoly (COEFF(GS_COEFF(sf)), x, y, zfit, npts,
GS_XTERMS(sf), GS_XORDER(sf), GS_YORDER(sf), GS_XMAXMIN(sf),
GS_XRANGE(sf), GS_YMAXMIN(sf), GS_YRANGE(sf))
case GS_CHEBYSHEV:
if (GS_XORDER(sf) == 1) {
call dgs_1devcheb (COEFF(GS_COEFF(sf)), y, zfit, npts,
GS_YORDER(sf), GS_YMAXMIN(sf), GS_YRANGE(sf))
} else if (GS_YORDER(sf) == 1) {
call dgs_1devcheb (COEFF(GS_COEFF(sf)), x, zfit, npts,
GS_XORDER(sf), GS_XMAXMIN(sf), GS_XRANGE(sf))
} else
call dgs_evcheb (COEFF(GS_COEFF(sf)), x, y, zfit, npts,
GS_XTERMS(sf), GS_XORDER(sf), GS_YORDER(sf), GS_XMAXMIN(sf),
GS_XRANGE(sf), GS_YMAXMIN(sf), GS_YRANGE(sf))
case GS_LEGENDRE:
if (GS_XORDER(sf) == 1) {
call dgs_1devleg (COEFF(GS_COEFF(sf)), y, zfit, npts,
GS_YORDER(sf), GS_YMAXMIN(sf), GS_YRANGE(sf))
} else if (GS_YORDER(sf) == 1) {
call dgs_1devleg (COEFF(GS_COEFF(sf)), x, zfit, npts,
GS_XORDER(sf), GS_XMAXMIN(sf), GS_XRANGE(sf))
} else
call dgs_evleg (COEFF(GS_COEFF(sf)), x, y, zfit, npts,
GS_XTERMS(sf), GS_XORDER(sf), GS_YORDER(sf), GS_XMAXMIN(sf),
GS_XRANGE(sf), GS_YMAXMIN(sf), GS_YRANGE(sf))
default:
call error (0, "GSVECTOR: Unknown surface type.")
}
end
|
