Changeset 17031

Timestamp:

Mar 17, 2008, 4:07:18 PM (18 years ago)

Author:

Paul Price

Message:

Normalisation wasn't set correctly for the case where we're not convolving.

File:

• trunk/ppStack/src/ppStackMatch.c (modified) (1 diff)

trunk/ppStack/src/ppStackMatch.c

@@ -133 +133 @@
         kernels->solution1 = psVectorAlloc(kernels->num + 2, PS_TYPE_F64);
         psVectorInit(kernels->solution1, 0.0);
-        for (int i = 0; i < kernels->num; i++) {
-            if (kernels->u->data.S32[i] == 0 && kernels->v->data.S32[i] == 0) {
-                kernels->solution1->data.F64[i] = 1.0;
-            }
-        }
+        int normIndex = PM_SUBTRACTION_INDEX_NORM(kernels); // Index for normalisation
+        kernels->solution1->data.F64[normIndex] = 1.0;
         psMetadataAddPtr(output->analysis, PS_LIST_TAIL, PM_SUBTRACTION_ANALYSIS_KERNEL,
                          PS_DATA_UNKNOWN | PS_META_DUPLICATE_OK, "Fake subtraction kernel", kernels);