aboutsummaryrefslogtreecommitdiff
path: root/noao/onedspec/identify/idrms.x
blob: 82916f1ad08a0fce49413d9e5e9ec68edf95c8be (plain) (blame)
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
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

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 + (FIT(id,i) - USER(id,i)) ** 2
		nrms = nrms + 1
	    }
	}

	if (nrms > 0)
	    rms = sqrt (rms / nrms)
	else
	    rms = INDEFD

	return (rms)
end