Changeset 9529 for trunk/psphot/src/models/pmModel_ZGAUSS.c
- Timestamp:
- Oct 12, 2006, 4:36:54 PM (20 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/models/pmModel_ZGAUSS.c (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/models/pmModel_ZGAUSS.c
r8882 r9529 32 32 psF32 p = pow(z, PAR[PM_PAR_7] - 1.0); 33 33 psF32 r = 1.0 / (1 + z*p + SQ(SQ(pr))); 34 psF32 f = PAR[PM_PAR_ FLUX]*r + PAR[PM_PAR_SKY];34 psF32 f = PAR[PM_PAR_I0]*r + PAR[PM_PAR_SKY]; 35 35 36 36 if (deriv != NULL) { 37 37 // note difference from a pure gaussian: q = params->data.F32[1]*r 38 psF32 t = PAR[PM_PAR_ FLUX]*r*r;38 psF32 t = PAR[PM_PAR_I0]*r*r; 39 39 psF32 q = t*(PAR[PM_PAR_7]*p + 4*PAR8*pr*pr*pr); 40 40 … … 72 72 norm *= 0.01; 73 73 74 psF64 Flux = PAR[PM_PAR_ FLUX] * Area * norm;74 psF64 Flux = PAR[PM_PAR_I0] * Area * norm; 75 75 76 76 return(Flux); … … 89 89 90 90 if (flux <= 0) return (1.0); 91 if (PAR[PM_PAR_ FLUX] <= 0) return (1.0);92 if (flux >= PAR[PM_PAR_ FLUX]) return (1.0);91 if (PAR[PM_PAR_I0] <= 0) return (1.0); 92 if (flux >= PAR[PM_PAR_I0]) return (1.0); 93 93 94 94 // convert Sx,Sy,Sxy to major/minor axes … … 99 99 axes = EllipseShapeToAxes (shape); 100 100 psF64 dr = 1.0 / axes.major; 101 psF64 limit = flux / PAR[PM_PAR_ FLUX];101 psF64 limit = flux / PAR[PM_PAR_I0]; 102 102 103 103 // XXX : we can do this faster with an intelligent starting choice
Note:
See TracChangeset
for help on using the changeset viewer.
