Changeset 14814 for trunk/ppSim/src/ppSimSetPSF.c
- Timestamp:
- Sep 11, 2007, 12:31:46 PM (19 years ago)
- File:
-
- 1 edited
-
trunk/ppSim/src/ppSimSetPSF.c (modified) (6 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/ppSim/src/ppSimSetPSF.c
r14657 r14814 11 11 psf = psMetadataLookupPtr (&status, chip->analysis, "PSPHOT.PSF"); 12 12 if (psf) { 13 return true;13 return true; 14 14 } 15 15 … … 20 20 char *psfModelName = psMetadataLookupStr(&status, config->arguments, "PSF.MODEL"); // Seeing sigma (pix) 21 21 if (psfModelName == NULL) { 22 psfModelName = defaultModel;22 psfModelName = defaultModel; 23 23 } 24 24 25 25 pmModelType type = pmModelClassGetType (psfModelName); 26 26 if (type == -1) { 27 psError (PS_ERR_UNKNOWN, false, "invalid model name");27 psError (PS_ERR_UNKNOWN, false, "invalid model name"); 28 28 return false; 29 29 } … … 32 32 psf = pmPSFAlloc (type, true, NULL); 33 33 for (int i = 0; i < psf->params->n; i++) { 34 34 35 35 if (i == PM_PAR_SKY) continue; 36 36 if (i == PM_PAR_I0) continue; … … 38 38 if (i == PM_PAR_YPOS) continue; 39 39 40 // generate a 'constant' polynomial, unmasked40 // generate a 'constant' polynomial, unmasked 41 41 param = psPolynomial2DAlloc(PS_POLYNOMIAL_ORD, 0, 0); 42 param->mask[0][0] = 0;42 param->mask[0][0] = 0; 43 43 44 44 psf->params->data[i] = param; … … 47 47 psEllipseAxes axes; 48 48 psEllipsePol pol; 49 49 50 50 // supply the semi-major axis 51 51 axes.major = seeing / 2.0; … … 65 65 66 66 if (!strcasecmp (psfModelName, "PS_MODEL_QGAUSS")) { 67 param = psf->params->data[PM_PAR_7];68 param->coeff[0][0] = 1.0;67 param = psf->params->data[PM_PAR_7]; 68 param->coeff[0][0] = 1.0; 69 69 } 70 70 71 71 if (!strcasecmp (psfModelName, "PS_MODEL_RGAUSS")) { 72 param = psf->params->data[PM_PAR_7];73 param->coeff[0][0] = 1.0;72 param = psf->params->data[PM_PAR_7]; 73 param->coeff[0][0] = 1.0; 74 74 } 75 75 76 76 psMetadataAdd (chip->analysis, PS_LIST_TAIL, "PSPHOT.PSF", PS_DATA_UNKNOWN, "psphot psf", psf); 77 psFree(psf); // Drop reference 78 77 79 return true; 78 80 }
Note:
See TracChangeset
for help on using the changeset viewer.
