source/libr2c2/trainrouteplanner.cpp

   1 #include <msp/core/maputils.h>
   2 #include "catalogue.h"
   3 #include "layout.h"
   4 #include "route.h"
   5 #include "train.h"
   6 #include "trainroutemetric.h"
   7 #include "trainrouteplanner.h"
   8 #include "trainrouter.h"
   9 #include "vehicle.h"
  10
  11 using namespace std;
  12 using namespace Msp;
  13
  14 namespace R2C2 {
  15
  16 TrainRoutePlanner::TrainRoutePlanner(Layout &layout):
  17         goal(0),
  18         result(PENDING),
  19         thread(0)
  20 {
  21         const map<unsigned, Train *> &trains = layout.get_trains();
  22         for(map<unsigned, Train *>::const_iterator i=trains.begin(); i!=trains.end(); ++i)
  23         {
  24                 TrainRoutingInfo info(*i->second);
  25                 if(info.router && info.router->get_destination())
  26                         routed_trains.push_back(info);
  27         }
  28 }
  29
  30 TrainRoutePlanner::~TrainRoutePlanner()
  31 {
  32         delete thread;
  33 }
  34
  35 TrainRoutePlanner::Result TrainRoutePlanner::plan()
  36 {
  37         prepare_plan();
  38         create_plan();
  39         if(result==PENDING)
  40                 finalize_plan();
  41
  42         return result;
  43 }
  44
  45 void TrainRoutePlanner::plan_async()
  46 {
  47         if(thread)
  48                 throw logic_error("already planning");
  49
  50         prepare_plan();
  51         thread = new PlanningThread(*this);
  52 }
  53
  54 TrainRoutePlanner::Result TrainRoutePlanner::check()
  55 {
  56         if(result==PENDING && goal)
  57         {
  58                 finalize_plan();
  59                 delete thread;
  60                 thread = 0;
  61         }
  62
  63         return result;
  64 }
  65
  66 const list<Route *> &TrainRoutePlanner::get_routes_for(const Train &train) const
  67 {
  68         return get_train_info(train).routes;
  69 }
  70
  71 const list<TrainRouter::SequencePoint> &TrainRoutePlanner::get_sequence_for(const Train &train) const
  72 {
  73         return get_train_info(train).sequence;
  74 }
  75
  76 const TrainRoutePlanner::TrainRoutingInfo &TrainRoutePlanner::get_train_info(const Train &train) const
  77 {
  78         for(vector<TrainRoutingInfo>::const_iterator i=routed_trains.begin(); i!=routed_trains.end(); ++i)
  79                 if(i->train==&train)
  80                         return *i;
  81
  82         throw key_error(train.get_name());
  83 }
  84
  85 const TrainRoutePlanner::RoutingStep &TrainRoutePlanner::get_step()
  86 {
  87         steps.splice(steps.end(), queue, queue.begin());
  88         return steps.back();
  89 }
  90
  91 void TrainRoutePlanner::prepare_plan()
  92 {
  93         steps.clear();
  94         queue.clear();
  95         goal = 0;
  96         result = PENDING;
  97
  98         queue.push_back(RoutingStep());
  99         RoutingStep &start = queue.back();
 100         for(vector<TrainRoutingInfo>::iterator i=routed_trains.begin(); i!=routed_trains.end(); ++i)
 101                 start.trains.push_back(TrainRoutingState(*i));
 102         start.update_estimate();
 103 }
 104
 105 void TrainRoutePlanner::create_plan()
 106 {
 107         while(!queue.empty())
 108         {
 109                 const RoutingStep &step = get_step();
 110                 if(step.is_goal())
 111                 {
 112                         goal = &step;
 113                         return;
 114                 }
 115
 116                 add_steps(step);
 117         }
 118
 119         result = FAILED;
 120 }
 121
 122 void TrainRoutePlanner::add_steps(const RoutingStep &step)
 123 {
 124         list<RoutingStep> new_steps;
 125         step.create_successors(new_steps);
 126         new_steps.sort();
 127         queue.merge(new_steps);
 128 }
 129
 130 void TrainRoutePlanner::finalize_plan()
 131 {
 132         for(vector<TrainRoutingInfo>::iterator i=routed_trains.begin(); i!=routed_trains.end(); ++i)
 133         {
 134                 i->routes.clear();
 135                 i->sequence.clear();
 136                 for(unsigned j=0; j<2; ++j)
 137                         i->track_history[j] = 0;
 138         }
 139
 140         map<Track *, TrainRouter::SequencePoint *> sequenced_tracks;
 141         unsigned sequence = steps.size();
 142         for(const RoutingStep *i=goal; i; i=i->prev)
 143                 for(vector<TrainRoutingState>::const_iterator j=i->trains.begin(); j!=i->trains.end(); ++j)
 144                 {
 145                         Track **history = j->info->track_history;
 146                         if(j->track.track()==history[0])
 147                                 continue;
 148
 149                         Route *route = 0;
 150                         bool start_new_route = true;
 151                         if(!j->info->routes.empty())
 152                         {
 153                                 route = j->info->routes.front();
 154                                 start_new_route = route->has_track(*j->track);
 155                                 if(!start_new_route)
 156                                 {
 157                                         unsigned nls = j->track->get_n_link_slots();
 158                                         for(unsigned k=0; (!start_new_route && k<nls); ++k)
 159                                         {
 160                                                 Track *link = j->track->get_link(k);
 161                                                 start_new_route = (link && link!=history[0] && route->has_track(*link));
 162                                         }
 163                                 }
 164                         }
 165
 166                         if(start_new_route)
 167                         {
 168                                 route = new Route(j->info->train->get_layout());
 169                                 route->set_name("Router");
 170                                 route->set_temporary(true);
 171                                 for(unsigned k=0; (k<2 && history[k]); ++k)
 172                                         route->add_track(*history[k]);
 173                                 j->info->routes.push_front(route);
 174                         }
 175
 176                         route->add_track(*j->track.track());
 177                         history[1] = history[0];
 178                         history[0] = j->track.track();
 179
 180                         bool waitable = j->track.endpoint().paths!=j->track->get_type().get_paths();
 181                         map<Track *, TrainRouter::SequencePoint *>::iterator k = sequenced_tracks.find(j->track.track());
 182                         if(k!=sequenced_tracks.end())
 183                         {
 184                                 if(!k->second->preceding_train)
 185                                 {
 186                                         k->second->preceding_train = j->info->train;
 187                                         k->second->sequence_in = sequence;
 188                                 }
 189                                 j->info->sequence.push_front(TrainRouter::SequencePoint(j->track->get_block(), sequence));
 190                                 if(waitable)
 191                                         k->second = &j->info->sequence.front();
 192                                 --sequence;
 193                         }
 194                         else if(waitable)
 195                         {
 196                                 j->info->sequence.push_front(TrainRouter::SequencePoint(j->track->get_block(), sequence));
 197                                 sequenced_tracks[j->track.track()] = &j->info->sequence.front();
 198                                 --sequence;
 199                         }
 200                 }
 201
 202         result = COMPLETE;
 203 }
 204
 205
 206 TrainRoutePlanner::TrainRoutingInfo::TrainRoutingInfo(Train &t):
 207         train(&t),
 208         speed(train->get_maximum_speed()),
 209         router(train->get_ai_of_type<TrainRouter>())
 210 {
 211         // If no maximum speed is specified, use a sensible default
 212         if(!speed)
 213                 speed = 20*train->get_layout().get_catalogue().get_scale();
 214 }
 215
 216
 217 TrainRoutePlanner::OccupiedTrack::OccupiedTrack(Track &t, unsigned p, OccupiedTrack *n):
 218         track(&t),
 219         path_length(track->get_type().get_path_length(p)),
 220         next(n),
 221         n_tracks(next ? next->n_tracks+1 : 1),
 222         refcount(1)
 223 {
 224         if(next)
 225                 ++next->refcount;
 226 }
 227
 228 TrainRoutePlanner::OccupiedTrack::OccupiedTrack(const OccupiedTrack &other):
 229         track(other.track),
 230         path_length(other.path_length),
 231         next(other.next),
 232         n_tracks(other.n_tracks),
 233         refcount(1)
 234 {
 235         if(next)
 236                 ++next->refcount;
 237 }
 238
 239 TrainRoutePlanner::OccupiedTrack::~OccupiedTrack()
 240 {
 241         if(next && !--next->refcount)
 242                 delete next;
 243 }
 244
 245
 246 TrainRoutePlanner::TrainRoutingState::TrainRoutingState(TrainRoutingInfo &inf):
 247         info(&inf),
 248         occupied_tracks(0),
 249         state(MOVING),
 250         delay(info->router->get_departure_delay()),
 251         duration(info->router->get_trip_duration()),
 252         waypoint(info->router->get_n_waypoints() ? 0 : -1),
 253         blocked_by(-1)
 254 {
 255         const Vehicle *veh = &info->train->get_vehicle(0);
 256         // TODO margins
 257         TrackOffsetIter track_and_offs = veh->get_placement().get_position(VehiclePlacement::FRONT_BUFFER);
 258         track = track_and_offs.track_iter();
 259         offset = track_and_offs.offset();
 260         path = track->get_active_path();
 261
 262         while(Vehicle *next = veh->get_link(1))
 263                 veh = next;
 264         track_and_offs = veh->get_placement().get_position(VehiclePlacement::BACK_BUFFER);
 265         back_offset = track_and_offs.offset();
 266
 267         TrackIter iter = track_and_offs.track_iter();
 268         while(1)
 269         {
 270                 occupied_tracks = new OccupiedTrack(*iter, iter->get_active_path(), occupied_tracks);
 271                 if(iter.track()==track.track())
 272                         break;
 273                 iter = iter.next();
 274         }
 275
 276         update_estimate();
 277 }
 278
 279 TrainRoutePlanner::TrainRoutingState::TrainRoutingState(const TrainRoutingState &other):
 280         info(other.info),
 281         track(other.track),
 282         path(other.path),
 283         occupied_tracks(other.occupied_tracks),
 284         offset(other.offset),
 285         back_offset(other.back_offset),
 286         state(other.state),
 287         delay(other.delay),
 288         duration(other.duration),
 289         waypoint(other.waypoint),
 290         distance_traveled(other.distance_traveled),
 291         remaining_estimate(other.remaining_estimate),
 292         wait_time(other.wait_time),
 293         blocked_by(other.blocked_by)
 294 {
 295         ++occupied_tracks->refcount;
 296 }
 297
 298 TrainRoutePlanner::TrainRoutingState::~TrainRoutingState()
 299 {
 300         if(occupied_tracks && !--occupied_tracks->refcount)
 301                 delete occupied_tracks;
 302 }
 303
 304 Time::TimeDelta TrainRoutePlanner::TrainRoutingState::get_time_to_next_track() const
 305 {
 306         return ((track->get_type().get_path_length(path)-offset)/info->speed)*Time::sec+delay;
 307 }
 308
 309 bool TrainRoutePlanner::TrainRoutingState::is_occupying(Track &trk) const
 310 {
 311         if(state==ARRIVED && !duration && info->router->get_trip_duration())
 312                 return false;
 313
 314         OccupiedTrack *occ = occupied_tracks;
 315         for(unsigned n=occ->n_tracks; n>0; --n, occ=occ->next)
 316                 if(occ->track==&trk)
 317                         return true;
 318         return false;
 319 }
 320
 321 bool TrainRoutePlanner::TrainRoutingState::check_arrival()
 322 {
 323         TrainRouter &router = *info->router;
 324         TrackIter next_track = track.next(path);
 325
 326         if(waypoint<0 && router.is_destination(*track) && !router.is_destination(*next_track))
 327         {
 328                 state = ARRIVED;
 329                 return true;
 330         }
 331         else if(waypoint>=0 && router.is_waypoint(waypoint, *track) && !router.is_waypoint(waypoint, *next_track))
 332         {
 333                 ++waypoint;
 334                 if(waypoint>=static_cast<int>(router.get_n_waypoints()))
 335                         waypoint = -1;
 336         }
 337
 338         return false;
 339 }
 340
 341 void TrainRoutePlanner::TrainRoutingState::advance(float distance)
 342 {
 343         offset += distance;
 344         back_offset += distance;
 345
 346         OccupiedTrack *last_occ = occupied_tracks;
 347         for(unsigned n=occupied_tracks->n_tracks; n>1; --n)
 348                 last_occ = last_occ->next;
 349
 350         // XXX What if there's multiple tracks to remove?
 351         if(back_offset>last_occ->path_length)
 352         {
 353                 back_offset -= last_occ->path_length;
 354                 if(occupied_tracks->refcount>1)
 355                 {
 356                         --occupied_tracks->refcount;
 357                         occupied_tracks = new OccupiedTrack(*occupied_tracks);
 358                 }
 359                 --occupied_tracks->n_tracks;
 360         }
 361
 362         distance_traveled += distance;
 363         remaining_estimate -= distance;
 364 }
 365
 366 void TrainRoutePlanner::TrainRoutingState::advance(const Time::TimeDelta &dt)
 367 {
 368         if(delay>=dt)
 369         {
 370                 delay -= dt;
 371                 return;
 372         }
 373
 374         float secs = dt/Time::sec;
 375         if(delay)
 376         {
 377                 secs -= delay/Time::sec;
 378                 delay = Time::zero;
 379         }
 380
 381         if(duration)
 382                 duration = max(duration-secs*Time::sec, Time::zero);
 383
 384         if(state==MOVING)
 385                 advance(info->speed*secs);
 386         else if(state!=ARRIVED)
 387                 wait_time += secs*Time::sec;
 388 }
 389
 390 void TrainRoutePlanner::TrainRoutingState::advance_track(unsigned next_path)
 391 {
 392         float distance = occupied_tracks->path_length-offset;
 393         track = track.next(path);
 394         path = next_path;
 395         occupied_tracks = new OccupiedTrack(*track, path, occupied_tracks);
 396         advance(distance);
 397         offset = 0;
 398 }
 399
 400 void TrainRoutePlanner::TrainRoutingState::update_estimate()
 401 {
 402         TrackIter iter = track.reverse(path);
 403         float distance = info->router->get_metric(waypoint).get_distance_from(*iter.track(), iter.entry());
 404         distance += track->get_type().get_path_length(path)-offset;
 405         remaining_estimate = distance;
 406 }
 407
 408
 409 TrainRoutePlanner::RoutingStep::RoutingStep():
 410         prev(0)
 411 { }
 412
 413 TrainRoutePlanner::RoutingStep::RoutingStep(const RoutingStep *p):
 414         time(p->time),
 415         cost_estimate(p->cost_estimate),
 416         trains(p->trains),
 417         prev(p)
 418 { }
 419
 420 void TrainRoutePlanner::RoutingStep::create_successors(list<RoutingStep> &new_steps) const
 421 {
 422         RoutingStep next(this);
 423         if(next.update_states())
 424         {
 425                 if(next.check_deadlocks())
 426                         return;
 427
 428                 new_steps.push_back(next);
 429                 return;
 430         }
 431
 432         int train_index = find_next_train();
 433         if(train_index<0)
 434                 return;
 435
 436         TrainRoutingState &train = next.trains[train_index];
 437
 438         Time::TimeDelta dt = train.get_time_to_next_track();
 439         next.advance(dt);
 440
 441         if(train.check_arrival())
 442         {
 443                 new_steps.push_back(next);
 444                 return;
 445         }
 446
 447         TrackIter next_track = train.track.next(train.path);
 448         train.advance_track(0);
 449
 450         const TrackType::Endpoint &next_entry_ep = next_track.endpoint();
 451         for(unsigned i=0; next_entry_ep.paths>>i; ++i)
 452                 if(next_entry_ep.has_path(i))
 453                 {
 454                         train.path = i;
 455                         train.update_estimate();
 456                         next.update_estimate();
 457                         if(next.is_viable())
 458                                 new_steps.push_back(next);
 459                 }
 460
 461         if(next_entry_ep.paths!=next_track->get_type().get_paths())
 462         {
 463                 RoutingStep wait(this);
 464                 wait.advance(dt);
 465                 wait.trains[train_index].state = WAITING;
 466                 if(wait.is_viable())
 467                         new_steps.push_back(wait);
 468         }
 469 }
 470
 471 bool TrainRoutePlanner::RoutingStep::update_states()
 472 {
 473         bool changes = false;
 474         for(vector<TrainRoutingState>::iterator i=trains.begin(); i!=trains.end(); ++i)
 475         {
 476                 if(i->state==ARRIVED)
 477                         continue;
 478
 479                 TrainState old_state = i->state;
 480
 481                 TrackIter next_track = i->track.next(i->path);
 482                 if(next_track)
 483                 {
 484                         i->blocked_by = get_occupant(*next_track);
 485                         if(i->blocked_by>=0)
 486                                 i->state = BLOCKED;
 487                         else if(i->state==BLOCKED)
 488                                 i->state = MOVING;
 489                 }
 490                 else
 491                         i->state = BLOCKED;
 492
 493                 if(i->state!=old_state)
 494                         changes = true;
 495         }
 496
 497         return changes;
 498 }
 499
 500 bool TrainRoutePlanner::RoutingStep::check_deadlocks() const
 501 {
 502         for(vector<TrainRoutingState>::const_iterator i=trains.begin(); i!=trains.end(); ++i)
 503         {
 504                 if(i->state!=BLOCKED)
 505                         continue;
 506
 507                 if(i->blocked_by<0)
 508                         return true;
 509
 510                 int slow = i->blocked_by;
 511                 int fast = trains[slow].blocked_by;
 512                 while(fast>=0 && trains[fast].blocked_by>=0)
 513                 {
 514                         if(fast==slow)
 515                                 return true;
 516
 517                         slow = trains[slow].blocked_by;
 518                         fast = trains[trains[fast].blocked_by].blocked_by;
 519                 }
 520         }
 521
 522         return false;
 523 }
 524
 525 int TrainRoutePlanner::RoutingStep::get_occupant(Track &track) const
 526 {
 527         for(unsigned i=0; i<trains.size(); ++i)
 528                 if(trains[i].is_occupying(track))
 529                         return i;
 530
 531         return -1;
 532 }
 533
 534 int TrainRoutePlanner::RoutingStep::find_next_train() const
 535 {
 536         Time::TimeDelta min_dt;
 537         int next_train = -1;
 538         for(unsigned i=0; i<trains.size(); ++i)
 539                 if(trains[i].state==MOVING)
 540                 {
 541                         Time::TimeDelta dt = trains[i].get_time_to_next_track();
 542                         if(dt<min_dt || next_train<0)
 543                         {
 544                                 min_dt = dt;
 545                                 next_train = i;
 546                         }
 547                 }
 548
 549         return next_train;
 550 }
 551
 552 void TrainRoutePlanner::RoutingStep::advance(const Time::TimeDelta &dt)
 553 {
 554         time += dt;
 555         for(vector<TrainRoutingState>::iterator i=trains.begin(); i!=trains.end(); ++i)
 556                 i->advance(dt);
 557 }
 558
 559 void TrainRoutePlanner::RoutingStep::update_estimate()
 560 {
 561         cost_estimate = Time::zero;
 562         for(vector<TrainRoutingState>::const_iterator i=trains.begin(); i!=trains.end(); ++i)
 563                 if(i->remaining_estimate>=0)
 564                         cost_estimate += i->wait_time+((i->distance_traveled+i->remaining_estimate)/i->info->speed)*Time::sec;
 565 }
 566
 567 bool TrainRoutePlanner::RoutingStep::is_viable() const
 568 {
 569         for(vector<TrainRoutingState>::const_iterator i=trains.begin(); i!=trains.end(); ++i)
 570                 if(i->remaining_estimate<0)
 571                         return false;
 572
 573         for(vector<TrainRoutingState>::const_iterator i=trains.begin(); i!=trains.end(); ++i)
 574                 if(i->state==MOVING)
 575                         return true;
 576
 577         return false;
 578 }
 579
 580 bool TrainRoutePlanner::RoutingStep::is_goal() const
 581 {
 582         for(vector<TrainRoutingState>::const_iterator i=trains.begin(); i!=trains.end(); ++i)
 583                 if(i->state!=ARRIVED)
 584                         return false;
 585         return true;
 586 }
 587
 588 bool TrainRoutePlanner::RoutingStep::operator<(const RoutingStep &other) const
 589 {
 590         return cost_estimate<other.cost_estimate;
 591 }
 592
 593
 594 TrainRoutePlanner::PlanningThread::PlanningThread(TrainRoutePlanner &p):
 595         planner(p)
 596 {
 597         launch();
 598 }
 599
 600 void TrainRoutePlanner::PlanningThread::main()
 601 {
 602         planner.create_plan();
 603 }
 604
 605 } // namespace R2C2