1 | #include <stdio.h>
|
---|
2 | #include <tchar.h>
|
---|
3 |
|
---|
4 | #include <string>
|
---|
5 |
|
---|
6 | #include <boost/bind.hpp>
|
---|
7 | #include <boost/coroutine/all.hpp>
|
---|
8 | #include <boost/coroutine/attributes.hpp>
|
---|
9 | #include <boost/coroutine/detail/flags.hpp>
|
---|
10 | #include <boost/timer/timer.hpp>
|
---|
11 | #include <boost/lexical_cast.hpp>
|
---|
12 | #include <boost/chrono.hpp>
|
---|
13 | #include <boost/math/distributions.hpp>
|
---|
14 | #include <boost/random/variate_generator.hpp>
|
---|
15 | #include <boost/random/mersenne_twister.hpp>
|
---|
16 | #include <boost/random/uniform_int_distribution.hpp>
|
---|
17 | #include <boost/random/uniform_real_distribution.hpp>
|
---|
18 | #include <boost/random/triangle_distribution.hpp>
|
---|
19 | #include <boost/random/exponential_distribution.hpp>
|
---|
20 | #include <boost/random/gamma_distribution.hpp>
|
---|
21 | #include <boost/random/normal_distribution.hpp>
|
---|
22 | #include <boost/random/lognormal_distribution.hpp>
|
---|
23 |
|
---|
24 | using namespace std;
|
---|
25 |
|
---|
26 | ///////////////////////////////////////////////////////////////////////////////
|
---|
27 | class job
|
---|
28 | {
|
---|
29 |
|
---|
30 | public:
|
---|
31 | job() {};
|
---|
32 | job(string _name, void(*processFunction)(boost::coroutines::symmetric_coroutine<job&>::yield_type&)) : name(_name), coroProc(processFunction), ID(nextID++)
|
---|
33 | {
|
---|
34 | nextEventTime = 0;
|
---|
35 | priority = 1;
|
---|
36 | }
|
---|
37 | ~job() {};
|
---|
38 |
|
---|
39 | boost::coroutines::symmetric_coroutine<job&>::yield_type *yield;
|
---|
40 | boost::coroutines::symmetric_coroutine<job&>::call_type coroProc;
|
---|
41 |
|
---|
42 | int ID;
|
---|
43 | static int nextID;
|
---|
44 | string name;
|
---|
45 | string STATUS;
|
---|
46 | int priority;
|
---|
47 | void setNextEventTime(double net) { nextEventTime = net; }
|
---|
48 | double getNextEventTime() { return nextEventTime; }
|
---|
49 |
|
---|
50 | private:
|
---|
51 | double nextEventTime;
|
---|
52 |
|
---|
53 | };
|
---|
54 | ///////////////////////////////////////////////////////////////////////////////
|
---|
55 |
|
---|
56 | ///////////////////////////////////////////////////////////////////////////////
|
---|
57 | bool compByNextEventTime(job* a, job* b)
|
---|
58 | {
|
---|
59 | if (a->getNextEventTime() == b->getNextEventTime())
|
---|
60 | {
|
---|
61 | if (a->priority == b->priority)
|
---|
62 | {
|
---|
63 | return false;
|
---|
64 | }
|
---|
65 | else
|
---|
66 | {
|
---|
67 | return a->priority > b->priority;
|
---|
68 | }
|
---|
69 | }
|
---|
70 | else
|
---|
71 | {
|
---|
72 | return a->getNextEventTime() < b->getNextEventTime();
|
---|
73 | }
|
---|
74 | }
|
---|
75 | ///////////////////////////////////////////////////////////////////////////////
|
---|
76 |
|
---|
77 | ///////////////////////////////////////////////////////////////////////////////
|
---|
78 | class jobShop
|
---|
79 | {
|
---|
80 | public:
|
---|
81 |
|
---|
82 | jobShop()
|
---|
83 | {
|
---|
84 | modelTime = 0;;
|
---|
85 | }
|
---|
86 | ~jobShop() {};
|
---|
87 | jobShop(const jobShop& sM)
|
---|
88 | {
|
---|
89 | modelTime = sM.modelTime;
|
---|
90 | }
|
---|
91 |
|
---|
92 | vector<job*> jobVector;
|
---|
93 | void insertProcess(job* j)
|
---|
94 | {
|
---|
95 | auto it = std::upper_bound(jobVector.begin(), jobVector.end(), j, compByNextEventTime);
|
---|
96 | jobVector.insert(it, j);
|
---|
97 | }
|
---|
98 |
|
---|
99 | double simtime()
|
---|
100 | {
|
---|
101 | return modelTime;
|
---|
102 | }
|
---|
103 | void setModelTime(double mt) { modelTime = mt; }
|
---|
104 |
|
---|
105 |
|
---|
106 | private:
|
---|
107 | double modelTime;
|
---|
108 | };
|
---|
109 | ///////////////////////////////////////////////////////////////////////////////
|
---|
110 |
|
---|
111 | ///////////////////////////////////////////////////////////////////////////////
|
---|
112 | typedef boost::coroutines::symmetric_coroutine<void> coro_t;
|
---|
113 | typedef boost::coroutines::symmetric_coroutine<job&> coro_simProcess_t;
|
---|
114 |
|
---|
115 | int job::nextID = 0;
|
---|
116 |
|
---|
117 | jobShop *m;
|
---|
118 | boost::random::mt19937 gen;
|
---|
119 |
|
---|
120 | long uniform_int(long min, long max);
|
---|
121 | void hold(double time);
|
---|
122 | void createJob(string _name, void(*processFunction)(coro_simProcess_t::yield_type &yield));
|
---|
123 | void activateNextJob();
|
---|
124 | void jobCoroFn(coro_simProcess_t::yield_type &yield);
|
---|
125 | void jobCoroFn2(coro_simProcess_t::yield_type &yield);
|
---|
126 |
|
---|
127 | ///////////////////////////////////////////////////////////////////////////////
|
---|
128 | long uniform_int(long min, long max)
|
---|
129 | {
|
---|
130 | boost::random::uniform_int_distribution<> dist(min, max);
|
---|
131 | return dist(gen);
|
---|
132 | }
|
---|
133 | ///////////////////////////////////////////////////////////////////////////////
|
---|
134 |
|
---|
135 | ///////////////////////////////////////////////////////////////////////////////
|
---|
136 | void hold(double time)
|
---|
137 | {
|
---|
138 |
|
---|
139 | job* j = m->jobVector.front();
|
---|
140 | m->jobVector.erase(m->jobVector.begin());
|
---|
141 |
|
---|
142 | coro_simProcess_t::yield_type &yield = *j->yield;
|
---|
143 | j->setNextEventTime(time + m->simtime());
|
---|
144 | m->insertProcess(j);
|
---|
145 |
|
---|
146 | yield();
|
---|
147 | }
|
---|
148 | ///////////////////////////////////////////////////////////////////////////////
|
---|
149 |
|
---|
150 | ///////////////////////////////////////////////////////////////////////////////
|
---|
151 | void createJob(string _name, void(*processFunction)(coro_simProcess_t::yield_type &yield))
|
---|
152 | {
|
---|
153 | job* j = new job(_name, processFunction);
|
---|
154 | j->setNextEventTime(m->simtime());
|
---|
155 | m->insertProcess(j);
|
---|
156 | }
|
---|
157 | ///////////////////////////////////////////////////////////////////////////////
|
---|
158 |
|
---|
159 | ///////////////////////////////////////////////////////////////////////////////
|
---|
160 | void activateNextJob()
|
---|
161 | {
|
---|
162 |
|
---|
163 | while (!m->jobVector.empty())
|
---|
164 | {
|
---|
165 |
|
---|
166 | job* j = m->jobVector.front();
|
---|
167 |
|
---|
168 | m->setModelTime(j->getNextEventTime());
|
---|
169 | j->coroProc(*j);
|
---|
170 |
|
---|
171 | if (j->STATUS == "COMPLETE")
|
---|
172 | {
|
---|
173 | m->jobVector.erase(m->jobVector.begin());
|
---|
174 | cout << endl;
|
---|
175 | cout << "Finished job " << j->name << " " << m->simtime() << endl;
|
---|
176 | cout << endl;
|
---|
177 | delete j;
|
---|
178 | }
|
---|
179 |
|
---|
180 | }
|
---|
181 | }
|
---|
182 | ///////////////////////////////////////////////////////////////////////////////
|
---|
183 |
|
---|
184 | ///////////////////////////////////////////////////////////////////////////////
|
---|
185 | void jobCoroFn(coro_simProcess_t::yield_type &yield)
|
---|
186 | {
|
---|
187 |
|
---|
188 | job& j = yield.get();
|
---|
189 | j.yield = &yield;
|
---|
190 |
|
---|
191 | hold(0);
|
---|
192 |
|
---|
193 | cout << j.name << " 1 " << m->simtime() << endl;
|
---|
194 | if (j.name == "P1")
|
---|
195 | hold(uniform_int(0, 10));
|
---|
196 | else
|
---|
197 | hold(10.0);
|
---|
198 |
|
---|
199 | cout << j.name << " 2 " << m->simtime() << endl;
|
---|
200 | if (j.name == "P1")
|
---|
201 | hold(uniform_int(5, 10));
|
---|
202 | else
|
---|
203 | hold(10.0);
|
---|
204 |
|
---|
205 | hold(0);
|
---|
206 | cout << j.name << " 3 " << m->simtime() << endl;
|
---|
207 | if (j.name == "P1")
|
---|
208 | hold(uniform_int(0, 10));
|
---|
209 | else
|
---|
210 | hold(10.0);
|
---|
211 |
|
---|
212 | cout << j.name << " 4 " << m->simtime() << endl;
|
---|
213 | if (j.name == "P1")
|
---|
214 | hold(uniform_int(5, 10));
|
---|
215 | else
|
---|
216 | hold(10.00);
|
---|
217 |
|
---|
218 | cout << j.name << " 5 " << m->simtime() << endl;
|
---|
219 | if (j.name == "P1")
|
---|
220 | hold(uniform_int(5, 10));
|
---|
221 | else
|
---|
222 | {
|
---|
223 | hold(10.0);
|
---|
224 | }
|
---|
225 |
|
---|
226 | cout << j.name << " 6 " << m->simtime() << endl;
|
---|
227 |
|
---|
228 | j.STATUS = "COMPLETE";
|
---|
229 | }
|
---|
230 | ///////////////////////////////////////////////////////////////////////////////
|
---|
231 |
|
---|
232 | ///////////////////////////////////////////////////////////////////////////////
|
---|
233 | void jobCoroFn2(coro_simProcess_t::yield_type &yield)
|
---|
234 | {
|
---|
235 | job& j = yield.get();
|
---|
236 | j.yield = &yield;
|
---|
237 |
|
---|
238 | cout << j.name << " A " << m->simtime() << endl;
|
---|
239 | hold(uniform_int(1, 10));
|
---|
240 |
|
---|
241 | cout << j.name << " B " << m->simtime() << endl;
|
---|
242 | createJob("P3", jobCoroFn);
|
---|
243 | hold(uniform_int(1, 10));
|
---|
244 |
|
---|
245 | cout << j.name << " C " << m->simtime() << endl;
|
---|
246 | hold(uniform_int(5, 10));
|
---|
247 |
|
---|
248 | cout << j.name << " D " << m->simtime() << endl;
|
---|
249 | j.STATUS = "COMPLETE";
|
---|
250 |
|
---|
251 | j.yield;
|
---|
252 | }
|
---|
253 | ///////////////////////////////////////////////////////////////////////////////
|
---|
254 |
|
---|
255 | int main()
|
---|
256 | {
|
---|
257 | boost::timer::cpu_timer timer;
|
---|
258 |
|
---|
259 | m = new jobShop();
|
---|
260 | gen.seed(49);
|
---|
261 |
|
---|
262 | createJob("P1", jobCoroFn);
|
---|
263 | createJob("P2", jobCoroFn2);
|
---|
264 |
|
---|
265 | activateNextJob();
|
---|
266 |
|
---|
267 | cout << endl;
|
---|
268 | cout << "Simulation Finished" << endl;
|
---|
269 | std::cout << timer.format() << endl;
|
---|
270 |
|
---|
271 | delete m;
|
---|
272 |
|
---|
273 | getchar();
|
---|
274 | return 0;
|
---|
275 | }
|
---|
276 |
|
---|