Ignore:
Timestamp:
Sep 30, 2016 9:07:16 AM (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:
a807206
Parents:
6e5b2a7
Message:

Updating parameter names regarding #649

File:
1 edited

Legend:

Unmodified
Added
Removed

  • sasmodels/models/core_shell_ellipsoid.c

    r29172aa r2222134  
    1 double form_volume(double equat_core, 
    2                    double polar_core, 
    3                    double equat_shell, 
    4                    double polar_shell); 
     1double form_volume(double radius_equat_core, 
     2                   double radius_polar_core, 
     3                   double radius_equat_shell, 
     4                   double radius_polar_shell); 
    55double Iq(double q, 
    6           double equat_core, 
    7           double polar_core, 
    8           double equat_shell, 
    9           double polar_shell, 
     6          double radius_equat_core, 
     7          double radius_polar_core, 
     8          double radius_equat_shell, 
     9          double radius_polar_shell, 
    1010          double sld_core, 
    1111          double sld_shell, 
    … …  
    1414 
    1515double Iqxy(double qx, double qy, 
    16           double equat_core, 
    17           double polar_core, 
    18           double equat_shell, 
    19           double polar_shell, 
     16          double radius_equat_core, 
     17          double radius_polar_core, 
     18          double radius_equat_shell, 
     19          double radius_polar_shell, 
    2020          double sld_core, 
    2121          double sld_shell, 
    … …  
    2525 
    2626 
    27 double form_volume(double equat_core, 
    28                    double polar_core, 
    29                    double equat_shell, 
    30                    double polar_shell) 
     27double form_volume(double radius_equat_core, 
     28                   double radius_polar_core, 
     29                   double radius_equat_shell, 
     30                   double radius_polar_shell) 
    3131{ 
    32     double vol = 4.0*M_PI/3.0*equat_shell*equat_shell*polar_shell; 
     32    double vol = 4.0*M_PI/3.0*radius_equat_shell*radius_equat_shell*radius_polar_shell; 
    3333    return vol; 
    3434} 
    … …  
    3636static double 
    3737core_shell_ellipsoid_kernel(double q, 
    38           double equat_core, 
    39           double polar_core, 
    40           double equat_shell, 
    41           double polar_shell, 
     38          double radius_equat_core, 
     39          double radius_polar_core, 
     40          double radius_equat_shell, 
     41          double radius_polar_shell, 
    4242          double sld_core, 
    4343          double sld_shell, 
    … …  
    5757        double zi = ( Gauss76Z[i]*(uplim-lolim) + uplim + lolim )/2.0; 
    5858        double yyy = Gauss76Wt[i] * gfn4(zi, 
    59                                   equat_core, 
    60                                   polar_core, 
    61                                   equat_shell, 
    62                                   polar_shell, 
     59                                  radius_equat_core, 
     60                                  radius_polar_core, 
     61                                  radius_equat_shell, 
     62                                  radius_polar_shell, 
    6363                                  delpc, 
    6464                                  delps, 
    … …  
    7777static double 
    7878core_shell_ellipsoid_kernel_2d(double q, double q_x, double q_y, 
    79           double equat_core, 
    80           double polar_core, 
    81           double equat_shell, 
    82           double polar_shell, 
     79          double radius_equat_core, 
     80          double radius_polar_core, 
     81          double radius_equat_shell, 
     82          double radius_polar_shell, 
    8383          double sld_core, 
    8484          double sld_shell, 
    … …  
    105105    // Call the IGOR library function to get the kernel: MUST use gfn4 not gf2 because of the def of params. 
    106106    double answer = gfn4(cos_val, 
    107                   equat_core, 
    108                   polar_core, 
    109                   equat_shell, 
    110                   polar_shell, 
     107                  radius_equat_core, 
     108                  radius_polar_core, 
     109                  radius_equat_shell, 
     110                  radius_polar_shell, 
    111111                  sldcs, 
    112112                  sldss, 
    … …  
    120120 
    121121double Iq(double q, 
    122           double equat_core, 
    123           double polar_core, 
    124           double equat_shell, 
    125           double polar_shell, 
     122          double radius_equat_core, 
     123          double radius_polar_core, 
     124          double radius_equat_shell, 
     125          double radius_polar_shell, 
    126126          double sld_core, 
    127127          double sld_shell, 
    … …  
    129129{ 
    130130    double intensity = core_shell_ellipsoid_kernel(q, 
    131            equat_core, 
    132            polar_core, 
    133            equat_shell, 
    134            polar_shell, 
     131           radius_equat_core, 
     132           radius_polar_core, 
     133           radius_equat_shell, 
     134           radius_polar_shell, 
    135135           sld_core, 
    136136           sld_shell, 
    … …  
    142142 
    143143double Iqxy(double qx, double qy, 
    144           double equat_core, 
    145           double polar_core, 
    146           double equat_shell, 
    147           double polar_shell, 
     144          double radius_equat_core, 
     145          double radius_polar_core, 
     146          double radius_equat_shell, 
     147          double radius_polar_shell, 
    148148          double sld_core, 
    149149          double sld_shell, 
    … …  
    155155    q = sqrt(qx*qx+qy*qy); 
    156156    double intensity = core_shell_ellipsoid_kernel_2d(q, qx/q, qy/q, 
    157                        equat_core, 
    158                        polar_core, 
    159                        equat_shell, 
    160                        polar_shell, 
     157                       radius_equat_core, 
     158                       radius_polar_core, 
     159                       radius_equat_shell, 
     160                       radius_polar_shell, 
    161161                       sld_core, 
    162162                       sld_shell, 
Note: See TracChangeset for help on using the changeset viewer.