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


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/parallelepiped.py

    r9a4811a ra807206  
    190190              ["sld_solvent", "1e-6/Ang^2", 1, [-inf, inf], "sld", 
    191191               "Solvent scattering length density"], 
    192               ["a_side", "Ang", 35, [0, inf], "volume", 
     192              ["length_a", "Ang", 35, [0, inf], "volume", 
    193193               "Shorter side of the parallelepiped"], 
    194               ["b_side", "Ang", 75, [0, inf], "volume", 
     194              ["length_b", "Ang", 75, [0, inf], "volume", 
    195195               "Second side of the parallelepiped"], 
    196               ["c_side", "Ang", 400, [0, inf], "volume", 
     196              ["length_c", "Ang", 400, [0, inf], "volume", 
    197197               "Larger side of the parallelepiped"], 
    198198              ["theta", "degrees", 60, [-inf, inf], "orientation", 
     
    206206source = ["lib/gauss76.c", "parallelepiped.c"] 
    207207 
    208 def ER(a_side, b_side, c_side): 
     208def ER(length_a, length_b, length_c): 
    209209    """ 
    210210        Return effective radius (ER) for P(q)*S(q) 
     
    212212 
    213213    # surface average radius (rough approximation) 
    214     surf_rad = sqrt(a_side * b_side / pi) 
    215  
    216     ddd = 0.75 * surf_rad * (2 * surf_rad * c_side + (c_side + surf_rad) * (c_side + pi * surf_rad)) 
     214    surf_rad = sqrt(length_a * length_b / pi) 
     215 
     216    ddd = 0.75 * surf_rad * (2 * surf_rad * length_c + (length_c + surf_rad) * (length_c + pi * surf_rad)) 
    217217    return 0.5 * (ddd) ** (1. / 3.) 
    218218 
     
    222222demo = dict(scale=1, background=0, 
    223223            sld=6.3e-6, sld_solvent=1.0e-6, 
    224             a_side=35, b_side=75, c_side=400, 
     224            length_a=35, length_b=75, length_c=400, 
    225225            theta=45, phi=30, psi=15, 
    226             a_side_pd=0.1, a_side_pd_n=10, 
    227             b_side_pd=0.1, b_side_pd_n=1, 
    228             c_side_pd=0.1, c_side_pd_n=1, 
     226            length_a_pd=0.1, length_a_pd_n=10, 
     227            length_b_pd=0.1, length_b_pd_n=1, 
     228            length_c_pd=0.1, length_c_pd_n=1, 
    229229            theta_pd=10, theta_pd_n=1, 
    230230            phi_pd=10, phi_pd_n=1, 
Note: See TracChangeset for help on using the changeset viewer.