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
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
|
from SimPy.SimGUI import *
from SimPy.SimulationStep import *
from random import Random
__version__ = '$Revision$ $Date$ kgm'
if __name__ == '__main__':
class Source(Process):
""" Source generates customers randomly"""
def __init__(self,seed=333):
Process.__init__(self)
self.SEED = seed
def generate(self,number,interval):
rv = Random(self.SEED)
for i in range(number):
c = Customer(name = "Customer%02d"%(i,))
activate(c,c.visit(timeInBank=12.0))
t = rv.expovariate(1.0/interval)
yield hold,self,t
def NoInSystem(R):
""" The number of customers in the resource R
in waitQ and active Q"""
return (len(R.waitQ)+len(R.activeQ))
class Customer(Process):
""" Customer arrives, is served and leaves """
def __init__(self,name):
Process.__init__(self)
self.name = name
def visit(self,timeInBank=0):
arrive=now()
Qlength = [NoInSystem(counter[i]) for i in range(Nc)]
##print "%7.4f %s: Here I am. %s "%(now(),self.name,Qlength)
for i in range(Nc):
if Qlength[i] ==0 or Qlength[i]==min(Qlength): join =i ; break
yield request,self,counter[join]
wait=now()-arrive
waitMonitor.observe(wait,t=now()) #Lmond
##print "%7.4f %s: Waited %6.3f"%(now(),self.name,wait)
tib = counterRV.expovariate(1.0/timeInBank)
yield hold,self,tib
yield release,self,counter[join]
serviceMonitor.observe(now()-arrive,t=now())
if trace:
gui.writeConsole("Customer leaves at %.1d"%now())
##print "%7.4f %s: Finished "%(now(),self.name)
def showtime():
gui.topconsole.config(text="time = %s"%now())
gui.root.update()
def runStep():
if gui.noRunYet:
showwarning("SimPy warning","Run 'Start run (stepping)' first")
return
showtime()
a=simulateStep(until=gui.params.endtime)
if a[1]=="notResumable":
gui.writeConsole(text="Run ended. Status: %s"%a[0])
showtime()
def runNoStep():
showtime()
for i in range(gui.params.nrRuns):
simulate(until=gui.param.sendtime)
showtime()
gui.writeConsole("%s simulation run(s) completed\n"%(i+1))
def contStep():
return
def model():
global Nc,counter,counterRV,waitMonitor,serviceMonitor,trace,lastLeave,noRunYet,initialized
counterRV = Random(gui.params.counterseed)
sourceseed = gui.params.sourceseed
nrRuns=gui.params.nrRuns
lastLeave=0
gui.noRunYet=True
for runNr in range(nrRuns):
gui.noRunYet=False
trace=gui.params.trace
if trace:
gui.writeConsole(text='\n** Run %s'%(runNr+1))
Nc = 2
counter = [Resource(name="Clerk0"),Resource(name="Clerk1")]
gui.waitMoni=waitMonitor = Monitor(name='Waiting Times')
waitMonitor.xlab='Time'
waitMonitor.ylab='Customer waiting time'
gui.serviceMoni=serviceMonitor = Monitor(name='Service Times')
serviceMonitor.xlab='Time'
serviceMonitor.ylab='Total service time = wait+service'
initialize()
source = Source(seed = sourceseed)
activate(source,source.generate(gui.params.numberCustomers,gui.params.interval),0.0)
simulate(showtime,until=gui.params.endtime)
showtime()
lastLeave+=now()
gui.writeConsole("%s simulation run(s) completed\n"%nrRuns)
gui.writeConsole("Parameters:\n%s"%gui.params)
def modelstep():
global Nc,counter,counterRV,waitMonitor,serviceMonitor,trace,lastLeave,noRunYet
counterRV = Random(gui.params.counterseed)
sourceseed = gui.params.sourceseed
nrRuns=gui.params.nrRuns
lastLeave=0
gui.noRunYet=True
trace=gui.params.trace
if trace:
gui.writeConsole(text='\n** Run %s'%(runNr+1))
Nc = 2
counter = [Resource(name="Clerk0"),Resource(name="Clerk1")]
gui.waitMoni=waitMonitor = Monitor(name='Waiting Times')
waitMonitor.xlab='Time'
waitMonitor.ylab='Customer waiting time'
gui.serviceMoni=serviceMonitor = Monitor(name='Service Times')
serviceMonitor.xlab='Time'
serviceMonitor.ylab='Total service time = wait+service'
initialize()
source = Source(seed = sourceseed)
activate(source,source.generate(gui.params.numberCustomers,gui.params.interval),0.0)
simulateStep(until=gui.params.endtime)
gui.noRunYet=False
def statistics():
if gui.noRunYet:
showwarning(title='SimPy warning',message="Run simulation first -- no data available.")
return
aver=lastLeave/gui.params.nrRuns
gui.writeConsole(text="Average time for %s customers to get through bank: %.1f\n(%s runs)\n"\
%(gui.params.numberCustomers,aver,gui.params.nrRuns))
__doc__="""
Modified bank11.py (from Bank Tutorial) with GUI.
Model: Simulate customers arriving
at random, using a Source, requesting service
from two counters each with their own queue
random servicetime.
Uses Monitor objects to record waiting times
and total service times."""
def showAuthors():
gui.showTextBox(text="Tony Vignaux\nKlaus Muller",title="Author information")
class MyGUI(SimGUI):
def __init__(self,win,**p):
SimGUI.__init__(self,win,**p)
self.help.add_command(label="Author(s)",
command=showAuthors,underline=0)
self.view.add_command(label="Statistics",
command=statistics,underline=0)
self.run.add_command(label="Start run (event stepping)",
command=modelstep,underline=0)
self.run.add_command(label="Next event",
command=runStep,underline=0)
self.run.add_command(label="Complete run (no stepping)",
command=model,underline=0)
root=Tk()
gui=MyGUI(root,title="SimPy GUI example",doc=__doc__)
gui.params=Parameters(endtime=2000,
sourceseed=1133,
counterseed=3939393,
numberCustomers=50,
interval=10.0,
trace=0,
nrRuns=1)
gui.mainloop()
|
