Changeset 14799 for trunk/psphot/src/psphotModelTest.c
- Timestamp:
- Sep 10, 2007, 9:15:36 AM (19 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/psphotModelTest.c (modified) (6 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/psphotModelTest.c
r14655 r14799 6 6 7 7 bool status; 8 int modelType ;8 int modelType = -1; 9 9 float obsMag, fitMag, value; 10 10 char name[64]; … … 63 63 // in fitMode, psf sets the model type 64 64 if (fitMode == PM_SOURCE_FIT_PSF) { 65 psf = psphotLoadPSF (config, view, recipe);65 psf = psphotLoadPSF (config, view, recipe); 66 66 if (!psf) psAbort("PSF_INPUT_FILE not supplied"); 67 67 modelType = psf->type; 68 source->type = PM_SOURCE_TYPE_STAR;68 source->type = PM_SOURCE_TYPE_STAR; 69 69 } 70 70 if (fitMode == PM_SOURCE_FIT_EXT) { … … 91 91 modelType = pmModelClassGetType (modelName); 92 92 if (modelType < 0) psAbort("unknown model %s", modelName); 93 source->type = PM_SOURCE_TYPE_EXTENDED;93 source->type = PM_SOURCE_TYPE_EXTENDED; 94 94 } 95 95 if (fitMode == PM_SOURCE_FIT_PSF_X_EXT) { 96 // we need to load BOTH a psf and an ext model97 psf = psphotLoadPSF (config, view, recipe);96 // we need to load BOTH a psf and an ext model 97 psf = psphotLoadPSF (config, view, recipe); 98 98 if (!psf) psAbort("PSF_INPUT_FILE not supplied"); 99 99 … … 120 120 modelType = pmModelClassGetType (modelName); 121 121 if (modelType < 0) psAbort("unknown model %s", modelName); 122 source->type = PM_SOURCE_TYPE_EXTENDED;122 source->type = PM_SOURCE_TYPE_EXTENDED; 123 123 } 124 124 … … 150 150 int nParams = pmModelClassParameterCount (modelType); 151 151 psF32 *params = model->params->data.F32; 152 params[PM_PAR_XPOS] = xObj; // XXX use the user-supplied value, 152 params[PM_PAR_XPOS] = xObj; // XXX use the user-supplied value, 153 153 params[PM_PAR_YPOS] = yObj; // XXX or use the centroid 154 154 for (int i = 0; i < nParams; i++) { … … 199 199 200 200 if (fitMode == PM_SOURCE_FIT_PSF_X_EXT) { 201 // build the psf for the object202 source->modelPSF = pmModelFromPSF (model, psf);203 source->modelEXT = model;204 status = psphotPSFConvModel (source, recipe, maskVal);205 model = source->modelConv;206 params = model->params->data.F32;201 // build the psf for the object 202 source->modelPSF = pmModelFromPSF (model, psf); 203 source->modelEXT = model; 204 status = psphotPSFConvModel (source, recipe, maskVal); 205 model = source->modelConv; 206 params = model->params->data.F32; 207 207 } else { 208 status = pmSourceFitModel (source, model, fitMode, maskVal);208 status = pmSourceFitModel (source, model, fitMode, maskVal); 209 209 } 210 210
Note:
See TracChangeset
for help on using the changeset viewer.
