Changeset 8246 for trunk/psModules/src/detrend/pmFringeStats.c
- Timestamp:
- Aug 8, 2006, 4:37:08 PM (20 years ago)
- File:
-
- 1 edited
-
trunk/psModules/src/detrend/pmFringeStats.c (modified) (8 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psModules/src/detrend/pmFringeStats.c
r7887 r8246 3 3 * @author Eugene Magnier, IfA 4 4 * 5 * @version $Revision: 1.1 5$ $Name: not supported by cvs2svn $6 * @date $Date: 2006-0 7-13 02:21:25$5 * @version $Revision: 1.16 $ $Name: not supported by cvs2svn $ 6 * @date $Date: 2006-08-09 02:37:07 $ 7 7 * 8 8 * Copyright 2004 IfA … … 327 327 dfPt[i] = 1.0 / medianSd->sampleStdev; 328 328 329 psTrace( __func__, 7, "[%d:%d,%d:%d]: %f %f\n", (int)region.x0, (int)region.x1,329 psTrace("psModules.detrend", 7, "[%d:%d,%d:%d]: %f %f\n", (int)region.x0, (int)region.x1, 330 330 (int)region.y0, (int)region.y1, fPt[i], dfPt[i]); 331 331 } … … 671 671 int numClipped = 0; // Number clipped 672 672 for (int i = 0; i < diffs->n; i++) { 673 psTrace( __func__, 10, "Region %d (%d): %f\n", i, mask->data.U8[i], diffs->data.F32[i]);673 psTrace("psModules.detrend", 10, "Region %d (%d): %f\n", i, mask->data.U8[i], diffs->data.F32[i]); 674 674 if (!mask->data.U8[i] && fabs(diffs->data.F32[i]) > middle + thresh) { 675 psTrace( __func__, 5, "Masking %d: %f\n", i, diffs->data.F32[i]);675 psTrace("psModules.detrend", 5, "Masking %d: %f\n", i, diffs->data.F32[i]); 676 676 mask->data.U8[i] = 1; 677 677 numClipped++; … … 732 732 if (!finite(fringe->f->data.F32[j])) { 733 733 mask->data.U8[j] = 1; 734 psTrace( __func__, 9, "Masking region %d because not finite in fringe %d.\n", j, i);734 psTrace("psModules.detrend", 9, "Masking region %d because not finite in fringe %d.\n", j, i); 735 735 } 736 736 } … … 749 749 fringeScaleDiffs(diff, science, fringes, scale); 750 750 numClipped = clipRegions(diff, mask, 3.0*rej); 751 psTrace( __func__, 4, "%d regions clipped in initial pass.\n", numClipped);751 psTrace("psModules.detrend", 4, "%d regions clipped in initial pass.\n", numClipped); 752 752 753 753 unsigned int iter = 0; // Iteration number … … 756 756 iter++; 757 757 scaleMeasure(scale, science, fringes); // The scales 758 psTrace( __func__, 1, "Fringe scales after iteration %d:\n", iter);759 psTrace( __func__, 1, "Background: %f %f\n", scale->coeff->data.F32[0], scale->coeffErr->data.F32[0]);758 psTrace("psModules.detrend", 1, "Fringe scales after iteration %d:\n", iter); 759 psTrace("psModules.detrend", 1, "Background: %f %f\n", scale->coeff->data.F32[0], scale->coeffErr->data.F32[0]); 760 760 for (int i = 0; i < scale->nFringeFrames; i++) { 761 psTrace( __func__, 1, "%d: %f %f\n", i, scale->coeff->data.F32[i + 1],761 psTrace("psModules.detrend", 1, "%d: %f %f\n", i, scale->coeff->data.F32[i + 1], 762 762 scale->coeffErr->data.F32[i + 1]); 763 763 } … … 766 766 iterClip = clipRegions(diff, mask, rej); // Number clipped 767 767 numClipped += iterClip; 768 psTrace( __func__, 9, "Clipped: %d\tFrac: %f\n", iterClip, (float)numClipped/(float)numRegions);768 psTrace("psModules.detrend", 9, "Clipped: %d\tFrac: %f\n", iterClip, (float)numClipped/(float)numRegions); 769 769 } while (iterClip > 0 && iter < nIter && (float)numClipped/(float)numRegions <= 1.0 - keepFrac); 770 770 psFree(diff); … … 811 811 psFree(scienceStats); 812 812 813 psTrace( __func__, 7, "Fringe solution:\n");813 psTrace("psModules.detrend", 7, "Fringe solution:\n"); 814 814 for (int i = 0; i < fringeImages->n + 1; i++) { 815 psTrace( __func__, 7, "%d: %f %f\n", i, scale->coeff->data.F32[i],815 psTrace("psModules.detrend", 7, "%d: %f %f\n", i, scale->coeff->data.F32[i], 816 816 scale->coeffErr->data.F32[i]); 817 817 }
Note:
See TracChangeset
for help on using the changeset viewer.
