1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
|
""" bank13_OO: Wait for the doorman to give a signal: *waitevent*"""
from SimPy.Simulation import Simulation, Process, Resource, SimEvent, hold, request, release, waitevent
from random import *
## Model components ------------------------
class Doorman(Process):
""" Doorman opens the door"""
def openthedoor(self):
""" He will opens the door at fixed intervals"""
for i in range(5):
yield hold, self, 30.0
self.sim.dooropen.signal()
print("%7.4f You may enter" % (self.sim.now()))
class Source(Process):
""" Source generates customers randomly"""
def generate(self, number, rate):
for i in range(number):
c = Customer(name="Customer%02d" % (i), sim=self.sim)
self.sim.activate(c, c.visit(timeInBank=12.0))
yield hold, self, expovariate(rate)
class Customer(Process):
""" Customer arrives, is served and leaves """
def visit(self, timeInBank=10):
arrive = self.sim.now()
if self.sim.dooropen.occurred:
msg = '.'
else:
msg = ' but the door is shut.'
print("%7.4f %s: Here I am%s" % (self.sim.now(), self.name, msg))
yield waitevent, self, self.sim.dooropen
print("%7.4f %s: The door is open!" % (self.sim.now(), self.name))
wait = self.sim.now() - arrive
print("%7.4f %s: Waited %6.3f" % (self.sim.now(), self.name, wait))
yield request, self, self.sim.counter
tib = expovariate(1.0 / timeInBank)
yield hold, self, tib
yield release, self, self.sim.counter
print("%7.4f %s: Finished " % (self.sim.now(), self.name))
## Model ----------------------------------
class BankModel(Simulation):
def run(self, aseed):
""" PEM """
seed(aseed)
self.dooropen = SimEvent("Door Open", sim=self)
self.counter = Resource(1, name="Clerk", sim=self)
doorman = Doorman(sim=self)
self.activate(doorman, doorman.openthedoor())
source = Source(sim=self)
self.activate(source,
source.generate(number=5, rate=0.1), at=0.0)
self.simulate(until=maxTime)
## Experiment data -------------------------
maxTime = 400.0 # minutes
seedVal = 393939
## Experiment ----------------------------------
mymodel = BankModel()
mymodel.run(aseed=seedVal)
|
