Changeset 9529 for trunk/psphot/src/models/pmModel_PGAUSS.c
- Timestamp:
- Oct 12, 2006, 4:36:54 PM (20 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/models/pmModel_PGAUSS.c (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/models/pmModel_PGAUSS.c
r8882 r9529 23 23 psF32 t = 1 + z + z*z/2.0; 24 24 psF32 r = 1.0 / (t + z*z*z/6.0); /* exp (-Z) */ 25 psF32 f = PAR[PM_PAR_ FLUX]*r + PAR[PM_PAR_SKY];25 psF32 f = PAR[PM_PAR_I0]*r + PAR[PM_PAR_SKY]; 26 26 27 27 if (deriv != NULL) { 28 // note difference from a pure gaussian: q = PAR[PM_PAR_ FLUX]*r29 psF32 q = PAR[PM_PAR_ FLUX]*r*r*t;28 // note difference from a pure gaussian: q = PAR[PM_PAR_I0]*r 29 psF32 q = PAR[PM_PAR_I0]*r*r*t; 30 30 deriv->data.F32[0] = +1.0; 31 31 deriv->data.F32[1] = +r; … … 92 92 norm *= 0.01; 93 93 94 psF64 Flux = params->data.F32[PM_PAR_ FLUX] * Area * norm;94 psF64 Flux = params->data.F32[PM_PAR_I0] * Area * norm; 95 95 96 96 return(Flux); … … 160 160 status = true; 161 161 status &= (dP < 0.5); 162 status &= (PAR[PM_PAR_ FLUX] > 0);163 status &= ((dPAR[PM_PAR_ FLUX]/PAR[PM_PAR_FLUX]) < 0.5);162 status &= (PAR[PM_PAR_I0] > 0); 163 status &= ((dPAR[PM_PAR_I0]/PAR[PM_PAR_I0]) < 0.5); 164 164 165 165 if (status) return true;
Note:
See TracChangeset
for help on using the changeset viewer.
