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


Ignore:
Timestamp:
Oct 1, 2016 12:42:06 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:
caddb14, 5031ca3
Parents:
2222134
Message:

updating more parameter names addressing #649

File:
1 edited

Legend:

Unmodified
Added
Removed
  • sasmodels/models/parallelepiped.c

    rdeb7ee0 ra807206  
    1 double form_volume(double a_side, double b_side, double c_side); 
    2 double Iq(double q, double sld, double solvent_sld, double a_side, double b_side, double c_side); 
     1double form_volume(double length_a, double length_b, double length_c); 
     2double Iq(double q, double sld, double solvent_sld, double length_a, double length_b, double length_c); 
    33double Iqxy(double qx, double qy, double sld, double solvent_sld, 
    4     double a_side, double b_side, double c_side, double theta, double phi, double psi); 
     4    double length_a, double length_b, double length_c, double theta, double phi, double psi); 
    55 
    66// From Igor library 
     
    3232 
    3333 
    34 double form_volume(double a_side, double b_side, double c_side) 
     34double form_volume(double length_a, double length_b, double length_c) 
    3535{ 
    36     return a_side * b_side * c_side; 
     36    return length_a * length_b * length_c; 
    3737} 
    3838 
     
    4141    double sld, 
    4242    double solvent_sld, 
    43     double a_side, 
    44     double b_side, 
    45     double c_side) 
     43    double length_a, 
     44    double length_b, 
     45    double length_c) 
    4646{ 
    4747    double tmp1, tmp2; 
    4848     
    49     double mu = q * b_side; 
     49    double mu = q * length_b; 
    5050     
    5151    // Scale sides by B 
    52     double a_scaled = a_side / b_side; 
    53     double c_scaled = c_side / b_side; 
     52    double a_scaled = length_a / length_b; 
     53    double c_scaled = length_c / length_b; 
    5454         
    5555    //Order of integration 
     
    102102    }    
    103103    
    104     const double vd = (sld-solvent_sld) * form_volume(a_side, b_side, c_side); 
     104    const double vd = (sld-solvent_sld) * form_volume(length_a, length_b, length_c); 
    105105     
    106106    // convert from [1e-12 A-1] to [cm-1] and 0.5 factor for outer integral 
     
    113113    double sld, 
    114114    double solvent_sld, 
    115     double a_side, 
    116     double b_side, 
    117     double c_side, 
     115    double length_a, 
     116    double length_b, 
     117    double length_c, 
    118118    double theta, 
    119119    double phi, 
     
    161161     
    162162    // Call the IGOR library function to get the kernel 
    163     double form = _pkernel( q*a_side, q*b_side, q*c_side, cos_val_a, cos_val_b, cos_val_c); 
     163    double form = _pkernel( q*length_a, q*length_b, q*length_c, cos_val_a, cos_val_b, cos_val_c); 
    164164   
    165165    // Multiply by contrast^2 
    166     const double vd = (sld - solvent_sld) * form_volume(a_side, b_side, c_side); 
     166    const double vd = (sld - solvent_sld) * form_volume(length_a, length_b, length_c); 
    167167    return 1.0e-4 * vd * vd * form; 
    168168} 
Note: See TracChangeset for help on using the changeset viewer.