blob: 257675debf6ec5fe17d2f4df0289d1c627eac402 (
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, 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
|
