Changes in sasmodels/models/cylinder.py [5428233:9474dda] in sasmodels
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sasmodels/models/cylinder.py
r5428233 r9474dda 43 43 .. figure:: img/orientation2.jpg 44 44 45 Examples of the angles for oriented ppagainst the detector plane.45 Examples of the angles for oriented cylinders against the detector plane. 46 46 47 47 NB: The 2nd virial coefficient of the cylinder is calculated based on the … … 113 113 title = "Right circular cylinder with uniform scattering length density." 114 114 description = """ 115 P(q)= 2*(sld - solvent_sld)*V*sin(qLcos(alpha/2))116 /[qLcos(alpha/2)]*J1(qRsin(alpha/2))/[qRsin(alpha)]115 f(q,alpha) = 2*(sld - solvent_sld)*V*sin(qLcos(alpha/2)) 116 /[qLcos(alpha/2)]*J1(qRsin(alpha/2))/[qRsin(alpha)] 117 117 118 P(q,alpha)= scale/V*f(q )^(2)+background118 P(q,alpha)= scale/V*f(q,alpha)^(2)+background 119 119 V: Volume of the cylinder 120 120 R: Radius of the cylinder … … 125 125 :the ouput is P(q)=scale/V*integral 126 126 from pi/2 to zero of... 127 f(q )^(2)*sin(alpha)*dalpha + background127 f(q,alpha)^(2)*sin(alpha)*dalpha + background 128 128 """ 129 129 … … 169 169 oldpars=dict(theta='cyl_theta', phi='cyl_phi', sld='sldCyl', solvent_sld='sldSolv') 170 170 171 172 # test values:173 # [174 # [ {parameters}, q, I(q)],175 # [ {parameters}, [q], [I(q)] ],176 # [ {parameters}, [q1, q2, ...], [I(q1), I(q2), ...]],177 178 # [ {parameters}, (qx, qy), I(qx, Iqy)],179 # [ {parameters}, [(qx1, qy1), (qx2, qy2), ...], [I(qx1,qy1), I(qx2,qy2), ...],180 # ...181 # ]182 # Precision defaults to 7 digits (relative) for single, 14 for double183 # May want a limited precision version that checks for 8-n or 15-n digits respectively184 tests = [185 [{}, 0.2, 0.041761386790780453],186 [{'theta':10.0, 'phi':10.0}, [0.2], [0.041761386790780453]],187 [{'theta':10.0, 'phi':10.0}, (0.2, 2.5), 0.038176446608393366],188 [{'theta':10.0, 'phi':10.0}, [(0.2, 2.5)], [0.038176446608393366]],189 ]190
Note: See TracChangeset
for help on using the changeset viewer.