Changeset 1e7b0db0 in sasmodels for sasmodels/models/core_shell_bicelle_elliptical.c
- Timestamp:
- Jan 11, 2017 5:57:04 AM (7 years ago)
- Branches:
- master, core_shell_microgels, costrafo411, magnetic_model, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests
- Children:
- 925ad6e
- Parents:
- 473a9f1
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sasmodels/models/core_shell_bicelle_elliptical.c
rfcb33e4 r1e7b0db0 76 76 double sinarg1 = qq*halfheight*cos_alpha; 77 77 double sinarg2 = qq*(halfheight+facthick)*cos_alpha; 78 si1 = s inc(sinarg1);79 si2 = s inc(sinarg2);78 si1 = sas_sinx_x(sinarg1); 79 si2 = sas_sinx_x(sinarg2); 80 80 for(int j=0;j<76;j++) { 81 81 //76 gauss points for the inner integral (WAS 20 points,so this may make unecessarily slow, but playing safe) … … 131 131 const double be1 = sas_J1c(qq*r); 132 132 const double be2 = sas_J1c( qq*(r + radthick ) ); 133 const double si1 = s inc( qq*halfheight*cos_val );134 const double si2 = s inc( qq*(halfheight + facthick)*cos_val );133 const double si1 = sas_sinx_x( qq*halfheight*cos_val ); 134 const double si2 = sas_sinx_x( qq*(halfheight + facthick)*cos_val ); 135 135 const double Aq = square( vol1*dr1*si1*be1 + vol2*dr2*si2*be2 + vol3*dr3*si2*be1); 136 136 //const double vol = form_volume(radius_minor, r_ratio, length);
Note: See TracChangeset
for help on using the changeset viewer.