Changes in sasmodels/models/vesicle.c [2c74c11:3a48772] in sasmodels
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sasmodels/models/vesicle.c
r2c74c11 r3a48772 8 8 { 9 9 //note that for the vesicle model, the volume is ONLY the shell volume 10 double volume; 11 volume =4.*M_PI*(radius+thickness)*(radius+thickness)*(radius+thickness)/3; 12 volume -=4.*M_PI*radius*radius*radius/3.; 13 return volume; 10 return M_4PI_3*(cube(radius+thickness) - cube(radius)); 14 11 } 15 12 … … 32 29 // core first, then add in shell 33 30 contrast = sld_solvent-sld; 34 vol = 4.0*M_PI/3.0*radius*radius*radius;35 f = vol *sph_j1c(q*radius)*contrast;31 vol = M_4PI_3*cube(radius); 32 f = vol * sph_j1c(q*radius) * contrast; 36 33 37 34 //now the shell. No volume normalization as this is done by the caller 38 35 contrast = sld-sld_solvent; 39 vol = 4.0*M_PI/3.0*(radius+thickness)*(radius+thickness)*(radius+thickness);40 f += vol *sph_j1c(q*(radius+thickness))*contrast;36 vol = M_4PI_3*cube(radius+thickness); 37 f += vol * sph_j1c(q*(radius+thickness)) * contrast; 41 38 42 39 //rescale to [cm-1]. 43 f2 = volfraction *f*f*1.0e-4;40 f2 = volfraction * f*f*1.0e-4; 44 41 45 return (f2);42 return f2; 46 43 }
Note: See TracChangeset
for help on using the changeset viewer.