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
|
include <pkg/inlfit.h>
# ING_RESULTS -- Print the results of the fit.
procedure ing_results$t (in, file, nl, x, y, wts, names, npts, nvars, len_name)
pointer in # INLFIT pointer
char file[ARB] # Output file name
pointer nl # NLFIT pointer
PIXEL x[ARB] # Ordinates (npts * nvars)
PIXEL y[ARB] # Abscissas
PIXEL wts[ARB] # Weights
char names[ARB] # Object names
int npts # Number of data points
int nvars # Number of variables
int len_name # Length of a name
int i, fd, rejected
pointer sp, fit, wts1, rejpts
int open(), in_geti()
pointer in_getp()
errchk open
begin
# Open the output file.
if (file[1] == EOS)
return
fd = open (file, APPEND, TEXT_FILE)
# Test the number of points.
if (npts == 0) {
call eprintf ("Incomplete output - no data points for fit\n")
return
}
# Allocate memory.
call smark (sp)
call salloc (fit, npts, TY_PIXEL)
call salloc (wts1, npts, TY_PIXEL)
# Evaluate the fit.
call nlvector$t (nl, x, Mem$t[fit], npts, nvars)
# Assign a zero weight to the rejected points.
rejected = in_geti (in, INLNREJPTS)
rejpts = in_getp (in, INLREJPTS)
call amov$t (wts, Mem$t[wts1], npts)
if (rejected > 0) {
do i = 1, npts {
if (Memi[rejpts+i-1] == YES)
Mem$t[wts1+i-1] = PIXEL (0.0)
}
}
# Print the title.
call fprintf (fd, "\n#%14.14s %14.14s %14.14s")
call pargstr ("objectid")
call pargstr ("function")
call pargstr ("fit")
call fprintf (fd, " %14.14s %14.14s\n")
call pargstr ("residuals")
call pargstr ("sigma")
# List function value, fit value, residual and error values.
do i = 1, npts {
call fprintf (fd, " %14.14s %14.7g %14.7g %14.7g %14.7g\n")
call pargstr (names[(i-1)*len_name+1])
if (Mem$t[wts1+i-1] <= 0.0) {
call parg$t (INDEF)
call parg$t (INDEF)
call parg$t (INDEF)
call parg$t (INDEF)
} else {
call parg$t (y[i])
call parg$t (Mem$t[fit+i-1])
call parg$t (y[i] - Mem$t[fit+i-1])
call parg$t (sqrt (PIXEL (1.0) / Mem$t[wts1+i-1]))
}
}
call fprintf (fd, "\n")
# Free allocated memory, and close output file.
call sfree (sp)
call close (fd)
end
|
