Changeset 2c74c11 in sasmodels for sasmodels/models/stacked_disks.c


Ignore:
Timestamp:
Jul 24, 2016 10:56:45 PM (8 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:
a4280bd
Parents:
f1765a2
Message:

implicit Iqxy; fix divide by 0 for q=0

File:
1 edited

Legend:

Unmodified
Added
Removed
  • sasmodels/models/stacked_disks.c

    r43b7eea r2c74c11  
    1313          double layer_sld, 
    1414          double solvent_sld); 
    15  
    16 double Iqxy(double qx, double qy, 
    17             double core_thick, 
    18             double layer_thick, 
    19             double radius, 
    20             double n_stacking, 
    21             double sigma_d, 
    22             double core_sld, 
    23             double layer_sld, 
    24             double solvent_sld, 
    25             double theta, 
    26             double phi); 
    2715 
    2816static 
     
    5240        const double sinarg2 = qq*(halfheight+layer_thick)*cos(zi); 
    5341 
    54     const double be1 = sas_J1c(besarg1); 
     42        const double be1 = sas_J1c(besarg1); 
    5543        const double be2 = sas_J1c(besarg2); 
    5644        const double si1 = sin(sinarg1)/sinarg1; 
     
    215203                    solvent_sld); 
    216204} 
    217  
    218 // Iqxy is never called since no orientation or magnetic parameters. 
    219 double Iqxy(double qx, double qy, 
    220             double core_thick, 
    221             double layer_thick, 
    222             double radius, 
    223             double n_stacking, 
    224             double sigma_d, 
    225             double core_sld, 
    226             double layer_sld, 
    227             double solvent_sld, 
    228             double theta, 
    229             double phi) 
    230 { 
    231     double q = sqrt(qx*qx + qy*qy); 
    232     return stacked_disks_kernel_2d(q, qx/q, qy/q, 
    233                     core_thick, 
    234                     layer_thick, 
    235                     radius, 
    236                     n_stacking, 
    237                     sigma_d, 
    238                     core_sld, 
    239                     layer_sld, 
    240                     solvent_sld, 
    241                     theta, 
    242                     phi); 
    243 } 
    244  
Note: See TracChangeset for help on using the changeset viewer.