- Timestamp:
- Feb 17, 2006, 7:13:42 AM (20 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/rel10_ifa/psModules/src/objects/models/pmModel_PGAUSS.c
r5257 r6448 26 26 27 27 if (deriv != NULL) { 28 // note difference from a pure gaussian: q = PAR[1]*r28 psF32 *dPAR = deriv->data.F32; 29 29 psF32 q = PAR[1]*r*r*t; 30 d eriv->data.F32[0] = +1.0;31 d eriv->data.F32[1] = +r;32 d eriv->data.F32[2] = q*(2.0*px*PAR[4] + PAR[6]*Y);33 d eriv->data.F32[3] = q*(2.0*py*PAR[5] + PAR[6]*X);34 d eriv->data.F32[4] = -2.0*q*px*X;35 d eriv->data.F32[5] = -2.0*q*py*Y;36 d eriv->data.F32[6] = -q*X*Y;30 dPAR[0] = +1.0; 31 dPAR[1] = +r; 32 dPAR[2] = q*(2.0*px*PAR[4] + PAR[6]*Y); 33 dPAR[3] = q*(2.0*py*PAR[5] + PAR[6]*X); 34 dPAR[4] = -2.0*q*px*X; 35 dPAR[5] = -2.0*q*py*Y; 36 dPAR[6] = -q*X*Y; 37 37 } 38 38 return(f); … … 47 47 48 48 beta_lim[0][0].data.F32[0] = 1000; 49 beta_lim[0][0].data.F32[1] = 10000;49 beta_lim[0][0].data.F32[1] = 3e6; 50 50 beta_lim[0][0].data.F32[2] = 5; 51 51 beta_lim[0][0].data.F32[3] = 5; … … 63 63 64 64 params_max[0][0].data.F32[0] = 1e5; 65 params_max[0][0].data.F32[1] = 1e 6;65 params_max[0][0].data.F32[1] = 1e8; 66 66 params_max[0][0].data.F32[2] = 1e4; // this should be set by image dimensions! 67 67 params_max[0][0].data.F32[3] = 1e4; // this should be set by image dimensions! … … 130 130 params[5] = 1.2 / moments->Sy; 131 131 params[6] = 0.0; 132 132 133 return(true); 133 134 }
Note:
See TracChangeset
for help on using the changeset viewer.
