def gsl_uniform_sequence(n,seed=None):
"""
Return sample sequence of length n from a uniform distribution.
"""
try:
import pygsl.rng
except ImportError:
print "Import error: not able to import pygsl"
return
rng=pygsl.rng.rng()
random._inst = random.Random()
if seed is None:
seed=random.randint(1,2**32-1)
rng.set(seed)
return rng.uniform(n)
def gsl_uniform_sequence(n, seed=None):
"""
Return sample sequence of length n from a uniform distribution.
"""
try:
import pygsl.rng
except ImportError:
print "Import error: not able to import pygsl"
return
rng = pygsl.rng.rng()
random._inst = random.Random()
if seed is None:
seed = random.randint(1, 2**32 - 1)
rng.set(seed)
return rng.uniform(n)
def gsl_powerlaw_sequence(n,exponent=2.0,scale=1.0,seed=None):
"""
Return sample sequence of length n from a power law distribution.
"""
try:
import pygsl.rng
except ImportError:
print "Import error: not able to import pygsl"
return
rng=pygsl.rng.rng()
random._inst = random.Random()
if seed is None:
seed=random.randint(1,2**32-1)
rng.set(seed)
return rng.pareto(exponent-1,scale,n)
def gsl_powerlaw_sequence(n, exponent=2.0, scale=1.0, seed=None):
"""
Return sample sequence of length n from a power law distribution.
"""
try:
import pygsl.rng
except ImportError:
print "Import error: not able to import pygsl"
return
rng = pygsl.rng.rng()
random._inst = random.Random()
if seed is None:
seed = random.randint(1, 2**32 - 1)
rng.set(seed)
return rng.pareto(exponent - 1, scale, n)
