Changeset c6652bb in sasmodels

Timestamp:

Mar 21, 2016 5:30:00 PM (9 years ago)

Author:

butler

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:

693570c

Parents:

f7930be

Message:

Review completeness of two_power_law and corrected equation in docs

File:

• sasmodels/models/two_power_law.py (modified) (5 diffs)

sasmodels/models/two_power_law.py

@@ -1 +1 @@
 r"""
-This model calculates an empirical functional form for SAS data characterized
-by two power laws.
-
 Definition
 ----------
@@ -11 +8 @@
 
     I(q) = \begin{cases}
-    A q^{-m1} + \text{background} & q <= qc \\
-    C q^{-m2} + \text{background} & q > qc
+    A q^{-m1} + \text{background} & q <= q_c \\
+    C q^{-m2} + \text{background} & q > q_c
     \end{cases}
 
-where $qc$ = the location of the crossover from one slope to the other,
-$A$ = the scaling coefficent that sets the overall intensity of the lower Q power law region,
-$m1$ = power law exponent at low Q,
-$m2$ = power law exponent at high Q
-
-The scaling of the second power law region (coefficent C) is then automatically scaled to
-match the first by following formula
+where $q_c$ = the location of the crossover from one slope to the other,
+$A$ = the scaling coefficent that sets the overall intensity of the lower Q
+power law region, $m1$ = power law exponent at low Q, and $m2$ = power law 
+exponent at high Q.  The scaling of the second power law region (coefficent C)
+is then automatically scaled to match the first by following formula:
 
 .. math::
-    C = \frac{A}{qc^{-m1} qc^{-m2}}
+    C = \frac{A q_c^{m2}}{q_c^{m1}}
 
 .. note::
@@ -42 +37 @@
 None.
 
+**Author:** NIST IGOR/DANSE **on:** pre 2010
+
+**Last Modified by:** Wojciech Wpotrzebowski **on:** February 18, 2016
+
+**Last Reviewed by:** Paul Butler **on:** March 21, 2016
 """
 
@@ -50 +50 @@
 
 name = "two_power_law"
-title = "Two Power Law intensity calculation"
+title = "This model calculates an empirical functional form for SAS data \
+characterized by two power laws."
 description = """
-            I(q) = coef_A*pow(qval,-1.0*power1) + background for q<=qc
-            =C*pow(qval,-1.0*power2) + background for q>qc
-            where C=coef_A*pow(qc,-1.0*power1)/pow(qc,-1.0*power2).
+            I(q) = coef_A*pow(qval,-1.0*power1) + background for q<=q_c
+            =C*pow(qval,-1.0*power2) + background for q>q_c
+            where C=coef_A*pow(q_c,-1.0*power1)/pow(q_c,-1.0*power2).
 
             coef_A = scaling coefficent
-            qc = crossover location [1/A]
+            q_c = crossover location [1/A]
             power_1 (=m1) = power law exponent at low Q
             power_2 (=m2) = power law exponent at high Q
@@ -64 +65 @@
 category = "shape-independent"
 
-# pylint: disable=bad-whitespace, line-too-long
 #            ["name", "units", default, [lower, upper], "type", "description"],
-parameters = [["coefficent_1",  "",         1.0, [-inf, inf], "", "coefficent A in low Q region"],
-              ["crossover",     "1/Ang",    0.04,[0, inf],    "", "crossover location"],
-              ["power_1",       "",         1.0, [0, inf],    "", "power law exponent at low Q"],
-              ["power_2",       "",         4.0, [0, inf],    "", "power law exponent at high Q"],
+parameters = [["coefficent_1",  "",         1.0, [-inf, inf], "",
+               "coefficent A in low Q region"],
+              ["crossover",     "1/Ang",    0.04,[0, inf],    "",
+               "crossover location"],
+              ["power_1",       "",         1.0, [0, inf],    "",
+               "power law exponent at low Q"],
+              ["power_2",       "",         4.0, [0, inf],    "",
+               "power law exponent at high Q"],
              ]
-# pylint: enable=bad-whitespace, line-too-long