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
|