Ignore:
Timestamp:
Sep 10, 2018 4:16:46 PM (6 years ago)
Author:
Paul Kienzle <pkienzle@…>
Branches:
master, core_shell_microgels, magnetic_model, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests
Children:
d299327
Parents:
3f818b2
Message:

clean up effective radius functions; improve mono_gauss_coil accuracy; start moving VR into C

File:
1 edited

Legend:

Unmodified
Added
Removed
  • sasmodels/models/hollow_rectangular_prism_thin_walls.py

    rc44b611 ree60aa7  
    109109source = ["lib/gauss76.c", "hollow_rectangular_prism_thin_walls.c"] 
    110110have_Fq = True 
    111 effective_radius_type = ["equivalent sphere","half length_a", "half length_b", "half length_c", 
    112                          "equivalent outer circular cross-section","half ab diagonal","half diagonal"] 
    113  
    114 #def ER(length_a, b2a_ratio, c2a_ratio): 
    115 #    """ 
    116 #        Return equivalent radius (ER) 
    117 #    """ 
    118 #    b_side = length_a * b2a_ratio 
    119 #    c_side = length_a * c2a_ratio 
    120 # 
    121 #    # surface average radius (rough approximation) 
    122 #    surf_rad = sqrt(length_a * b_side / pi) 
    123 # 
    124 #    ddd = 0.75 * surf_rad * (2 * surf_rad * c_side + (c_side + surf_rad) * (c_side + pi * surf_rad)) 
    125 #    return 0.5 * (ddd) ** (1. / 3.) 
    126  
    127 def VR(length_a, b2a_ratio, c2a_ratio): 
    128     """ 
    129         Return shell volume and total volume 
    130     """ 
    131     b_side = length_a * b2a_ratio 
    132     c_side = length_a * c2a_ratio 
    133     vol_total = length_a * b_side * c_side 
    134     vol_shell = 2.0 * (length_a*b_side + length_a*c_side + b_side*c_side) 
    135     return vol_shell, vol_total 
     111effective_radius_type = [ 
     112    "equivalent sphere", "half length_a", "half length_b", "half length_c", 
     113    "equivalent outer circular cross-section", 
     114    "half ab diagonal", "half diagonal", 
     115    ] 
    136116 
    137117 
Note: See TracChangeset for help on using the changeset viewer.