random.expovariate

All Samples(694) | Call(533) | Derive(0) | Import(161)

Exponential distribution.

lambd is 1.0 divided by the desired mean.  It should be
nonzero.  (The parameter would be called "lambda", but that is
a reserved word in Python.)  Returned values range from 0 to
positive infinity if lambd is positive, and from negative
infinity to 0 if lambd is negative.

src/x/g/xgcode-HEAD/PathSampler.py xgcode(Download)

        events.append((t, state))
    while True:
        rate = -rate_matrix[(state, state)]
        t += random.expovariate(rate)
        if t > path_length:
            if state == terminal_state:
                return events

    state = initial_state
    while True:
        rate_away = -rate_matrix[(state, state)]
        t += random.expovariate(rate_away)
        if t > path_length:
            if state == terminal_state:
                return events

src/x/g/xgcode-HEAD/BranchLengthSampler.py xgcode(Download)

    def __call__(self):
        """
        @return: a branch length drawn from an exponential distribution
        """
        mean = 0.1
        return random.expovariate(1/mean)

src/l/e/lepton-1.0b2/examples/fireworks.py lepton(Download)

 
import os
import math
from random import expovariate, uniform, gauss
from pyglet import image
from pyglet.gl import *

def fire(dt=None):
	Kaboom()
	pyglet.clock.schedule_once(fire, expovariate(1.0 / (MEAN_FIRE_INTERVAL - 1)) + 1)
 
fire()
win.set_visible(True)
pyglet.clock.schedule_interval(default_system.update, (1.0/30.0))

src/l/e/lepton-1.0b2/examples/magnet.py lepton(Download)

from lepton.emitter import StaticEmitter, PerParticleEmitter
from lepton.controller import Movement, Magnet, Collector, Lifetime, Fader
from lepton.domain import Sphere, Point, Disc
from random import expovariate
 
win = pyglet.window.Window(resizable=True, visible=False)
win.clear()

def summon(dt=None):
	if len(electrons) < max_electrons:
		electron_emitter.emit(1 ,electrons)
	pyglet.clock.schedule_once(summon, expovariate(1.0)+1.0)
 
summon()

src/s/p/speechresearch-HEAD/trunk/gmisclib/mcmc_cooperate.py speechresearch(Download)

def test0():
	import time
	test = connection(*test_args)
	test.set("a", "A")
	# test.barrier(Barrier(1))
	time.sleep(random.expovariate(1.0))
	test.set("a", "B")
	time.sleep(random.expovariate(1.0))
	test.barrier(Barrier(2))
	time.sleep(random.expovariate(1.0))
	for x in test.get_list("a"):
		assert x == "B"
	time.sleep(random.expovariate(1.0))

def test0s():
	import time
	test = connection(*test_args)
	test.set("a", "A")
	time.sleep(random.expovariate(1.0))
	test.barrier(Barrier(1))
	print "Barrier 1s"
	time.sleep(random.expovariate(1.0))
	test.barrier(Barrier(1, 3))
	print "Barrier1.3s"
	time.sleep(random.expovariate(1.0))

	print "test.spread 0s tmp=", tmp
	for x in test.get_list("a"):
		assert x == "x"
	time.sleep(random.expovariate(1.0))
	test.barrier(Barrier(3))
	test.set("a", "B")
	test.close()

