Changeset a807206 in sasmodels for sasmodels/models/stacked_disks.c


Ignore:
Timestamp:
Sep 30, 2016 10:42:06 PM (8 years ago)
Author:
butler
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:
caddb14, 5031ca3
Parents:
2222134
Message:

updating more parameter names addressing #649

File:
1 edited

Legend:

Unmodified
Added
Removed

  • sasmodels/models/stacked_disks.c

    r2c74c11 ra807206  
    1 double form_volume(double core_thick, 
    2                    double layer_thick, 
     1double form_volume(double thick_core, 
     2                   double thick_layer, 
    33                   double radius, 
    44                   double n_stacking); 
    55 
    66double Iq(double q, 
    7           double core_thick, 
    8           double layer_thick, 
     7          double thick_core, 
     8          double thick_layer, 
    99          double radius, 
    1010          double n_stacking, 
    11           double sigma_d, 
     11          double sigma_dnn, 
    1212          double core_sld, 
    1313          double layer_sld, 
    … …  
    2121               double solvent_sld, 
    2222               double halfheight, 
    23                double layer_thick, 
     23               double thick_layer, 
    2424               double zi, 
    25                double sigma_d, 
     25               double sigma_dnn, 
    2626               double d, 
    2727               double n_stacking) 
    … …  
    3838 
    3939        const double sinarg1 = qq*halfheight*cos(zi); 
    40         const double sinarg2 = qq*(halfheight+layer_thick)*cos(zi); 
     40        const double sinarg2 = qq*(halfheight+thick_layer)*cos(zi); 
    4141 
    4242        const double be1 = sas_J1c(besarg1); 
    … …  
    4848        const double dr2 = (layer_sld-solvent_sld); 
    4949        const double area = M_PI*radius*radius; 
    50         const double totald=2.0*(layer_thick+halfheight); 
     50        const double totald=2.0*(thick_layer+halfheight); 
    5151 
    5252        const double t1 = area*(2.0*halfheight)*dr1*(si1)*(be1); 
    … …  
    5959        double sqq=0.0; 
    6060        for(int kk=1;kk<n_stacking;kk+=1) { 
    61                 double dexpt=qq*cos(zi)*qq*cos(zi)*d*d*sigma_d*sigma_d*kk/2.0; 
     61                double dexpt=qq*cos(zi)*qq*cos(zi)*d*d*sigma_dnn*sigma_dnn*kk/2.0; 
    6262                sqq=sqq+(n_stacking-kk)*cos(qq*cos(zi)*d*kk)*exp(-1.*dexpt); 
    6363        } 
    … …  
    7474static 
    7575double stacked_disks_kernel(double q, 
    76                             double core_thick, 
    77                             double layer_thick, 
     76                            double thick_core, 
     77                            double thick_layer, 
    7878                            double radius, 
    7979                            double n_stacking, 
    80                             double sigma_d, 
     80                            double sigma_dnn, 
    8181                            double core_sld, 
    8282                            double layer_sld, 
    … …  
    8888        double summ = 0.0;      //initialize integral 
    8989 
    90         double d=2.0*layer_thick+core_thick; 
    91         double halfheight = core_thick/2.0; 
     90        double d=2.0*thick_layer+thick_core; 
     91        double halfheight = thick_core/2.0; 
    9292 
    9393        for(int i=0;i<N_POINTS_76;i++) { 
    … …  
    100100                                   solvent_sld, 
    101101                                   halfheight, 
    102                                    layer_thick, 
     102                                   thick_layer, 
    103103                                   zi, 
    104                                    sigma_d, 
     104                                   sigma_dnn, 
    105105                                   d, 
    106106                                   n_stacking); 
    … …  
    117117 
    118118static double stacked_disks_kernel_2d(double q, double q_x, double q_y, 
    119                             double core_thick, 
    120                             double layer_thick, 
     119                            double thick_core, 
     120                            double thick_layer, 
    121121                            double radius, 
    122122                            double n_stacking, 
    123                             double sigma_d, 
     123                            double sigma_dnn, 
    124124                            double core_sld, 
    125125                            double layer_sld, 
    … …  
    154154 
    155155    // Call the IGOR library function to get the kernel 
    156     double d = 2 * layer_thick + core_thick; 
    157     double halfheight = core_thick/2.0; 
     156    double d = 2 * thick_layer + thick_core; 
     157    double halfheight = thick_core/2.0; 
    158158    double answer = _kernel(q, 
    159159                     radius, 
    … …  
    162162                     solvent_sld, 
    163163                     halfheight, 
    164                      layer_thick, 
     164                     thick_layer, 
    165165                     alpha, 
    166                      sigma_d, 
     166                     sigma_dnn, 
    167167                     d, 
    168168                     n_stacking); 
    … …  
    175175} 
    176176 
    177 double form_volume(double core_thick, 
    178                    double layer_thick, 
     177double form_volume(double thick_core, 
     178                   double thick_layer, 
    179179                   double radius, 
    180180                   double n_stacking){ 
    181     double d = 2 * layer_thick + core_thick; 
     181    double d = 2 * thick_layer + thick_core; 
    182182    return acos(-1.0) * radius * radius * d * n_stacking; 
    183183} 
    184184 
    185185double Iq(double q, 
    186           double core_thick, 
    187           double layer_thick, 
     186          double thick_core, 
     187          double thick_layer, 
    188188          double radius, 
    189189          double n_stacking, 
    190           double sigma_d, 
     190          double sigma_dnn, 
    191191          double core_sld, 
    192192          double layer_sld, 
    … …  
    194194{ 
    195195    return stacked_disks_kernel(q, 
    196                     core_thick, 
    197                     layer_thick, 
     196                    thick_core, 
     197                    thick_layer, 
    198198                    radius, 
    199199                    n_stacking, 
    200                     sigma_d, 
     200                    sigma_dnn, 
    201201                    core_sld, 
    202202                    layer_sld, 
Note: See TracChangeset for help on using the changeset viewer.