Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • sasmodels/models/pearl_necklace.c

    r2126131 r4b541ac  
    3434    // But there is a 1/(1-sinc) term below which blows up so don't bother 
    3535    const double q_edge = q * edge_sep; 
    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); 
     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 = sas_3j1x_x(q*radius); 
    3939 
    4040    // Precomputed sinc terms 
    41     const double si = sinc(q*A_s); 
     41    const double si = sas_sinx_x(q*A_s); 
    4242    const double omsi = 1.0 - si; 
    4343    const double pow_si = pow(si, num_pearls); 
     
    5454        - 2.0 * (1.0 - pow_si/si)*beta*beta / (omsi*omsi) 
    5555        + 2.0 * num_strings*beta*beta / omsi 
    56         + num_strings * (2.0*gamma - square(sinc(q_edge/2.0))) 
     56        + num_strings * (2.0*gamma - square(sas_sinx_x(q_edge/2.0))) 
    5757        ); 
    5858 
Note: See TracChangeset for help on using the changeset viewer.