def test1():
	import time
	test = connection(*test_args)
	print "rank=", test.rank()
	test.set("key", "value")
	assert "value" in test.get_list("key")
	time.sleep(random.expovariate(100.0))
	test.set("k2", "2", "k", "wahoonie")
	test.barrier(Barrier(1))
	print 'Barrier 1'
	time.sleep(random.expovariate(100.0))
	test.barrier(Barrier(1, 1))
	time.sleep(random.expovariate(100.0))

	test.set("k2", "1", "k", "no_value")
	test.barrier(Barrier(2))
	print 'Barrier 2'
	time.sleep(random.expovariate(100.0))
	test.barrier(Barrier(2, 0, 1))
	tmp = test.get_list("k2")
	assert tmp[0] == "1", "tmp=%s" % tmp
	test.barrier(Barrier(2, 1))
	print 'Barrier 2.1'
	test.set("k2", "2", "k", "wahoonie")
	time.sleep(random.expovariate(100.0))

	assert tmp[0] == "wahoonie", "tmp=%s" % tmp
	test.set("k3", "x y")
	assert "x y" in test.get_list("k3")
	time.sleep(random.expovariate(100.0))
	assert test.get_combined("k2", 'float_median') == 2.0
	test.barrier(Barrier(3))
	time.sleep(random.expovariate(100.0))
	test.barrier(Barrier(3))
	time.sleep(random.expovariate(100.0))
	assert test.barrier(Barrier(2), exc=False) == Barrier(3)
	time.sleep(random.expovariate(100.0))
	assert test.barrier(Barrier(1), exc=False) == Barrier(3)
	time.sleep(random.expovariate(100.0))

	assert test.barrier(Barrier(1), exc=False) == Barrier(3)
	time.sleep(random.expovariate(100.0))
	assert test.barrier(Barrier(3)) is None
	time.sleep(random.expovariate(100.0))
	assert test.barrier(Barrier(4)) is None
	assert test.spread("k3", "x", Barrier(6))[0] == "x"
	time.sleep(random.expovariate(100.0))
	assert test.get_consensus("k3", "a", Barrier(6, 3), "string_median") == "a"
	test.barrier(Barrier(7))
	print "barrier 7"
	tmp = test.spread("k3", "x", Barrier(7, 3))
	assert tmp[0] == "x", "tmp= %s" % tmp
	time.sleep(random.expovariate(100.0))

src/s/p/speechresearch-HEAD/gmisclib/mcmc_cooperate.py speechresearch(Download)

def test0():
	import time
	test = connection(*test_args)
	test.set("a", "A")
	# test.barrier(Barrier(1))
	time.sleep(random.expovariate(1.0))
	test.set("a", "B")
	time.sleep(random.expovariate(1.0))
	test.barrier(Barrier(2))
	time.sleep(random.expovariate(1.0))
	for x in test.get_list("a"):
		assert x == "B"
	time.sleep(random.expovariate(1.0))

def test0s():
	import time
	test = connection(*test_args)
	test.set("a", "A")
	time.sleep(random.expovariate(1.0))
	test.barrier(Barrier(1))
	print "Barrier 1s"
	time.sleep(random.expovariate(1.0))
	test.barrier(Barrier(1, 3))
	print "Barrier1.3s"
	time.sleep(random.expovariate(1.0))

	print "test.spread 0s tmp=", tmp
	for x in test.get_list("a"):
		assert x == "x"
	time.sleep(random.expovariate(1.0))
	test.barrier(Barrier(3))
	test.set("a", "B")
	test.close()

def test1():
	import time
	test = connection(*test_args)
	print "rank=", test.rank()
	test.set("key", "value")
	assert "value" in test.get_list("key")
	time.sleep(random.expovariate(100.0))
	test.set("k2", "2", "k", "wahoonie")
	test.barrier(Barrier(1))
	print 'Barrier 1'
	time.sleep(random.expovariate(100.0))
	test.barrier(Barrier(1, 1))
	time.sleep(random.expovariate(100.0))

	test.set("k2", "1", "k", "no_value")
	test.barrier(Barrier(2))
	print 'Barrier 2'
	time.sleep(random.expovariate(100.0))
	test.barrier(Barrier(2, 0, 1))
	tmp = test.get_list("k2")
	assert tmp[0] == "1", "tmp=%s" % tmp
	test.barrier(Barrier(2, 1))
	print 'Barrier 2.1'
	test.set("k2", "2", "k", "wahoonie")
	time.sleep(random.expovariate(100.0))

	assert tmp[0] == "wahoonie", "tmp=%s" % tmp
	test.set("k3", "x y")
	assert "x y" in test.get_list("k3")
	time.sleep(random.expovariate(100.0))
	assert test.get_combined("k2", 'float_median') == 2.0
	test.barrier(Barrier(3))
	time.sleep(random.expovariate(100.0))
	test.barrier(Barrier(3))
	time.sleep(random.expovariate(100.0))
	assert test.barrier(Barrier(2), exc=False) == Barrier(3)
	time.sleep(random.expovariate(100.0))
	assert test.barrier(Barrier(1), exc=False) == Barrier(3)
	time.sleep(random.expovariate(100.0))

	assert test.barrier(Barrier(1), exc=False) == Barrier(3)
	time.sleep(random.expovariate(100.0))
	assert test.barrier(Barrier(3)) is None
	time.sleep(random.expovariate(100.0))
	assert test.barrier(Barrier(4)) is None
	assert test.spread("k3", "x", Barrier(6))[0] == "x"
	time.sleep(random.expovariate(100.0))
	assert test.get_consensus("k3", "a", Barrier(6, 3), "string_median") == "a"
	test.barrier(Barrier(7))
	print "barrier 7"
	tmp = test.spread("k3", "x", Barrier(7, 3))
	assert tmp[0] == "x", "tmp= %s" % tmp
	time.sleep(random.expovariate(100.0))

