Gecko [ goons4good.org ]

Name	Size	Permission	Date
tests	[ DIR ]	drwxr-xr-x	2021-09-16 10:54
__init__.py	827 B	-rw-r--r--	2013-04-07 01:04
__init__.pyc	1022 B	-rw-r--r--	2018-04-10 19:40
__init__.pyo	1022 B	-rw-r--r--	2018-04-10 19:40
alter_code1.py	8.3 KB	-rw-r--r--	2013-04-07 01:04
alter_code1.pyc	9.49 KB	-rw-r--r--	2018-04-10 19:40
alter_code1.pyo	9.49 KB	-rw-r--r--	2018-04-10 19:40
alter_code2.py	4.53 KB	-rw-r--r--	2013-04-07 01:04
alter_code2.pyc	5.96 KB	-rw-r--r--	2018-04-10 19:40
alter_code2.pyo	5.96 KB	-rw-r--r--	2018-04-10 19:40
array_printer.py	456 B	-rw-r--r--	2013-04-07 01:04
array_printer.pyc	667 B	-rw-r--r--	2018-04-10 19:40
array_printer.pyo	667 B	-rw-r--r--	2018-04-10 19:40
arrayfns.py	2.47 KB	-rw-r--r--	2013-04-07 01:04
arrayfns.pyc	4.32 KB	-rw-r--r--	2018-04-10 19:40
arrayfns.pyo	4.32 KB	-rw-r--r--	2018-04-10 19:40
compat.py	3.11 KB	-rw-r--r--	2013-04-07 01:04
compat.pyc	4.88 KB	-rw-r--r--	2018-04-10 19:40
compat.pyo	4.88 KB	-rw-r--r--	2018-04-10 19:40
fft.py	869 B	-rw-r--r--	2013-04-07 01:04
fft.pyc	1.03 KB	-rw-r--r--	2018-04-10 19:40
fft.pyo	1.03 KB	-rw-r--r--	2018-04-10 19:40
fix_default_axis.py	7.85 KB	-rw-r--r--	2013-04-07 01:04
fix_default_axis.pyc	7.68 KB	-rw-r--r--	2018-04-10 19:40
fix_default_axis.pyo	7.68 KB	-rw-r--r--	2018-04-10 19:40
functions.py	3.55 KB	-rw-r--r--	2013-04-07 01:04
functions.pyc	7.31 KB	-rw-r--r--	2018-04-10 19:40
functions.pyo	7.31 KB	-rw-r--r--	2018-04-10 19:40
linear_algebra.py	2.13 KB	-rw-r--r--	2013-04-07 01:04
linear_algebra.pyc	3.34 KB	-rw-r--r--	2018-04-10 19:40
linear_algebra.pyo	3.34 KB	-rw-r--r--	2018-04-10 19:40
ma.py	73.96 KB	-rw-r--r--	2013-04-07 01:04
ma.pyc	83.81 KB	-rw-r--r--	2018-04-10 19:40
ma.pyo	83.81 KB	-rw-r--r--	2018-04-10 19:40
matrix.py	1.57 KB	-rw-r--r--	2013-04-07 01:04
matrix.pyc	2.31 KB	-rw-r--r--	2018-04-10 19:40
matrix.pyo	2.31 KB	-rw-r--r--	2018-04-10 19:40
misc.py	1.14 KB	-rw-r--r--	2013-04-07 01:04
misc.pyc	1.53 KB	-rw-r--r--	2018-04-10 19:40
misc.pyo	1.53 KB	-rw-r--r--	2018-04-10 19:40
mlab.py	3.36 KB	-rw-r--r--	2013-04-07 01:04
mlab.pyc	5.32 KB	-rw-r--r--	2018-04-10 19:40
mlab.pyo	5.32 KB	-rw-r--r--	2018-04-10 19:40
precision.py	4.16 KB	-rw-r--r--	2013-04-07 01:04
precision.pyc	3.85 KB	-rw-r--r--	2018-04-10 19:40
precision.pyo	3.85 KB	-rw-r--r--	2018-04-10 19:40
random_array.py	11.18 KB	-rw-r--r--	2013-04-07 01:04
random_array.pyc	13.89 KB	-rw-r--r--	2018-04-10 19:40
random_array.pyo	13.89 KB	-rw-r--r--	2018-04-10 19:40
rng.py	3.67 KB	-rw-r--r--	2013-04-07 01:04
rng.pyc	6.74 KB	-rw-r--r--	2018-04-10 19:40
rng.pyo	6.74 KB	-rw-r--r--	2018-04-10 19:40
rng_stats.py	1.31 KB	-rw-r--r--	2013-04-07 01:04
rng_stats.pyc	1.83 KB	-rw-r--r--	2018-04-10 19:40
rng_stats.pyo	1.83 KB	-rw-r--r--	2018-04-10 19:40
setup.py	335 B	-rw-r--r--	2013-04-07 01:04
setup.pyc	664 B	-rw-r--r--	2018-04-10 19:40
setup.pyo	664 B	-rw-r--r--	2018-04-10 19:40
setupscons.py	282 B	-rw-r--r--	2013-04-07 01:04
setupscons.pyc	613 B	-rw-r--r--	2018-04-10 19:40
setupscons.pyo	613 B	-rw-r--r--	2018-04-10 19:40
typeconv.py	1.58 KB	-rw-r--r--	2013-04-07 01:04
typeconv.pyc	1.59 KB	-rw-r--r--	2018-04-10 19:40
typeconv.pyo	1.59 KB	-rw-r--r--	2018-04-10 19:40
ufuncs.py	1.2 KB	-rw-r--r--	2013-04-07 01:04
ufuncs.pyc	1.79 KB	-rw-r--r--	2018-04-10 19:40
ufuncs.pyo	1.79 KB	-rw-r--r--	2018-04-10 19:40
user_array.py	182 B	-rw-r--r--	2013-04-07 01:04
user_array.pyc	370 B	-rw-r--r--	2018-04-10 19:40
user_array.pyo	370 B	-rw-r--r--	2018-04-10 19:40

