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
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
|
include <pkg/gtools.h>
include <pkg/inlfit.h>
# ING_AXES -- Set axes data. The applications program may set additional
# axes types.
procedure ing_axesd (in, gt, nl, axis, x, y, z, npts, nvars)
pointer in # INLFIT pointer
pointer gt # GTOOLS pointer
pointer nl # NLFIT pointer
int axis # Output axis
double x[ARB] # Independent variables (npts * nvars)
double y[npts] # Dependent variable
double z[npts] # Output values
int npts # Number of points
int nvars # Number of variables
int i, j
int axistype, axisnum
int gtlabel[2], gtunits[2]
double a, b, xmin, xmax
pointer sp, label, units, minptr, maxptr
double nlevald()
double ing_dvzd()
errchk adivd()
extern ing_dvzd()
data gtlabel/GTXLABEL, GTYLABEL/
data gtunits/GTXUNITS, GTYUNITS/
int in_geti()
pointer in_getp()
begin
# Allocate string space.
call smark (sp)
call salloc (label, SZ_LINE, TY_CHAR)
call salloc (units, SZ_LINE, TY_CHAR)
# Get the appropiate axis type and variable number.
call in_gkey (in, in_geti (in, INLGKEY), axis, axistype, axisnum)
# Get and set axes labels and units.
call ing_getlabel (in, axistype, axisnum, Memc[label], Memc[units],
SZ_LINE)
call gt_sets (gt, gtlabel[axis], Memc[label])
call gt_sets (gt, gtunits[axis], Memc[units])
# Branch on axis type.
switch (axistype) {
case KEY_VARIABLE: # Independent variable
do i = 1, npts
z[i] = x[(i-1)*nvars+axisnum]
case KEY_FUNCTION: # Function variable
call amovd (y, z, npts)
case KEY_FIT: # Fitted values
call nlvectord (nl, x, z, npts, nvars)
case KEY_RESIDUALS: # Residuals
call nlvectord (nl, x, z, npts, nvars)
call asubd (y, z, z, npts)
case KEY_RATIO: # Ratio
call nlvectord (nl, x, z, npts, nvars)
call advzd (y, z, z, npts, ing_dvzd)
case KEY_NONLINEAR: # Linear component removed
call aclrd (z, npts)
minptr = in_getp (in, INLXMIN)
maxptr = in_getp (in, INLXMAX)
a = nlevald (nl, Memd[minptr], nvars)
do i = 1, nvars {
xmin = Memd[minptr+i-1]
xmax = Memd[maxptr+i-1]
Memd[minptr+i-1] = xmax
b = (nlevald (nl, Memd[minptr], nvars) - a) /
(xmax - xmin)
Memd[minptr+i-1] = xmin
do j = 1, npts
z[j] = z[j] + y[j] - a - b * (x[(j-1)*nvars+i] - xmin)
}
case KEY_UAXIS: # User axes plots.
if (axis == 1) {
do i = 1, npts
z[i] = x[(i-1)*nvars+axisnum]
} else
call amovd (y, z, npts)
call ing_uaxesd (axisnum, in, nl, x, y, z, npts, nvars)
default:
call error (0, "ing_axes: Unknown axis type")
}
# Free memory.
call sfree (sp)
end
# ING_DVZ -- Error action to take on zero division.
double procedure ing_dvzd (x)
double x # Numerator
begin
return (double (1.0))
end
|
