| 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 |
|
|---|