Changeset 9aac25d in sasmodels for sasmodels/models/fcc.c
- Timestamp:
- Feb 8, 2016 7:27:30 AM (8 years ago)
- Branches:
- master, core_shell_microgels, costrafo411, magnetic_model, release_v0.94, release_v0.95, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests
- Children:
- 98f3053
- Parents:
- 3882eeb
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sasmodels/models/fcc.c
reeb8bac r9aac25d 7 7 double _FCC_Integrand(double q, double dnn, double d_factor, double theta, double phi); 8 8 double _FCCeval(double Theta, double Phi, double temp1, double temp3); 9 double _sphereform(double q, double radius, double sld, double solvent_sld);10 9 11 10 … … 39 38 40 39 return (result); 41 }42 43 double _sphereform(double q, double radius, double sld, double solvent_sld){44 const double qr = q*radius;45 double sn, cn;46 SINCOS(qr, sn, cn);47 const double bes = (qr == 0.0 ? 1.0 : 3.0*(sn-qr*cn)/(qr*qr*qr));48 const double fq = bes * (sld - solvent_sld)*form_volume(radius);49 return 1.0e-4*fq*fq;50 40 } 51 41 … … 88 78 89 79 answer = (vb-va)/2.0*summ; 90 answer = answer* _sphereform(q,radius,sld,solvent_sld)*latticescale;80 answer = answer*sphere_form(q,radius,sld,solvent_sld)*latticescale; 91 81 92 82 return answer; … … 175 165 176 166 // Use SphereForm directly from libigor 177 answer = _sphereform(q,radius,sld,solvent_sld)*Zq*latticescale;167 answer = sphere_form(q,radius,sld,solvent_sld)*Zq*latticescale; 178 168 179 169 return answer;
Note: See TracChangeset
for help on using the changeset viewer.