src/s/i/simpy-HEAD/SimPy/SimPy/trunk/SimPyModels/bank08renege_OO.py simpy(Download)

 from TheBank tutorial. (KGM)
"""
from SimPy.Simulation  import *
from random import expovariate, seed, uniform
 
## Model components ------------------------

    def generate(self,number,interval,counter):       
        for i in range(number):
            c = Customer(name = "Customer%02d"%(i,),sim=self.sim)
            self.sim.activate(c,c.visit(counter,timeInBank=12.0))
            t = expovariate(1.0/interval)
            yield hold,self,t

        if self.acquired(counter):
            wait = self.sim.now()-arrive
            print "%7.4f %s: Waited %6.3f"%(self.sim.now(),self.name,wait)
            tib = expovariate(1.0/timeInBank)            
            yield hold,self,tib                          
            yield release,self,counter
            print "%7.4f %s: Finished"%(self.sim.now(),self.name)

src/s/i/simpy-HEAD/SimPy/SimPy/trunk/SimPyModels/bank08renege.py simpy(Download)

 from TheBank tutorial. (KGM)
"""
from SimPy.Simulation  import *
from random import expovariate, seed, uniform
 
## Model components ------------------------

    def generate(self,number,interval,counter):       
        for i in range(number):
            c = Customer(name = "Customer%02d"%(i,))
            activate(c,c.visit(counter,timeInBank=12.0))
            t = expovariate(1.0/interval)
            yield hold,self,t

        if self.acquired(counter):
            wait=now()-arrive
            print "%7.4f %s: Waited %6.3f"%(now(),self.name,wait)
            tib = expovariate(1.0/timeInBank)            
            yield hold,self,tib                          
            yield release,self,counter
            print "%7.4f %s: Finished"%(now(),self.name)

src/s/i/simpy-HEAD/SimPy/Python3/SimPyModels/bank08renege_OO.py simpy(Download)

 from TheBank tutorial. (KGM)
"""
from SimPy.Simulation  import *
from random import expovariate, seed, uniform
 
## Model components ------------------------

    def generate(self,number,interval,counter):       
        for i in range(number):
            c = Customer(name = "Customer%02d"%(i,),sim=self.sim)
            self.sim.activate(c,c.visit(counter,timeInBank=12.0))
            t = expovariate(1.0/interval)
            yield hold,self,t

        if self.acquired(counter):
            wait = self.sim.now()-arrive
            print("%7.4f %s: Waited %6.3f"%(self.sim.now(),self.name,wait))
            tib = expovariate(1.0/timeInBank)            
            yield hold,self,tib                          
            yield release,self,counter
            print("%7.4f %s: Finished"%(self.sim.now(),self.name))

src/s/i/simpy-HEAD/SimPy/Python3/SimPyModels/bank08renege.py simpy(Download)

 from TheBank tutorial. (KGM)
"""
from SimPy.Simulation  import *
from random import expovariate, seed, uniform
 
## Model components ------------------------

    def generate(self,number,interval,counter):       
        for i in range(number):
            c = Customer(name = "Customer%02d"%(i,))
            activate(c,c.visit(counter,timeInBank=12.0))
            t = expovariate(1.0/interval)
            yield hold,self,t

        if self.acquired(counter):
            wait=now()-arrive
            print("%7.4f %s: Waited %6.3f"%(now(),self.name,wait))
            tib = expovariate(1.0/timeInBank)            
            yield hold,self,tib                          
            yield release,self,counter
            print("%7.4f %s: Finished"%(now(),self.name))

1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 Next