Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • sasmodels/models/elliptical_cylinder.c

    rd277229 rfbaef04  
    1313 
    1414static double 
    15 radius_from_min_dimension(double radius_minor, double r_ratio, double length) 
     15radius_from_min_dimension(double radius_minor, double r_ratio, double hlength) 
    1616{ 
    1717    const double rad_min = (r_ratio > 1.0 ? radius_minor : r_ratio*radius_minor); 
    18     return (rad_min < length ? rad_min : length); 
     18    return (rad_min < hlength ? rad_min : hlength); 
    1919} 
    2020 
    2121static double 
    22 radius_from_max_dimension(double radius_minor, double r_ratio, double length) 
     22radius_from_max_dimension(double radius_minor, double r_ratio, double hlength) 
    2323{ 
    2424    const double rad_max = (r_ratio < 1.0 ? radius_minor : r_ratio*radius_minor); 
    25     return (rad_max > length ? rad_max : length); 
     25    return (rad_max > hlength ? rad_max : hlength); 
    2626} 
    2727 
     
    3535static double 
    3636effective_radius(int mode, double radius_minor, double r_ratio, double length) 
     37//effective_radius_type = ["equivalent sphere","average radius","min radius","max radius", 
     38//                         "equivalent circular cross-section","half length","half min dimension","half max dimension","half diagonal"] 
    3739{ 
    3840    if (mode == 1) { 
     
    4951        return 0.5*length; 
    5052    } else if (mode == 7) { 
    51         return radius_from_min_dimension(radius_minor,r_ratio,length); 
     53        return radius_from_min_dimension(radius_minor,r_ratio,0.5*length); 
    5254    } else if (mode == 8) { 
    53         return radius_from_max_dimension(radius_minor,r_ratio,length); 
     55        return radius_from_max_dimension(radius_minor,r_ratio,0.5*length); 
    5456    } else { 
    5557        return radius_from_diagonal(radius_minor,r_ratio,length); 
Note: See TracChangeset for help on using the changeset viewer.