Initial commit

1 files changed, 47 insertions, 0 deletions

diff --git a/noao/imred/dtoi/hdicfit/hdicadd.x b/noao/imred/dtoi/hdicfit/hdicadd.x
new file mode 100644
index 00000000..1fae152d
--- /dev/null
+++ b/noao/imred/dtoi/hdicfit/hdicadd.x
@@ -0,0 +1,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