Changeset 20307 for trunk/psModules/src/objects/pmSourcePhotometry.c
- Timestamp:
- Oct 21, 2008, 4:11:36 PM (18 years ago)
- File:
-
- 1 edited
-
trunk/psModules/src/objects/pmSourcePhotometry.c (modified) (5 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psModules/src/objects/pmSourcePhotometry.c
r20077 r20307 3 3 * @author EAM, IfA; GLG, MHPCC 4 4 * 5 * @version $Revision: 1.4 5$ $Name: not supported by cvs2svn $6 * @date $Date: 2008-10- 13 01:59:05$5 * @version $Revision: 1.46 $ $Name: not supported by cvs2svn $ 6 * @date $Date: 2008-10-22 02:11:08 $ 7 7 * 8 8 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 122 122 if (source->modelFits) { 123 123 for (int i = 0; i < source->modelFits->n; i++) { 124 pmModel *model = source->modelFits->data[i];125 status = pmSourcePhotometryModel (&model->mag, model);124 pmModel *model = source->modelFits->data[i]; 125 status = pmSourcePhotometryModel (&model->mag, model); 126 126 } 127 127 if (source->modelEXT) { 128 source->extMag = source->modelEXT->mag;128 source->extMag = source->modelEXT->mag; 129 129 } 130 130 } else { 131 131 if (source->modelEXT) { 132 status = pmSourcePhotometryModel (&source->extMag, source->modelEXT);132 status = pmSourcePhotometryModel (&source->extMag, source->modelEXT); 133 133 } 134 134 } … … 177 177 178 178 if (!psImageShiftMask(&flux, &mask, source->pixels, source->maskObj, maskVal, dx, dy, 179 NAN, 0xff, PS_INTERPOLATE_BI CUBE)) {179 NAN, 0xff, PS_INTERPOLATE_BIQUADRATIC)) { 180 180 // Not much we can do about it 181 181 psErrorClear(); … … 256 256 } 257 257 if (mode & PM_SOURCE_PHOT_APCORR) { 258 rflux = pow (10.0, 0.4*source->psfMag);259 source->apMag -= PS_SQR(model->radiusFit)*rflux * psf->skyBias + psf->skySat / rflux;258 rflux = pow (10.0, 0.4*source->psfMag); 259 source->apMag -= PS_SQR(model->radiusFit)*rflux * psf->skyBias + psf->skySat / rflux; 260 260 } 261 261 } … … 333 333 } 334 334 if (apSum <= 0) { 335 *apMag = NAN;335 *apMag = NAN; 336 336 return true; 337 337 }
Note:
See TracChangeset
for help on using the changeset viewer.
