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
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
# procedure to evaluate interpolator at a sequence of ordered points
# aasumes that all x_values land within [1,n].
procedure asieva(x, y, n, coeff)
include "interpdef.h"
include "asidef.h"
real x[ARB] # ordered x array
int n # no. of points in x
real coeff[ARB]
real y[ARB] # interpolated values
begin
switch (ITYPEI) { # switch on interpolator type
case IT_NEAREST :
call iievne(x, y, n, coeff[COFF + 1])
case IT_LINEAR :
call iievli(x, y, n, coeff[COFF + 1])
case IT_POLY3 :
call iievp3(x, y, n, coeff[COFF + 1])
case IT_POLY5 :
call iievp5(x, y, n, coeff[COFF + 1])
case IT_SPLINE3 :
call iievs3(x, y, n, coeff[COFF + 1])
}
return
end
|
