Changeset a807206 in sasmodels for sasmodels/models/polymer_micelle.c
- Timestamp:
- Sep 30, 2016 10:42:06 PM (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:
- caddb14, 5031ca3
- Parents:
- 2222134
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sasmodels/models/polymer_micelle.c
r0bef47b ra807206 7 7 double corona_sld, 8 8 double radius_core, 9 double r adius_gyr,9 double rg, 10 10 double d_penetration, 11 11 double n_aggreg); … … 19 19 double corona_sld, 20 20 double radius_core, 21 double r adius_gyr,21 double rg, 22 22 double d_penetration, 23 23 double n_aggreg) … … 35 35 36 36 // Self-correlation term of the chains 37 const double qrg2 = q*r adius_gyr*q*radius_gyr;37 const double qrg2 = q*rg*q*rg; 38 38 const double debye_chain = (qrg2 == 0.0) ? 1.0 : 2.0*(expm1(-qrg2)+qrg2)/(qrg2*qrg2); 39 39 const double term2 = n_aggreg * beta_corona * beta_corona * debye_chain; … … 41 41 // Interference cross-term between core and chains 42 42 const double chain_ampl = (qrg2 == 0.0) ? 1.0 : -expm1(-qrg2)/qrg2; 43 const double bes_corona = sinc(q*(radius_core + d_penetration * r adius_gyr));43 const double bes_corona = sinc(q*(radius_core + d_penetration * rg)); 44 44 const double term3 = 2 * n_aggreg * n_aggreg * beta_core * beta_corona * 45 45 bes_core * chain_ampl * bes_corona; … … 70 70 double corona_sld, 71 71 double radius_core, 72 double r adius_gyr,72 double rg, 73 73 double d_penetration, 74 74 double n_aggreg) … … 82 82 corona_sld, 83 83 radius_core, 84 r adius_gyr,84 rg, 85 85 d_penetration, 86 86 n_aggreg);
Note: See TracChangeset
for help on using the changeset viewer.