Index: /trunk/psLib/src/math/psPolynomialMD.c
===================================================================
--- /trunk/psLib/src/math/psPolynomialMD.c	(revision 17074)
+++ /trunk/psLib/src/math/psPolynomialMD.c	(revision 17075)
@@ -3,13 +3,14 @@
 #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,5 +399,6 @@
             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,4 +409,7 @@
             }
             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,4 +417,5 @@
             }
         }
+        psTrace("psLib.math", 7, "Rejected %d points (%d remaining)\n", numClipped, numGood);
     }