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
|
include <error.h>
include "../lib/parser.h"
# FT_RFEVAL - Evaluate reference equation (left hand side) for all observations
# and store them into a one column table, called the reference table (rtable).
procedure ft_rfeval (code, otable, rtable)
pointer code # equation code
pointer otable # 2d observation table
pointer rtable # 1d reference table (output)
int n
real rval
real dummy[1]
#bool clgetb()
int mct_nrows()
real pr_eval()
pointer mct_getrow()
begin
# Debug ?
#if (clgetb ("debug.fitcode")) {
#call eprintf ("ft_rfeval (code=%d) (otable=%d)\n")
#call pargi (code)
#call pargi (otable)
#}
# Allocate space for reference table.
# call mct_alloc (rtable, mct_nrows (otable), 1, TY_REAL)
# Loop over all data in the table.
do n = 1, mct_nrows (otable) {
# Evaluate the equation.
iferr (rval = pr_eval (code, Memr[mct_getrow (otable, n)], dummy)) {
call eprintf ("ft_ref (%d)\n")
call pargi (n)
call erract (EA_ERROR)
}
# Put data into reference table.
call mct_putr (rtable, n, 1, rval)
}
# Debug ?
#call dg_dref ("from ft_rfeval", rtable)
end
|
