Changeset ce27e21 in sasmodels for sasmodels/models/cylinder.c


Ignore:
Timestamp:
Aug 24, 2014 5:18:14 PM (10 years ago)
Author:
Paul Kienzle <pkienzle@…>
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:
1780d59
Parents:
14de349
Message:

first pass for sasview wrapper around opencl models

File:
1 edited

Legend:

Unmodified
Added
Removed
  • sasmodels/models/cylinder.c

    r14de349 rce27e21  
    140140DOCUMENTATION END 
    141141*/ 
    142 real volume(real radius, real length); 
     142real form_volume(real radius, real length); 
    143143real Iq(real q, real sld, real solvent_sld, real radius, real length); 
    144144real Iqxy(real qx, real qy, real sld, real solvent_sld, real radius, real length, real theta, real phi); 
    145145 
    146146 
    147 real volume(real radius, real length) 
     147real form_volume(real radius, real length) 
    148148{ 
    149149    return M_PI*radius*radius*length; 
    150150} 
    151  
    152151 
    153152real Iq(real q, 
     
    170169    // NOTE that for this (Fournet) definition of the integral, one must MULTIPLY by Vcyl 
    171170    // The additional volume factor is for polydisperse volume normalization. 
    172     const real s = (sld - solvent_sld) * volume(radius, length); 
     171    const real s = (sld - solvent_sld) * form_volume(radius, length); 
    173172    return REAL(1.0e-4) * form * s * s; 
    174173} 
     
    203202    const real bj = (besarg == REAL(0.0) ? REAL(0.5) : J1(besarg)/besarg); 
    204203    const real si = (siarg == REAL(0.0) ? REAL(1.0) : sin(siarg)/siarg); 
    205     const real form = REAL(4.0)*bj*bj*si*si; 
     204    const real form = bj*bj*si*si; 
    206205 
    207206    // Multiply by contrast^2, normalize by cylinder volume and convert to cm-1 
    208207    // NOTE that for this (Fournet) definition of the integral, one must MULTIPLY by Vcyl 
    209208    // The additional volume factor is for polydisperse volume normalization. 
    210     const real s = (sld - solvent_sld) * volume(radius, length); 
     209    const real s = (sld - solvent_sld) * form_volume(radius, length); 
    211210    return REAL(1.0e-4) * form * s * s; // * correction; 
    212211} 
Note: See TracChangeset for help on using the changeset viewer.