Changeset a151caa in sasmodels for sasmodels/models/core_multi_shell.py

Timestamp:

Jul 28, 2017 10:49:09 PM (7 years ago)

Author:

Paul Kienzle <pkienzle@…>

Branches:

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

Children:

404ebbd

Parents:

0bdddc2

Message:

tuned random model generator for be_polyelectrolyte, barbell, core multishell, core-shell bicelle

File:

• sasmodels/models/core_multi_shell.py (modified) (2 diffs)

sasmodels/models/core_multi_shell.py

@@ -70 +70 @@
         sld_shell: the SLD of the shell#
         A_shell#: the coefficient in the exponential function
-        
-        
+
+
     scale: 1.0 if data is on absolute scale
     volfraction: volume fraction of spheres
@@ -102 +102 @@
 
 source = ["lib/sas_3j1x_x.c", "core_multi_shell.c"]
+
+def random():
+    import numpy as np
+    num_shells = np.minimum(np.random.poisson(3)+1, 10)
+    total_radius = 10**np.random.uniform(1.7, 4)
+    thickness = np.random.exponential(size=num_shells+1)
+    thickness *= total_radius/np.sum(thickness)
+    volume_fraction = 10**np.random.uniform(-4, -1)
+    pars = dict(
+        #background=0,
+        scale=volume_fraction/total_radius**3*1e10,
+        n=num_shells,
+        radius=thickness[0],
+    )
+    for k, v in enumerate(thickness[1:]):
+        pars['thickness%d'%(k+1)] = v
+    return pars
 
 def profile(sld_core, radius, sld_solvent, n, sld, thickness):