Changeset 99658f6 in sasmodels for sasmodels/models/pearl_necklace.py

Timestamp:

Nov 6, 2018 2:10:43 PM (5 years ago)

Author:

grethevj

Branches:

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

Children:

cf3d0ce

Parents:

5024a56

Message:

updated ER functions including cylinder excluded volume, to match 4.x

File:

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

sasmodels/models/pearl_necklace.py

@@ -53 +53 @@
 R Schweins and K Huber, *Particle Scattering Factor of Pearl Necklace Chains*,
 *Macromol. Symp.* 211 (2004) 25-42 2004
+L. Onsager, Ann. New York Acad. Sci. 51, 627-659 (1949). 
 """
 
@@ -95 +96 @@
 source = ["lib/sas_Si.c", "lib/sas_3j1x_x.c", "pearl_necklace.c"]
 single = False  # use double precision unless told otherwise
-
-def volume(radius, edge_sep, thick_string, num_pearls):
-    """
-    Calculates the total particle volume of the necklace.
-    Redundant with form_volume.
-    """
-    num_pearls = int(num_pearls + 0.5)
-    number_of_strings = num_pearls - 1.0
-    string_vol = edge_sep * pi * pow((thick_string / 2.0), 2.0)
-    pearl_vol = 4.0 /3.0 * pi * pow(radius, 3.0)
-    total_vol = number_of_strings * string_vol
-    total_vol += num_pearls * pearl_vol
-    return total_vol
-
-def ER(radius, edge_sep, thick_string, num_pearls):
-    """
-    Calculation for effective radius.
-    """
-    num_pearls = int(num_pearls + 0.5)
-    tot_vol = volume(radius, edge_sep, thick_string, num_pearls)
-    rad_out = (tot_vol/(4.0/3.0*pi)) ** (1./3.)
-    return rad_out
+effective_radius_type = [
+    "equivalent volume sphere", 
+    ]
+    
+#def volume(radius, edge_sep, thick_string, num_pearls):
+#    """
+#    Calculates the total particle volume of the necklace.
+#    Redundant with form_volume.
+#    """
+#    num_pearls = int(num_pearls + 0.5)
+#    number_of_strings = num_pearls - 1.0
+#    string_vol = edge_sep * pi * pow((thick_string / 2.0), 2.0)
+#    pearl_vol = 4.0 /3.0 * pi * pow(radius, 3.0)
+#    total_vol = number_of_strings * string_vol
+#    total_vol += num_pearls * pearl_vol
+#    return total_vol
+#
+#def ER(radius, edge_sep, thick_string, num_pearls):
+#    """
+#    Calculation for effective radius.
+#    """
+#    num_pearls = int(num_pearls + 0.5)
+#    tot_vol = volume(radius, edge_sep, thick_string, num_pearls)
+#    rad_out = (tot_vol/(4.0/3.0*pi)) ** (1./3.)
+#    return rad_out
 
 def random():