Changeset 2d81cfe in sasmodels for sasmodels/models/core_shell_bicelle.py
- Timestamp:
- Nov 29, 2017 11:13:23 AM (6 years ago)
- Branches:
- master, core_shell_microgels, magnetic_model, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests
- Children:
- 237b800f
- Parents:
- a839b22
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sasmodels/models/core_shell_bicelle.py
r30b60d2 r2d81cfe 50 50 \begin{align*} 51 51 F(Q,\alpha) = &\bigg[ 52 (\rho_c - \rho_f) V_c \frac{2J_1(QRsin \alpha)}{QRsin\alpha}\frac{sin(QLcos\alpha/2)}{Q(L/2)cos\alpha} \\ 53 &+(\rho_f - \rho_r) V_{c+f} \frac{2J_1(QRsin\alpha)}{QRsin\alpha}\frac{sin(Q(L/2+t_f)cos\alpha)}{Q(L/2+t_f)cos\alpha} \\ 54 &+(\rho_r - \rho_s) V_t \frac{2J_1(Q(R+t_r)sin\alpha)}{Q(R+t_r)sin\alpha}\frac{sin(Q(L/2+t_f)cos\alpha)}{Q(L/2+t_f)cos\alpha} 52 (\rho_c - \rho_f) V_c 53 \frac{2J_1(QRsin \alpha)}{QRsin\alpha} 54 \frac{sin(QLcos\alpha/2)}{Q(L/2)cos\alpha} \\ 55 &+(\rho_f - \rho_r) V_{c+f} 56 \frac{2J_1(QRsin\alpha)}{QRsin\alpha} 57 \frac{sin(Q(L/2+t_f)cos\alpha)}{Q(L/2+t_f)cos\alpha} \\ 58 &+(\rho_r - \rho_s) V_t 59 \frac{2J_1(Q(R+t_r)sin\alpha)}{Q(R+t_r)sin\alpha} 60 \frac{sin(Q(L/2+t_f)cos\alpha)}{Q(L/2+t_f)cos\alpha} 55 61 \bigg] 56 62 \end{align*} … … 64 70 cylinders is then given by integrating over all possible $\theta$ and $\phi$. 65 71 66 For oriented bicelles the *theta*, and *phi* orientation parameters will appear when fitting 2D data,67 see the :ref:`cylinder` model for further information.72 For oriented bicelles the *theta*, and *phi* orientation parameters will appear 73 when fitting 2D data, see the :ref:`cylinder` model for further information. 68 74 Our implementation of the scattering kernel and the 1D scattering intensity 69 75 use the c-library from NIST. … … 92 98 """ 93 99 100 import numpy as np 94 101 from numpy import inf, sin, cos, pi 95 102 … … 149 156 150 157 def random(): 151 import numpy as np152 158 pars = dict( 153 159 radius=10**np.random.uniform(1.3, 3), … … 173 179 qx = q*cos(pi/6.0) 174 180 qy = q*sin(pi/6.0) 175 tests = [[{}, 0.05, 7.4883545957], 176 [{'theta':80., 'phi':10.}, (qx, qy), 2.81048892474 ] 177 ] 181 tests = [ 182 [{}, 0.05, 7.4883545957], 183 [{'theta':80., 'phi':10.}, (qx, qy), 2.81048892474] 184 ] 178 185 del qx, qy # not necessary to delete, but cleaner 179
Note: See TracChangeset
for help on using the changeset viewer.