Changeset 14540 for trunk/psModules/src/imcombine/pmSubtractionKernels.c

Timestamp:

Aug 16, 2007, 3:50:50 PM (19 years ago)

Author:

Paul Price

Message:

Moving photometric scaling into the kernel for RINGS, and fixing up normalisation for odd functions.

File:

• trunk/psModules/src/imcombine/pmSubtractionKernels.c (modified) (2 diffs)

trunk/psModules/src/imcombine/pmSubtractionKernels.c

@@ -526 +526 @@
     for (int i = 0, index = 0; i < numRings + 1; i++) {
 
-        float lower2;           // Lower limit of radius^2
-        float upper2;           // Upper limit of radius^2
+        float lower2;                   // Lower limit of radius^2
+        float upper2;                   // Upper limit of radius^2
         if (i == 0) {
             lower2 = 0;
@@ -599 +599 @@
                     if (uOrder % 2 == 0 && vOrder % 2 == 0) {
                         psBinaryOp(poly, poly, "*", psScalarAlloc(1.0 / norm, PS_TYPE_F32));
+                        uCoords->data.S32[j] = 0;
+                        vCoords->data.S32[j] = 0;
+                        poly->data.F32[j] = -1.0;
+                        psVectorExtend(uCoords, RINGS_BUFFER, 1);
+                        psVectorExtend(vCoords, RINGS_BUFFER, 1);
+                        psVectorExtend(poly, RINGS_BUFFER, 1);
+                    } else {
+                        psBinaryOp(poly, poly, "*", psScalarAlloc(1.0 / (float)j, PS_TYPE_F32));
                     }
                 }