Changeset 4e96703 in sasmodels for example

Timestamp:

Nov 9, 2018 12:21:36 PM (5 years ago)

Author:

Paul Kienzle <pkienzle@…>

Branches:

master

Children:

e2da671

Parents:

12f77e9 (diff), cf3d0ce (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the (diff) links above to see all the changes relative to each parent.

Message:

Merge branch 'beta_approx' into ticket-608-user-defined-weights

Location:

example

Files:

• cylinder_eval.py (modified) (1 diff)
• weights/cyclic_gaussian.py (added)
• weights/gaussian_eq.py (added)
• weights/laplace.py (added)
• weights/maier_saupe.py (added)
• weights/maier_saupe_eq.py (added)

example/cylinder_eval.py

@@ -3 +3 @@
 """
 
-from numpy import logspace
+from numpy import logspace, sqrt
 from matplotlib import pyplot as plt
 from sasmodels.core import load_model
-from sasmodels.direct_model import call_kernel
+from sasmodels.direct_model import call_kernel, call_Fq
 
 model = load_model('cylinder')
 q = logspace(-3, -1, 200)
 kernel = model.make_kernel([q])
-Iq = call_kernel(kernel, dict(radius=200.))
-plt.loglog(q, Iq)
+pars = {'radius': 200, 'radius_pd': 0.1, 'scale': 2}
+Iq = call_kernel(kernel, pars)
+F, Fsq, Reff, V, Vratio = call_Fq(kernel, pars)
+plt.loglog(q, Iq, label='2 I(q)')
+plt.loglog(q, F**2/V, label='<F(q)>^2/V')
+plt.loglog(q, Fsq/V, label='<F^2(q)>/V')
 plt.xlabel('q (1/A)')
-plt.ylabel('I(q)')
+plt.ylabel('I(q) (1/cm)')
 plt.title('Cylinder with radius 200.')
+plt.legend()
 plt.show()