Changeset a807206 in sasmodels for sasmodels/models/parallelepiped.py
- Timestamp:
- Sep 30, 2016 10:42:06 PM (8 years ago)
- 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
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sasmodels/models/parallelepiped.py
r9a4811a ra807206 190 190 ["sld_solvent", "1e-6/Ang^2", 1, [-inf, inf], "sld", 191 191 "Solvent scattering length density"], 192 [" a_side", "Ang", 35, [0, inf], "volume",192 ["length_a", "Ang", 35, [0, inf], "volume", 193 193 "Shorter side of the parallelepiped"], 194 [" b_side", "Ang", 75, [0, inf], "volume",194 ["length_b", "Ang", 75, [0, inf], "volume", 195 195 "Second side of the parallelepiped"], 196 [" c_side", "Ang", 400, [0, inf], "volume",196 ["length_c", "Ang", 400, [0, inf], "volume", 197 197 "Larger side of the parallelepiped"], 198 198 ["theta", "degrees", 60, [-inf, inf], "orientation", … … 206 206 source = ["lib/gauss76.c", "parallelepiped.c"] 207 207 208 def ER( a_side, b_side, c_side):208 def ER(length_a, length_b, length_c): 209 209 """ 210 210 Return effective radius (ER) for P(q)*S(q) … … 212 212 213 213 # 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)) 217 217 return 0.5 * (ddd) ** (1. / 3.) 218 218 … … 222 222 demo = dict(scale=1, background=0, 223 223 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, 225 225 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, 229 229 theta_pd=10, theta_pd_n=1, 230 230 phi_pd=10, phi_pd_n=1,
Note: See TracChangeset
for help on using the changeset viewer.