Rename

# This module re-creates the RNG interface from Numeric
# Replace import RNG with import numpy.oldnumeric.rng as RNG
#
# It is for backwards compatibility only.

__all__ = ['CreateGenerator','ExponentialDistribution','LogNormalDistribution',
           'NormalDistribution', 'UniformDistribution', 'error', 'ranf',
           'default_distribution', 'random_sample', 'standard_generator']

import numpy.random.mtrand as mt
import math

class error(Exception):
    pass

class Distribution(object):
    def __init__(self, meth, *args):
        self._meth = meth
        self._args = args

def density(self,x):
        raise NotImplementedError

def __call__(self, x):
        return self.density(x)

def _onesample(self, rng):
        return getattr(rng, self._meth)(*self._args)

def _sample(self, rng, n):
        kwds = {'size' : n}
        return getattr(rng, self._meth)(*self._args, **kwds)

class ExponentialDistribution(Distribution):
    def __init__(self, lambda_):
        if (lambda_ <= 0):
            raise error, "parameter must be positive"
        Distribution.__init__(self, 'exponential', lambda_)

def density(x):
        if x < 0:
            return 0.0
        else:
            lambda_ = self._args[0]
            return lambda_*math.exp(-lambda_*x)

class LogNormalDistribution(Distribution):
    def __init__(self, m, s):
        m = float(m)
        s = float(s)
        if (s <= 0):
            raise error, "standard deviation must be positive"
        Distribution.__init__(self, 'lognormal', m, s)
        sn = math.log(1.0+s*s/(m*m));
        self._mn = math.log(m)-0.5*sn
        self._sn = math.sqrt(sn)
        self._fac = 1.0/math.sqrt(2*math.pi)/self._sn

def density(x):
        m,s = self._args
        y = (math.log(x)-self._mn)/self._sn
        return self._fac*math.exp(-0.5*y*y)/x

class NormalDistribution(Distribution):
    def __init__(self, m, s):
        m = float(m)
        s = float(s)
        if (s <= 0):
            raise error, "standard deviation must be positive"
        Distribution.__init__(self, 'normal', m, s)
        self._fac = 1.0/math.sqrt(2*math.pi)/s

def density(x):
        m,s = self._args
        y = (x-m)/s
        return self._fac*math.exp(-0.5*y*y)

class UniformDistribution(Distribution):
    def __init__(self, a, b):
        a = float(a)
        b = float(b)
        width = b-a
        if (width <=0):
            raise error, "width of uniform distribution must be > 0"
        Distribution.__init__(self, 'uniform', a, b)
        self._fac = 1.0/width

def density(x):
        a, b = self._args
        if (x < a) or (x >= b):
            return 0.0
        else:
            return self._fac

default_distribution = UniformDistribution(0.0,1.0)

class CreateGenerator(object):
    def __init__(self, seed, dist=None):
        if seed <= 0:
            self._rng = mt.RandomState()
        elif seed > 0:
            self._rng = mt.RandomState(seed)
        if dist is None:
            dist = default_distribution
        if not isinstance(dist, Distribution):
            raise error, "Not a distribution object"
        self._dist = dist

def ranf(self):
        return self._dist._onesample(self._rng)

def sample(self, n):
        return self._dist._sample(self._rng, n)

standard_generator = CreateGenerator(-1)

def ranf():
    "ranf() = a random number from the standard generator."
    return standard_generator.ranf()

def random_sample(*n):
    """random_sample(n) = array of n random numbers;

random_sample(n1, n2, ...)= random array of shape (n1, n2, ..)"""

if not n:
        return standard_generator.ranf()
    m = 1
    for i in n:
        m = m * i
    return standard_generator.sample(m).reshape(*n)