Initial commit

1 files changed, 28 insertions, 0 deletions

diff --git a/noao/rv/rvidlines/idrms.x b/noao/rv/rvidlines/idrms.x
new file mode 100644
index 00000000..257675de
--- /dev/null
+++ b/noao/rv/rvidlines/idrms.x
@@ -0,0 +1,28 @@
+include	"identify.h"
+
+# ID_RMS -- Compute RMS of fit about the user coordinates
+
+double procedure id_rms (id)
+
+pointer	id			# ID pointer
+
+int	i, nrms
+double	rms, id_zshiftd()
+
+begin
+	rms = 0.
+	nrms = 0
+	for (i=1; i<=ID_NFEATURES(id); i=i+1) {
+	    if (!IS_INDEFD (USER(id,i)) && WTS(id,i) != 0.) {
+		rms = rms + (id_zshiftd (id, FIT(id,i), 0) - USER(id,i)) ** 2
+		nrms = nrms + 1
+	    }
+	}
+
+	if (nrms > 0)
+	    rms = sqrt (rms / nrms)
+	else
+	    rms = INDEFD
+
+	return (rms)
+end