Changeset e187b25 in sasmodels for sasmodels/models/onion.py

Timestamp:

Aug 5, 2016 8:48:52 AM (8 years ago)

Author:

Paul Kienzle <pkienzle@…>

Branches:

master, core_shell_microgels, costrafo411, magnetic_model, release_v0.94, release_v0.95, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests

Children:

785cbec

Parents:

50ec515

Message:

onion, core multi shell: fix profile plots

File:

• sasmodels/models/onion.py (modified) (4 diffs)

sasmodels/models/onion.py

@@ -306 +306 @@
 parameters = [["sld_core", "1e-6/Ang^2", 1.0, [-inf, inf], "sld",
                "Core scattering length density"],
-              ["core_radius", "Ang", 200., [0, inf], "volume",
+              ["radius_core", "Ang", 200., [0, inf], "volume",
                "Radius of the core"],
               ["sld_solvent", "1e-6/Ang^2", 6.4, [-inf, inf], "sld",
@@ -325 +325 @@
 single = False
 
-#def Iq(q, *args, **kw):
-#    return q
-
 profile_axes = ['Radius (A)', 'SLD (1e-6/A^2)']
-def profile(core_sld, core_radius, solvent_sld, n_shells,
-            in_sld, out_sld, thickness, A):
+def profile(sld_core, radius_core, sld_solvent, n_shells,
+            sld_in, sld_out, thickness, A):
     """
     Returns shape profile with x=radius, y=SLD.
     """
 
-    total_radius = 1.25*(sum(thickness[:n_shells]) + core_radius + 1)
+    total_radius = 1.25*(sum(thickness[:n_shells]) + radius_core + 1)
     dr = total_radius/400  # 400 points for a smooth plot
 
     r = []
-    beta = []
+    rho = []
 
     # add in the core
     r.append(0)
-    beta.append(core_sld)
-    r.append(core_radius)
-    beta.append(core_sld)
+    rho.append(sld_core)
+    r.append(radius_core)
+    rho.append(sld_core)
 
     # add in the shells
@@ -352 +349 @@
         r0 = r[-1]
         r.append(r0)
-        beta.append(in_sld[k])
+        rho.append(sld_in[k])
 
         if fabs(A[k]) < 1.0e-16:
             # flat shell
             r.append(r0 + thickness[k])
-            beta.append(out_sld[k])
+            rho.append(sld_out[k])
         else:
             # exponential shell
@@ -363 +360 @@
             # to protect against a thickness0.
             num_steps = np.floor(thickness[k]/dr) + 1
-            slope = (out_sld[k] - in_sld[k])/expm1(A[k])
-            const = (in_sld[k] - slope)
+            slope = (sld_out[k] - sld_in[k]) / expm1(A[k])
+            const = (sld_in[k] - slope)
             for rk in np.linspace(0, thickness[k], num_steps+1):
                 r.append(r0+rk)
-                beta.append(slope*exp(A[k]*rk/thickness[k]) + const)
+                rho.append(slope*exp(A[k]*rk/thickness[k]) + const)
 
     # add in the solvent
     r.append(r[-1])
-    beta.append(solvent_sld)
+    rho.append(sld_solvent)
     r.append(total_radius)
-    beta.append(solvent_sld)
-
-    return np.asarray(r), np.asarray(beta)*1e-6
+    rho.append(sld_solvent)
+
+    return np.asarray(r), np.asarray(rho)
 
 def ER(core_radius, n, thickness):