blob: de84ae2612f89b4888384abcf559a6c2d523ad52 (
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 "ecidentify.h"
# EC_RMS -- Compute RMS of fit about the user coordinates
double procedure ec_rms (ec)
pointer ec # ID pointer
int i, nrms
double rms
begin
rms = 0.
nrms = 0
for (i=1; i<=EC_NFEATURES(ec); i=i+1) {
if (!IS_INDEFD (USER(ec,i)) && FTYPE(ec,i) > 0) {
rms = rms + (FIT(ec,i) - USER(ec,i)) ** 2
nrms = nrms + 1
}
}
if (nrms > 0)
rms = sqrt (rms / nrms)
else
rms = INDEFD
return (rms)
end
|
