Changeset 36033 for branches/eam_branches/ipp-20130711/psphot/src/psphotSourceFits.c

Timestamp:

Aug 26, 2013, 6:03:17 AM (13 years ago)

Author:

eugene

Message:

iterate over Reff and Sidx

File:

• branches/eam_branches/ipp-20130711/psphot/src/psphotSourceFits.c (modified) (3 diffs)

branches/eam_branches/ipp-20130711/psphot/src/psphotSourceFits.c

@@ -558 +558 @@
 }
 
-# define TIMING 1
+# define TIMING 0
 
 bool psphotSersicModelGuessPCM (pmPCMdata *pcm, pmSource *source, psImageMaskType maskVal, float psfSize);
@@ -759 +759 @@
 
 // float indexGuessInv[] = {0.5, 0.33, 0.25, 0.167, 0.125, 0.083};
-float indexGuessInv[] = {0.5, 0.25, 0.167, 0.125};
-# define N_INDEX_GUESS_INV 4
-
-float reffGuess[] = {3.0, 10.0, 30.0};
-# define N_REFF_GUESS 3
+float indexGuessInv[] = {0.5, 0.4, 0.3, 0.25, 0.20, 0.15, 0.125};
+# define N_INDEX_GUESS_INV 7
+
+float reffGuess[] = {3.0, 10.0, 20.0, 30.0, 40.0};
+# define N_REFF_GUESS 5
 
 // A sersic model is very sensitive to the index.  attempt to find the index first by grid search in just the index
@@ -882 +882 @@
                 rMin = reffGuess[j];
             }
-            fprintf (stderr, "%d | %f %f %f %f | %f %f %f %f", i, indexGuessInv[i], reffGuess[j], Io, Chisq, sMin, rMin, iMin, xMin);
-            fprintf (stderr, "\n");
+            // fprintf (stderr, "%d | %f %f %f %f | %f %f %f %f", i, indexGuessInv[i], reffGuess[j], Io, Chisq, sMin, rMin, iMin, xMin);
+            // fprintf (stderr, "\n");
         }
     }