Changeset 29490

Timestamp:

Oct 20, 2010, 8:54:42 AM (16 years ago)

Author:

eugene

Message:

update the source size analysis : CR now detected using M_minor and Kron mag; EXT now detected using only psf - kron mags; update visualization of source sizes; fix apresid visualization (kron vs apmag options); for stand-alone psphot, do not mask as SAT NAN pixels which are already masked; function for dynamic calculation of CR detection parameters

Location:

branches/eam_branches/ipp-20100823/psModules/src/imcombine

Files:

• pmSubtractionEquation.c (modified) (2 diffs)
• pmSubtractionKernels.c (modified) (1 diff)

branches/eam_branches/ipp-20100823/psModules/src/imcombine/pmSubtractionEquation.c

@@ -1104 +1104 @@
         // XXX TEST: try some constraint on the svd solution
         // solution = psMatrixSolveSVD(solution, sumMatrix, sumVector, NAN);
-        solution = psMatrixSolveSVD(solution, sumMatrix, sumVector, NAN);
+        // SINGLE solution
+        if (1) {
+            solution = psMatrixSolveSVD(solution, sumMatrix, sumVector, NAN);
+        } else {
+            psVector *PERM = NULL;
+            psImage *LU = psMatrixLUDecomposition(NULL, &PERM, sumMatrix);
+            solution = psMatrixLUSolution(solution, LU, sumVector, PERM);
+            psFree (LU);
+            psFree (PERM);
+        }
 
         for (int i = 0; i < solution->n; i++) {
@@ -1175 +1184 @@
         psVectorInit(solution, 0);
 
+        // DUAL solution
         solution = psMatrixSolveSVD(solution, sumMatrix, sumVector, NAN);
 

branches/eam_branches/ipp-20100823/psModules/src/imcombine/pmSubtractionKernels.c

@@ -182 +182 @@
     return true;
 }
+
+// XXX *** this code used the central pixel to force zero net flux,
+// Alard actually uses kernel(0) for this purpose (for even order, kernel[i] = kernel'[i] - kernel[0])
 
 static bool pmSubtractionKernelPreCalcNormalize(pmSubtractionKernels *kernels, pmSubtractionKernelPreCalc *preCalc,