Changeset 823e620 in sasmodels for sasmodels/weights.py


Ignore:
Timestamp:
Jan 29, 2016 10:46:54 PM (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:
eafc9fa
Parents:
803f835
Message:

delint

File:
1 edited

Legend:

Unmodified
Added
Removed
  • sasmodels/weights.py

    r3c56da87 r823e620  
    3636        *center* is the center of the distribution 
    3737 
    38         *lb*,*ub* are the min and max allowed values 
     38        *lb*, *ub* are the min and max allowed values 
    3939 
    4040        *relative* is True if the distribution width is proportional to the 
     
    8484    default = dict(npts=80, width=0, nsigmas=8) 
    8585    def _weights(self, center, sigma, lb, ub): 
    86         x = self._linspace(center, sigma, max(lb,1e-8), max(ub,1e-8)) 
     86        x = self._linspace(center, sigma, max(lb, 1e-8), max(ub, 1e-8)) 
    8787        px = np.exp(-0.5*(np.log(x)-center)**2)/sigma**2/(x*sigma) 
    8888        return x, px 
     
    9393    default = dict(npts=80, width=0, nsigmas=8) 
    9494    def _weights(self, center, sigma, lb, ub): 
    95         x = self._linspace(center, sigma, max(lb,1e-8), max(ub,1e-8)) 
    96         R= x/center 
     95        x = self._linspace(center, sigma, max(lb, 1e-8), max(ub, 1e-8)) 
     96        R = x/center 
    9797        z = (center/sigma)**2 
    9898        arg = z*np.log(z) + (z-1)*np.log(R) - R*z - np.log(center) - gammaln(z) 
     
    117117        # TODO: interpolate into the array dispersion using npts 
    118118        x = center + self.values*sigma 
    119         idx = (x>=lb)&(x<=ub) 
     119        idx = (x >= lb) & (x <= ub) 
    120120        x = x[idx] 
    121121        px = self.weights[idx] 
     
    124124 
    125125# dispersion name -> disperser lookup table. 
    126 models = dict((d.type,d) for d in ( 
     126models = dict((d.type, d) for d in ( 
    127127    GaussianDispersion, RectangleDispersion, 
    128128    ArrayDispersion, SchulzDispersion, LogNormalDispersion 
     
    149149    of the parameter, and false if it is an absolute width. 
    150150 
    151     Returns *(value,weight)*, where *value* and *weight* are vectors. 
     151    Returns *(value, weight)*, where *value* and *weight* are vectors. 
    152152    """ 
    153153    cls = models[disperser] 
    154154    obj = cls(n, width, nsigmas) 
    155     v,w = obj.get_weights(value, limits[0], limits[1], relative) 
    156     return v,w 
     155    v, w = obj.get_weights(value, limits[0], limits[1], relative) 
     156    return v, w 
    157157 
    158158# Hack to allow sasview dispersion objects to interoperate with sasmodels 
    159 dispersers = dict((v.__name__,k) for k,v in models.items()) 
     159dispersers = dict((v.__name__, k) for k, v in models.items()) 
    160160dispersers['DispersionModel'] = RectangleDispersion.type 
    161161 
Note: See TracChangeset for help on using the changeset viewer.