Changes in sasmodels/models/pearl_necklace.c [4b541ac:2126131] in sasmodels
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sasmodels/models/pearl_necklace.c
r4b541ac r2126131 34 34 // But there is a 1/(1-sinc) term below which blows up so don't bother 35 35 const double q_edge = q * edge_sep; 36 const double beta = ( sas_Si(q*(A_s-radius)) - sas_Si(q*radius)) / q_edge;37 const double gamma = sas_Si(q_edge) / q_edge;38 const double psi = s as_3j1x_x(q*radius);36 const double beta = (Si(q*(A_s-radius)) - Si(q*radius)) / q_edge; 37 const double gamma = Si(q_edge) / q_edge; 38 const double psi = sph_j1c(q*radius); 39 39 40 40 // Precomputed sinc terms 41 const double si = s as_sinx_x(q*A_s);41 const double si = sinc(q*A_s); 42 42 const double omsi = 1.0 - si; 43 43 const double pow_si = pow(si, num_pearls); … … 54 54 - 2.0 * (1.0 - pow_si/si)*beta*beta / (omsi*omsi) 55 55 + 2.0 * num_strings*beta*beta / omsi 56 + num_strings * (2.0*gamma - square(s as_sinx_x(q_edge/2.0)))56 + num_strings * (2.0*gamma - square(sinc(q_edge/2.0))) 57 57 ); 58 58
Note: See TracChangeset
for help on using the changeset viewer.