1 /* calcsd.c */
2 #include "ransurf.h"
3 #include "local_proto.h"
4
CalcSD(void)5 void CalcSD(void)
6 {
7 int Row, Col, DoFilter;
8 double Effect;
9
10 G_debug(2, "CalcSD()");
11
12 FilterSD = 0.0;
13 for (DoFilter = 0; DoFilter < NumFilters; DoFilter++) {
14 CopyFilter(&Filter, AllFilters[DoFilter]);
15 if (Filter.Mult < 0.0)
16 Filter.Mult *= -1.0;
17
18 MakeBigF();
19 for (Row = 0; Row < BigF.NumR; Row++) {
20 for (Col = 0; Col < BigF.NumC; Col++) {
21 DistDecay(&Effect, Row - BigF.RowPlus, Col - BigF.ColPlus);
22 FilterSD += Effect * Effect;
23 }
24 }
25 }
26
27 G_debug(3, "(FilterSD):%.12lf", FilterSD);
28 FilterSD = sqrt(FilterSD);
29 G_debug(3, "(FilterSD):%.12lf", FilterSD);
30
31 return;
32 }
33