- Timestamp:
- Dec 22, 2005, 4:26:41 PM (21 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/eam_rel9_b0/psModules/src/objects/models/pmModel_QGAUSS.c
r5832 r5833 31 31 32 32 if (deriv != NULL) { 33 psF32 *dPAR = deriv->data.F32; 34 33 35 // note difference from a pure gaussian: q = params->data.F32[1]*r 34 36 psF32 t = PAR[1]*r*r; 35 37 psF32 q = t*(PAR[7] + 2.25*pow(z, 1.25)); 36 38 37 d eriv->data.F32[0] = +1.0;38 d eriv->data.F32[1] = +r;39 d eriv->data.F32[2] = q*(2.0*px*PAR[4] + PAR[6]*Y);40 d eriv->data.F32[3] = q*(2.0*py*PAR[5] + PAR[6]*X);41 d eriv->data.F32[4] = -2.0*q*px*X;42 d eriv->data.F32[5] = -2.0*q*py*Y;43 d eriv->data.F32[6] = -q*X*Y;44 d eriv->data.F32[7] = -t*z;39 dPAR[0] = +1.0; 40 dPAR[1] = +r; 41 dPAR[2] = q*(2.0*px*PAR[4] + PAR[6]*Y); 42 dPAR[3] = q*(2.0*py*PAR[5] + PAR[6]*X); 43 dPAR[4] = -2.0*q*px*X; 44 dPAR[5] = -2.0*q*py*Y; 45 dPAR[6] = -q*X*Y; 46 dPAR[7] = -t*z; 45 47 } 46 48 return(f);
Note:
See TracChangeset
for help on using the changeset viewer.
