Changes in / [8b9adb5:c9f885e] in sasmodels
- Location:
- sasmodels/models
- Files:
-
- 2 deleted
- 17 edited
Legend:
- Unmodified
- Added
- Removed
-
sasmodels/models/be_polyelectrolyte.py
r0e86967 r168052c 127 127 :return: 2D-Intensity 128 128 """ 129 i ntensity= Iq(sqrt(qx**2 + qy**2), *args)130 return i ntensity129 iq = Iq(sqrt(qx**2 + qy**2), *args) 130 return iq 131 131 132 132 Iqxy.vectorized = True # Iqxy accepts an array of qx, qy values -
sasmodels/models/core_shell_cylinder.py
rf0aa7f8 reb69cce 155 155 156 156 def ER(radius, thickness, length): 157 """158 Returns the effective radius used in the S*P calculation159 """160 157 radius = radius + thickness 161 158 length = length + 2 * thickness … … 164 161 165 162 def VR(radius, thickness, length): 166 """167 Returns volume ratio168 """169 163 whole = pi * (radius + thickness) ** 2 * (length + 2 * thickness) 170 164 core = pi * radius ** 2 * length -
sasmodels/models/correlation_length.py
r3e6c5c1 r5054e80 78 78 # names and the target sasview model name. 79 79 oldname = 'CorrLengthModel' 80 # pylint: disable=bad-continuation 81 oldpars = dict( 82 lorentz_scale='scale_l', porod_scale='scale_p', 83 cor_length='length_l', exponent_p='exponent_p', 84 exponent_l='exponent_l' 85 ) 80 86 81 oldpars = dict(lorentz_scale='scale_l', porod_scale='scale_p', 82 cor_length='length_l', exponent_p='exponent_p', 83 exponent_l='exponent_l') 84 85 tests = [[{}, 0.001, 1009.98], 87 tests = [ 88 [{}, 0.001, 1009.98], 86 89 [{}, 0.150141, 0.174645], 87 [{}, 0.442528, 0.0203957]] 90 [{}, 0.442528, 0.0203957] 91 ] 92 # pylint: enable=bad-continuation -
sasmodels/models/cylinder.py
r0c3a226 reb69cce 144 144 145 145 def ER(radius, length): 146 """147 Return equivalent radius (ER)148 """149 146 ddd = 0.75 * radius * (2 * radius * length + (length + radius) * (length + pi * radius)) 150 147 return 0.5 * (ddd) ** (1. / 3.) -
sasmodels/models/dab.py
r94bd809 reb69cce 1 1 r""" 2 2 3 Calculates the scattering from a randomly distributed, two-phase system based on 3 4 the Debye-Anderson-Brumberger (DAB) model for such systems. The two-phase system -
sasmodels/models/fcc.py
rc0ccea8 r9aac25d 114 114 single = False 115 115 116 # pylint: disable=bad-whitespace, line-too-long117 116 # ["name", "units", default, [lower, upper], "type","description"], 118 117 parameters = [["dnn", "Ang", 220, [-inf, inf], "", "Nearest neighbour distance"], … … 125 124 ["psi", "degrees", 60, [-inf, inf], "orientation", "Out of plane angle"] 126 125 ] 127 # pylint: enable=bad-whitespace, line-too-long128 126 129 127 source = ["lib/J1.c", "lib/gauss150.c", "lib/sphere_form.c", "fcc.c"] -
sasmodels/models/guinier.py
r723cebe reb69cce 36 36 37 37 # ["name", "units", default, [lower, upper], "type","description"], 38 parameters = [["rg", "Ang", 60.0, [0, inf], "", "Radius of Gyration"]] 38 parameters = [ 39 ["rg", "Ang", 60.0, [0, inf], "", "Radius of Gyration"], 40 ] 39 41 40 42 Iq = """ … … 49 51 50 52 # parameters for demo 51 demo = dict(scale=1.0, 53 demo = dict(scale=1.0,rg=60.0) 52 54 53 55 # For testing against the old sasview models, include the converted parameter … … 57 59 58 60 # parameters for unit tests 59 tests = [[{'rg' : 31.5}, 0.005, 0.991756]] 61 tests = [ 62 [{'rg' : 31.5}, 0.005, 0.991756] 63 ] -
sasmodels/models/hollow_cylinder.py
re0fd913 rec2ca99 93 93 source = ["lib/J1.c", "lib/gauss76.c", "hollow_cylinder.c"] 94 94 95 # pylint: disable=W061396 95 def ER(radius, core_radius, length): 97 96 """ -
sasmodels/models/lamellarCailleHG.py
r652a78a r13ed84c 146 146 solvent_sld='sld_solvent') 147 147 # 148 tests = [[{'scale': 1.0, 'background': 0.0, 'tail_length': 10.0, 'head_length': 2.0, 149 'Nlayers': 30.0, 'spacing': 40., 'Caille_parameter': 0.001, 'sld': 0.4, 150 'head_sld': 2.0, 'solvent_sld': 6.0, 'tail_length_pd': 0.0, 151 'head_length_pd': 0.0, 'spacing_pd': 0.0}, [0.001], [6838238.571488]]] 148 tests = [ 149 [ {'scale': 1.0, 'background' : 0.0, 'tail_length' : 10.0, 'head_length' : 2.0,'Nlayers' : 30.0, 'spacing' : 40., 150 'Caille_parameter' : 0.001, 'sld' : 0.4, 'head_sld' : 2.0, 'solvent_sld' : 6.0, 151 'tail_length_pd' : 0.0, 'head_length_pd' : 0.0, 'spacing_pd' : 0.0 }, [0.001], [6838238.571488]] 152 ] -
sasmodels/models/lamellarFFHG.py
r3eb6b90 r7f47777 68 68 category = "shape:lamellae" 69 69 70 # pylint: disable=bad-whitespace, line-too-long71 70 # ["name", "units", default, [lower, upper], "type","description"], 72 parameters = [["tail_length", "Ang", 15, [0, inf], "volume", "Tail thickness"], 73 ["head_length", "Ang", 10, [0, inf], "volume", "head thickness"], 74 ["sld", "1e-6/Ang^2", 0.4, [-inf,inf], "", "Tail scattering length density"], 75 ["head_sld", "1e-6/Ang^2", 3.0, [-inf,inf], "", "Head scattering length density"], 76 ["solvent_sld", "1e-6/Ang^2", 6, [-inf,inf], "", "Solvent scattering length density"]] 77 # pylint: enable=bad-whitespace, line-too-long 71 parameters = [["tail_length", "Ang", 15, [0, inf], "volume", 72 "Tail thickness"], 73 ["head_length", "Ang", 10, [0, inf], "volume", 74 "head thickness"], 75 ["sld", "1e-6/Ang^2", 0.4, [-inf,inf], "", 76 "Tail scattering length density"], 77 ["head_sld", "1e-6/Ang^2", 3.0, [-inf,inf], "", 78 "Head scattering length density"], 79 ["solvent_sld", "1e-6/Ang^2", 6, [-inf,inf], "", 80 "Solvent scattering length density"], 81 ] 78 82 79 83 # No volume normalization despite having a volume parameter … … 109 113 110 114 demo = dict(scale=1, background=0, 111 tail_length=15, 115 tail_length=15,head_length=10, 112 116 sld=0.4, head_sld=3.0, solvent_sld=6.0, 113 tail_length_pd= 0.2, tail_length_pd_n=40,114 head_length_pd= 0.01, head_length_pd_n=40)117 tail_length_pd= 0.2, tail_length_pd_n=40, 118 head_length_pd= 0.01, head_length_pd_n=40) 115 119 116 120 oldname = 'LamellarFFHGModel' … … 118 122 sld='sld_tail', head_sld='sld_head', solvent_sld='sld_solvent') 119 123 # 120 tests = [[{'scale': 1.0, 'background': 0.0, 'tail_length': 15.0, 'head_length': 10.0, 121 'sld': 0.4, 'head_sld': 3.0, 'solvent_sld': 6.0}, [0.001], [653143.9209]]] 124 tests = [ 125 [ {'scale': 1.0, 'background' : 0.0, 'tail_length' : 15.0, 'head_length' : 10.0,'sld' : 0.4, 126 'head_sld' : 3.0, 'solvent_sld' : 6.0, }, [0.001], [653143.9209]] 127 ] 128 129 -
sasmodels/models/lorentz.py
reb5901b reb69cce 56 56 57 57 # parameters for demo 58 demo = dict(scale=1.0, background=0.0,cor_length=50.0)58 demo = dict(scale=1.0,background=0.0,cor_length=50.0) 59 59 60 60 # For testing against the old sasview models, include the converted parameter … … 64 64 65 65 # parameters for unit tests 66 tests = [[{'cor_length': 250}, 0.01, 0.137931]] 66 tests = [ 67 [{'cor_length' : 250},0.01,0.137931] 68 ] -
sasmodels/models/peak_lorentz.py
r04b0b30 r14ba6f6 11 11 I(q) = \frac{scale}{\bigl(1+\bigl(\frac{q-q_0}{B}\bigr)^2\bigr)} + background 12 12 13 with the peak having height of $I_0$ centered at $q_0$ and having 14 a HWHM (half-width half-maximum) of B. 13 with the peak having height of $I_0$ centered at $q_0$ and having a HWHM (half-width half-maximum) of B. 15 14 16 15 For 2D data the scattering intensity is calculated in the same way as 1D, … … 56 55 57 56 def Iq(q, peak_pos, peak_hwhm): 58 """59 Return I(q)60 """61 57 inten = (1/(1+((q-peak_pos)/peak_hwhm)**2)) 62 58 return inten … … 64 60 65 61 def Iqxy(qx, qy, *args): 66 """67 Return I(qx, qy)68 """69 62 return Iq(sqrt(qx ** 2 + qy ** 2), *args) 70 63 Iqxy.vectorized = True # Iqxy accepts an array of qx, qy values -
sasmodels/models/porod.py
r1ce6e82 r16bb433 13 13 of the sample, and $\Delta\rho$ is the contrast factor. 14 14 15 For 2D data: The 2D scattering intensity is calculated in the same way as 1D, 15 For 2D data: The 2D scattering intensity is calculated in the same way as 1D, 16 16 where the q vector is defined as 17 17 -
sasmodels/models/rpa.py
r8dd6914 r82c299f 89 89 # ["name", "units", default, [lower, upper], "type","description"], 90 90 parameters = [ 91 ["case_num", CASES, 0, [0, 10], "", "Component organization" ],91 ["case_num", CASES, 0, [0, 10], "", "Component organization" ], 92 92 93 93 ["Na", "", 1000.0, [1, inf], "", "Degree of polymerization"], -
sasmodels/models/sphere.py
rc691551 r9c461c7 100 100 101 101 def ER(radius): 102 """103 Return equivalent radius (ER)104 """105 102 return radius 106 103 -
sasmodels/models/triaxial_ellipsoid.py
r469e763 r9c461c7 82 82 ---------- 83 83 84 L A Feigin and D I Svergun, *Structure Analysis by Small-Angle X-Ray 85 and Neutron Scattering*, Plenum,New York, 1987.84 L A Feigin and D I Svergun, *Structure Analysis by Small-Angle X-Ray and Neutron Scattering*, Plenum, 85 New York, 1987. 86 86 """ 87 87 … … 120 120 121 121 def ER(req_minor, req_major, rpolar): 122 """123 Returns the effective radius used in the S*P calculation124 """125 122 import numpy as np 126 123 from .ellipsoid import ER as ellipsoid_ER -
sasmodels/models/vesicle.py
r068cebd r216fa6d 60 60 """ 61 61 62 import numpy as np 62 63 from numpy import pi, inf 63 64 … … 130 131 # NOTE: test results taken from values returned by SasView 3.1.2 131 132 tests = [[{}, 0.0010005303255, 17139.8268799], 132 [{}, 0.200027832249, 0.130387268704 ],133 [{}, 0.200027832249, 0.130387268704 ], 133 134 [{}, 'ER', 130.], 134 135 [{}, 'VR', 0.54483386436],
Note: See TracChangeset
for help on using the changeset viewer.