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
|
include "hdicfit.h"
# HDIC_ADDPOINT -- Add a density, exposure, weights point into the sample.
procedure hdic_addpoint (ic,
nden, nexp, nwts, den, y, wts, userwts, x, wd, sdev, npts)
pointer ic # Pointer to ic data structure
real nden # New density value to be added
real nexp # New exposure value to be added
double nwts # New weight value to be added
pointer den # Pointer to existing density array
pointer y # Pointer to existing exposure array
pointer wts # Pointer to existing wts array
pointer userwts # Pointer to existing userwts array
pointer x # Pointer to existing array of ind variables
pointer wd # Pointer to flag array for deletion reasons
pointer sdev # Pointer to standard deviation array
int npts # Number of points, incremented on output
begin
npts = npts + 1
call realloc (den, npts, TY_DOUBLE)
call realloc (y, npts, TY_DOUBLE)
call realloc (wts, npts, TY_DOUBLE)
call realloc (userwts, npts, TY_DOUBLE)
call realloc (x, npts, TY_DOUBLE)
call realloc (wd, npts, TY_INT)
call realloc (sdev, npts, TY_DOUBLE)
Memd[den+npts-1] = double (nden)
Memd[y +npts-1] = double (nexp)
Memd[wts+npts-1] = double (nwts)
Memd[userwts+npts-1] = double (nwts)
Memi[wd +npts-1] = NDELETE
Memd[sdev+npts-1] = ADDED_PT
# Sort the data and then update the reference vector.
call hdic_sort (Memd[den], Memd[y], Memd[wts], Memd[userwts],
Memi[wd], Memd[sdev], npts)
call hdic_init (Memd[den], npts, Memd[den+npts-1])
IC_NEWX(ic) = YES
IC_NEWY(ic) = YES
IC_NEWWTS(ic) = YES
end
|
