Changeset 9529 for trunk/psphot/src/models/pmModel_RGAUSS.c
- Timestamp:
- Oct 12, 2006, 4:36:54 PM (20 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/models/pmModel_RGAUSS.c (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/models/pmModel_RGAUSS.c
r8882 r9529 30 30 psF32 p = pow(z, PAR[PM_PAR_7] - 1.0); 31 31 psF32 r = 1.0 / (1 + z + z*p); 32 psF32 f = PAR[PM_PAR_ FLUX]*r + PAR[PM_PAR_SKY];32 psF32 f = PAR[PM_PAR_I0]*r + PAR[PM_PAR_SKY]; 33 33 34 34 if (deriv != NULL) { 35 35 // note difference from a pure gaussian: q = params->data.F32[1]*r 36 psF32 t = PAR[PM_PAR_ FLUX]*r*r;36 psF32 t = PAR[PM_PAR_I0]*r*r; 37 37 psF32 q = t*(1 + PAR[PM_PAR_7]*p); 38 38 … … 87 87 88 88 if (flux <= 0) return (1.0); 89 if (PAR[PM_PAR_ FLUX] <= 0) return (1.0);90 if (flux >= PAR[PM_PAR_ FLUX]) return (1.0);89 if (PAR[PM_PAR_I0] <= 0) return (1.0); 90 if (flux >= PAR[PM_PAR_I0]) return (1.0); 91 91 92 92 // if Sx == Sy, sigma = Sx == Sy 93 93 psF64 sigma = hypot (1.0 / PAR[PM_PAR_SXX], 1.0 / PAR[PM_PAR_SYY]) / sqrt(2.0); 94 94 psF64 dz = 1.0 / (2.0 * sigma*sigma); 95 psF64 limit = flux / PAR[PM_PAR_ FLUX];95 psF64 limit = flux / PAR[PM_PAR_I0]; 96 96 97 97 // we can do this much better with intelligent choices here
Note:
See TracChangeset
for help on using the changeset viewer.
