Changeset 40a87fa in sasmodels for sasmodels/models/two_power_law.py


Ignore:
Timestamp:
Aug 8, 2016 11:24:11 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:
2472141
Parents:
2d65d51
Message:

lint and latex cleanup

File:
1 edited

Legend:

Unmodified
Added
Removed
  • sasmodels/models/two_power_law.py

    r2c74c11 r40a87fa  
    1414where $q_c$ = the location of the crossover from one slope to the other, 
    1515$A$ = the scaling coefficent that sets the overall intensity of the lower Q 
    16 power law region, $m1$ = power law exponent at low Q, and $m2$ = power law  
     16power law region, $m1$ = power law exponent at low Q, and $m2$ = power law 
    1717exponent at high Q.  The scaling of the second power law region (coefficent C) 
    1818is then automatically scaled to match the first by following formula: 
     
    6262category = "shape-independent" 
    6363 
    64 #            ["name", "units", default, [lower, upper], "type", "description"], 
    65 parameters = [["coefficent_1",  "",         1.0, [-inf, inf], "", 
    66                "coefficent A in low Q region"], 
    67               ["crossover",     "1/Ang",    0.04,[0, inf],    "", 
    68                "crossover location"], 
    69               ["power_1",       "",         1.0, [0, inf],    "", 
    70                "power law exponent at low Q"], 
    71               ["power_2",       "",         4.0, [0, inf],    "", 
    72                "power law exponent at high Q"], 
    73              ] 
     64# pylint: disable=bad-whitespace, line-too-long 
     65#   ["name", "units", default, [lower, upper], "type", "description"], 
     66parameters = [ 
     67    ["coefficent_1", "",       1.0, [-inf, inf], "", "coefficent A in low Q region"], 
     68    ["crossover",    "1/Ang",  0.04,[0, inf],    "", "crossover location"], 
     69    ["power_1",      "",       1.0, [0, inf],    "", "power law exponent at low Q"], 
     70    ["power_2",      "",       4.0, [0, inf],    "", "power law exponent at high Q"], 
     71    ] 
     72# pylint: enable=bad-whitespace, line-too-long 
    7473 
    7574 
     
    8988    :return:                    Calculated intensity 
    9089    """ 
    91     iq = empty(q.shape, 'd') 
    92     idx = (q <= crossover) 
     90    result= empty(q.shape, 'd') 
     91    index = (q <= crossover) 
    9392    with errstate(divide='ignore'): 
    9493        coefficent_2 = coefficent_1 * power(crossover, power_2 - power_1) 
    95         iq[idx] = coefficent_1 * power(q[idx], -power_1) 
    96         iq[~idx] = coefficent_2 * power(q[~idx], -power_2) 
    97     return iq 
     94        result[index] = coefficent_1 * power(q[index], -power_1) 
     95        result[~index] = coefficent_2 * power(q[~index], -power_2) 
     96    return result 
    9897 
    9998Iq.vectorized = True  # Iq accepts an array of q values 
Note: See TracChangeset for help on using the changeset viewer.