Changeset e5d7a60 in sasmodels

Timestamp:

Jun 29, 2017 1:14:24 AM (7 years ago)

Author:

Paul Kienzle <pkienzle@…>

Branches:

master, core_shell_microgels, magnetic_model, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests

Children:

a2fcbd8

Parents:

4f611f1

Message:

play with cylinder integral

File:

• explore/symint.py (modified) (2 diffs)

explore/symint.py

@@ -103 +103 @@
     Zq = kernel(q=q, theta=theta)
     Zq *= abs(sin(theta))
-    pylab.semilogy(degrees(theta), np.fmax(Zq, 1.e-6))
-    pylab.title("%s Z(q) for q=%g" % (KERNEL.__name__, q))
+    pylab.semilogy(degrees(theta), np.fmax(Zq, 1.e-6), label="Q=%g"%q)
+    pylab.title("%s I(q, theta) sin(theta)" % (KERNEL.__name__,))
     pylab.xlabel("theta (degrees)")
     pylab.ylabel("Iq 1/cm")
-    pylab.show()
 
 NORM, KERNEL = make_cylinder(radius=10., length=100000.)
@@ -125 +124 @@
     gridded_integrals(Q, n=2**19+1)
     #scipy_romberg(Q)
-    plot(Q, n=2000)
+    plot(0.5, n=2000)
+    plot(0.6, n=2000)
+    plot(0.8, n=2000)
+    pylab.legend()
+    pylab.show()