Changeset 6379 for trunk/psphot/src/psphotApResid.c
- Timestamp:
- Feb 7, 2006, 8:52:03 PM (20 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/psphotApResid.c (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/psphotApResid.c
r6056 r6379 1 1 # include "psphot.h" 2 static char DEFAULT_OPTION[] = "SKYBIAS"; 2 3 3 4 // measure the aperture residual statistics … … 69 70 // APTREND options : NONE SKYBIAS XY_LIN XY_QUAD SKY_XY_LIN FULL 70 71 char *ApTrendOption = psMetadataLookupPtr (&status, config, "APTREND"); 71 if (!status) ApTrendOption = psStringCopy ("SKYBIAS");72 if (!status) ApTrendOption = DEFAULT_OPTION; 72 73 73 74 // 3hi/1lo sigma clipping on the rflux vs metric fit … … 190 191 psf->ApTrend = psVectorChiClipFitPolynomial4D (psf->ApTrend, stats, mask, PSFTRY_MASK_ALL, apResid, dMag, xPos, yPos, r2rflux, flux); 191 192 } 192 193 # if (1)194 psPolynomial4D *poly = psf->ApTrend;195 for (int nt = 0; nt <= poly->nT; nt++) {196 for (int nz = 0; nz <= poly->nZ; nz++) {197 for (int ny = 0; ny <= poly->nY; ny++) {198 for (int nx = 0; nx <= poly->nX; nx++) {199 if (poly->mask[nx][ny][nz][nt]) continue;200 fprintf (stderr, "%d %d %d %d : %22.15g\n", nx, ny, nz, nt, poly->coeff[nx][ny][nz][nt]);201 }202 }203 }204 }205 # endif206 193 207 194 // construct the fitted values and the residuals … … 256 243 1e6*psf->ApTrend->coeff[0][2][0][0]); 257 244 258 // psFree (stats); 259 // psFree (mask); 260 // psFree (rflux); 261 // psFree (apResid); 262 245 psFree (mask); 246 psFree (xPos); 247 psFree (yPos); 248 psFree (flux); 249 psFree (r2rflux); 250 psFree (apResid); 251 psFree (dMag); 252 253 psFree (fit); 254 psFree (resid); 255 psFree (stats); 256 psFree (residStats); 263 257 return true; 264 258 }
Note:
See TracChangeset
for help on using the changeset viewer.
