Changeset 26547 for branches/eam_branches/20091201/psModules/src/imcombine/pmSubtraction.c

Timestamp:

Jan 8, 2010, 2:24:01 PM (17 years ago)

Author:

Paul Price

Message:

Dual convolution seems to be working again!

File:

• branches/eam_branches/20091201/psModules/src/imcombine/pmSubtraction.c (modified) (1 diff)

branches/eam_branches/20091201/psModules/src/imcombine/pmSubtraction.c

@@ -117 +117 @@
     for (int i = 0; i < numKernels; i++) {
         double value = p_pmSubtractionSolutionCoeff(kernels, polyValues, i, wantDual); // Polynomial value
+        if (wantDual) {
+            // The model is built with the dual convolution terms added, so to produce zero residual the
+            // equation results in negative coefficients which we must undo.
+            value *= -1.0;
+        }
 
         switch (kernels->type) {