Changeset 17075 for trunk/psLib/src/math/psPolynomialMD.c

Timestamp:

Mar 20, 2008, 2:37:35 PM (18 years ago)

Author:

Paul Price

Message:

Some idiot forgot the definition of "RMS", resulting in many values being excluded.

File:

• trunk/psLib/src/math/psPolynomialMD.c (modified) (4 diffs)

trunk/psLib/src/math/psPolynomialMD.c

@@ -3 +3 @@
 #include <limits.h>
 
-#include <psMemory.h>
-#include <psVector.h>
-#include <psImage.h>
-#include <psPolynomial.h>
-#include <psArray.h>
-#include <psMatrix.h>
-#include <psError.h>
-#include <psAbort.h>
-#include <psAssert.h>
+#include "psMemory.h"
+#include "psVector.h"
+#include "psImage.h"
+#include "psPolynomial.h"
+#include "psArray.h"
+#include "psMatrix.h"
+#include "psError.h"
+#include "psAbort.h"
+#include "psAssert.h"
+#include "psTrace.h"
 
 #include "psPolynomialMD.h"
@@ -398 +399 @@
             rms += PS_SQR(diff);
         }
-        rms = sqrt(rms) / (double)numGood;
+        rms = sqrt(rms / (double)numGood);
+        psTrace("psLib.math", 7, "RMS from %d points is %lf\n", numGood, rms);
 
         // Reject
@@ -407 +409 @@
             }
             if (fabs(deviations->data.F32[i]) > limit) {
+                psTrace("psLib.math", 9, "Rejected point %d (%f,%f: %f > %f)\n",
+                        i, values->data.F32[i], psPolynomialMDEval(poly, coordsArray->data[i]),
+                        deviations->data.F32[i], limit);
                 ownMask->data.U8[i] = 0xff;
                 numClipped++;
@@ -412 +417 @@
             }
         }
+        psTrace("psLib.math", 7, "Rejected %d points (%d remaining)\n", numClipped, numGood);
     }