Changeset 10075 for trunk/psModules/src/objects/pmSource.c

Timestamp:

Nov 17, 2006, 9:31:44 PM (20 years ago)

Author:

magnier

Message:

working on new linear fit process

File:

• trunk/psModules/src/objects/pmSource.c (modified) (16 diffs)

trunk/psModules/src/objects/pmSource.c

@@ -6 +6 @@
  *  @author EAM, IfA: significant modifications.
  *
- *  @version $Revision: 1.18 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2006-11-17 23:01:30 $
+ *  @version $Revision: 1.19 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-11-18 07:31:44 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -31 +31 @@
 static void sourceFree(pmSource *tmp)
 {
-    psTrace(__func__, 5, "---- begin ----\n");
+    psTrace("psModules.objects", 5, "---- begin ----\n");
     psFree(tmp->peak);
     psFree(tmp->pixels);
@@ -40 +40 @@
     psFree(tmp->modelEXT);
     psFree(tmp->blends);
-    psTrace(__func__, 5, "---- end ----\n");
+    psTrace("psModules.objects", 5, "---- end ----\n");
 }
 
@@ -65 +65 @@
 pmSource *pmSourceAlloc()
 {
-    psTrace(__func__, 5, "---- begin ----\n");
+    psTrace("psModules.objects", 5, "---- begin ----\n");
     static int id = 1;
     pmSource *source = (pmSource *) psAlloc(sizeof(pmSource));
@@ -90 +90 @@
     source->pixWeight = NAN;
 
-    psTrace(__func__, 5, "---- end ----\n");
+    psTrace("psModules.objects", 5, "---- end ----\n");
     return(source);
 }
@@ -205 +205 @@
 pmPSFClump pmSourcePSFClump(psArray *sources, psMetadata *recipe)
 {
-    psTrace(__func__, 5, "---- begin ----\n");
+    psTrace("psModules.objects", 5, "---- begin ----\n");
 
     # define NPIX 10
@@ -291 +291 @@
         stats = psImageStats (stats, splane, NULL, 0);
         peaks = pmFindImagePeaks (splane, stats[0].max / 2);
-        psTrace (__func__, 2, "clump threshold is %f\n", stats[0].max/2);
+        psTrace ("psModules.objects", 2, "clump threshold is %f\n", stats[0].max/2);
 
         const bool keep_psf_clump = psMetadataLookupBool(NULL, recipe, "KEEP_PSF_CLUMP");
@@ -323 +323 @@
         psArraySort (peaks, pmPeaksCompareDescend);
         clump = peaks->data[0];
-        psTrace (__func__, 2, "clump is at %d, %d (%f)\n", clump->x, clump->y, clump->counts);
+        psTrace ("psModules.objects", 2, "clump is at %d, %d (%f)\n", clump->x, clump->y, clump->counts);
 
         // define section window for clump
@@ -369 +369 @@
         psfClump.dY = stats->clippedStdev;
 
-        psTrace (__func__, 2, "clump  X,  Y: %f, %f\n", psfClump.X, psfClump.Y);
-        psTrace (__func__, 2, "clump DX, DY: %f, %f\n", psfClump.dX, psfClump.dY);
+        psTrace ("psModules.objects", 2, "clump  X,  Y: %f, %f\n", psfClump.X, psfClump.Y);
+        psTrace ("psModules.objects", 2, "clump DX, DY: %f, %f\n", psfClump.dX, psfClump.dY);
         // these values should be pushed on the recipe somewhere
 
@@ -379 +379 @@
     }
 
-    psTrace(__func__, 5, "---- end ----\n");
+    psTrace("psModules.objects", 5, "---- end ----\n");
     return (psfClump);
 }
@@ -391 +391 @@
 bool pmSourceRoughClass(psArray *sources, psMetadata *recipe, pmPSFClump clump)
 {
-    psTrace(__func__, 5, "---- begin ----");
+    psTrace("psModules.objects", 5, "---- begin ----");
 
     psBool rc = true;
@@ -498 +498 @@
     psFree (starsn);
 
-    psTrace (__func__, 2, "Nstar:    %3d\n", Nstar);
-    psTrace (__func__, 2, "Npsf:     %3d\n", Npsf);
-    psTrace (__func__, 2, "Next:     %3d\n", Next);
-    psTrace (__func__, 2, "Nsatstar: %3d\n", Nsatstar);
-    psTrace (__func__, 2, "Nsat:     %3d\n", Nsat);
-    psTrace (__func__, 2, "Ncr:      %3d\n", Ncr);
-
-    psTrace(__func__, 5, "---- end ----\n");
+    psTrace ("psModules.objects", 2, "Nstar:    %3d\n", Nstar);
+    psTrace ("psModules.objects", 2, "Npsf:     %3d\n", Npsf);
+    psTrace ("psModules.objects", 2, "Next:     %3d\n", Next);
+    psTrace ("psModules.objects", 2, "Nsatstar: %3d\n", Nsatstar);
+    psTrace ("psModules.objects", 2, "Nsat:     %3d\n", Nsat);
+    psTrace ("psModules.objects", 2, "Ncr:      %3d\n", Ncr);
+
+    psTrace("psModules.objects", 5, "---- end ----\n");
     return(rc);
 }
@@ -531 +531 @@
                      psF32 radius)
 {
-    psTrace(__func__, 5, "---- begin ----\n");
+    psTrace("psModules.objects", 5, "---- begin ----\n");
     PS_ASSERT_PTR_NON_NULL(source, false);
     PS_ASSERT_PTR_NON_NULL(source->peak, false);
@@ -637 +637 @@
     // XXX EAM - the limit is a bit arbitrary.  make it user defined?
     if ((numPixels < 0.75*R2) || (Sum <= 0)) {
-        psTrace (__func__, 3, "no valid pixels for source\n");
-        psTrace(__func__, 5, "---- end (false) ----\n");
+        psTrace ("psModules.objects", 3, "no valid pixels for source\n");
+        psTrace("psModules.objects", 5, "---- end (false) ----\n");
         return (false);
     }
 
-    psTrace (__func__, 4, "sky: %f  Sum: %f  X1: %f  Y1: %f  X2: %f  Y2: %f  XY: %f  Npix: %d\n",
+    psTrace ("psModules.objects", 4, "sky: %f  Sum: %f  X1: %f  Y1: %f  X2: %f  Y2: %f  XY: %f  Npix: %d\n",
              sky, Sum, X1, Y1, X2, Y2, XY, numPixels);
 
@@ -653 +653 @@
     y = Y1/Sum;
     if ((fabs(x) > radius) || (fabs(y) > radius)) {
-        psTrace (__func__, 3, "large centroid swing; invalid peak %d, %d\n",
+        psTrace ("psModules.objects", 3, "large centroid swing; invalid peak %d, %d\n",
                  source->peak->x, source->peak->y);
-        psTrace(__func__, 5, "---- end(false)  ----\n");
+        psTrace("psModules.objects", 5, "---- end(false)  ----\n");
         return (false);
     }
@@ -673 +673 @@
     source->moments->Sy = sqrt(PS_MAX(Y2/Sum - PS_SQR(y), 0));
 
-    psTrace (__func__, 4,
+    psTrace ("psModules.objects", 4,
              "sky: %f  Sum: %f  x: %f  y: %f  Sx: %f  Sy: %f  Sxy: %f\n",
              sky, Sum, source->moments->x, source->moments->y,
              source->moments->Sx, source->moments->Sy, source->moments->Sxy);
 
-    psTrace(__func__, 5, "---- end ----\n");
+    psTrace("psModules.objects", 5, "---- end ----\n");
     return(true);
 }