1 /******************************************************************************
3 * Written by Durk Talsma, started May 5, 2004.
5 * This program is free software; you can redistribute it and/or
6 * modify it under the terms of the GNU General Public License as
7 * published by the Free Software Foundation; either version 2 of the
8 * License, or (at your option) any later version.
10 * This program is distributed in the hope that it will be useful, but
11 * WITHOUT ANY WARRANTY; without even the implied warranty of
12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13 * General Public License for more details.
15 * You should have received a copy of the GNU General Public License
16 * along with this program; if not, write to the Free Software
17 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
20 ****************************************************************************
22 *****************************************************************************/
40 #include <simgear/compiler.h>
41 #include <simgear/sg_inlines.h>
42 #include <simgear/math/sg_geodesy.hxx>
43 #include <simgear/props/props.hxx>
44 #include <simgear/route/waypoint.hxx>
45 #include <simgear/structure/subsystem_mgr.hxx>
46 #include <simgear/xml/easyxml.hxx>
48 #include <AIModel/AIFlightPlan.hxx>
49 #include <AIModel/AIManager.hxx>
50 #include <AIModel/AIAircraft.hxx>
51 #include <Airports/simple.hxx>
52 #include <Main/fg_init.hxx> // That's pretty ugly, but I need fgFindAirportID
55 #include "SchedFlight.hxx"
56 #include "TrafficMgr.hxx"
58 /******************************************************************************
59 * the FGAISchedule class contains data members and code to maintain a
60 * schedule of Flights for an articically controlled aircraft.
61 *****************************************************************************/
62 FGAISchedule::FGAISchedule()
77 FGAISchedule::FGAISchedule(string mdl,
88 FGScheduledFlightVec flt)*/
89 FGAISchedule::FGAISchedule(string model,
106 flightIdentifier = flightId;
115 /*for (FGScheduledFlightVecIterator i = flt.begin();
118 flights.push_back(new FGScheduledFlight((*(*i))));*/
129 FGAISchedule::FGAISchedule(const FGAISchedule &other)
131 modelPath = other.modelPath;
132 homePort = other.homePort;
133 livery = other.livery;
134 registration = other.registration;
136 flightIdentifier = other.flightIdentifier;
137 flights = other.flights;
138 AIManagerRef = other.AIManagerRef;
139 acType = other.acType;
140 airline = other.airline;
141 m_class = other.m_class;
142 firstRun = other.firstRun;
143 radius = other.radius;
144 groundOffset = other.groundOffset;
145 flightType = other.flightType;
147 distanceToUser = other.distanceToUser;
148 currentDestination = other.currentDestination;
149 firstRun = other.firstRun;
150 runCount = other.runCount;
152 lastRun = other.lastRun;
153 initialized = other.initialized;
159 FGAISchedule::~FGAISchedule()
161 /* for (FGScheduledFlightVecIterator flt = flights.begin(); flt != flights.end(); flt++)
168 bool FGAISchedule::init()
171 //SGTime* currTimeDate = globals->get_time_params();
173 //tm *temp = currTimeDate->getGmt();
175 //sgTimeFormatTime(&targetTimeDate, buffer);
176 //cout << "Scheduled Time " << buffer << endl;
177 //cout << "Time :" << time(NULL) << " SGTime : " << sgTimeGetGMT(temp) << endl;
178 /*for (FGScheduledFlightVecIterator i = flights.begin();
182 //i->adjustTime(now);
183 if (!((*i)->initializeAirports()))
186 //sort(flights.begin(), flights.end());
187 // Since time isn't initialized yet when this function is called,
188 // Find the closest possible airport.
189 // This should give a reasonable initialization order.
190 //setClosestDistanceToUser();
194 bool FGAISchedule::update(time_t now, const SGVec3d& userCart)
199 remainingTimeEnroute,
204 scheduleFlights(now);
205 if (flights.empty()) { // No flights available for this aircraft
210 // Sort all the scheduled flights according to scheduled departure time.
211 // Because this is done at every update, we only need to check the status
212 // of the first listed flight.
213 //sort(flights.begin(), flights.end(), compareScheduledFlights);
216 if (fgGetBool("/sim/traffic-manager/instantaneous-action") == true) {
217 deptime = now; // + rand() % 300; // Wait up to 5 minutes until traffic starts moving to prevent too many aircraft
218 // from cluttering the gate areas.
223 FGScheduledFlight* flight = flights.front();
225 deptime = flight->getDepartureTime();
226 //cerr << "Settiing departure time " << deptime << endl;
230 // Check if this aircraft has been released.
231 FGTrafficManager *tmgr = (FGTrafficManager *) globals->get_subsystem("traffic-manager");
232 if (tmgr->isReleased(AIManagerRef)) {
235 return true; // in visual range, let the AIManager handle it
239 // This flight entry is entirely in the past, do we need to
240 // push it forward in time to the next scheduled departure.
241 if (flight->getArrivalTime() < now) {
242 SG_LOG (SG_GENERAL, SG_BULK, "Traffic Manager: Flight is in the Past");
243 // Don't just update: check whether we need to load a new leg. etc.
244 // This update occurs for distant aircraft, so we can update the current leg
245 // and detach it from the current list of aircraft.
247 flights.erase(flights.begin()); // pop_front(), effectively
251 FGAirport* dep = flight->getDepartureAirport();
252 FGAirport* arr = flight->getArrivalAirport();
257 double speed = 450.0;
259 totalTimeEnroute = flight->getArrivalTime() - flight->getDepartureTime();
260 if (flight->getDepartureTime() < now) {
261 elapsedTimeEnroute = now - flight->getDepartureTime();
262 remainingTimeEnroute = totalTimeEnroute - elapsedTimeEnroute;
263 double x = elapsedTimeEnroute / (double) totalTimeEnroute;
265 // current pos is based on great-circle course between departure/arrival,
266 // with percentage of distance travelled, based upon percentage of time
268 double course, az2, distanceM;
269 SGGeodesy::inverse(dep->geod(), arr->geod(), course, az2, distanceM);
270 double coveredDistance = distanceM * x;
272 SGGeodesy::direct(dep->geod(), course, coveredDistance, position, az2);
274 SG_LOG (SG_GENERAL, SG_BULK, "Traffic Manager: Flight is in progress, %=" << x);
275 speed = ((distanceM - coveredDistance) * SG_METER_TO_NM) / 3600.0;
278 remainingTimeEnroute = totalTimeEnroute;
279 elapsedTimeEnroute = 0;
280 position = dep->geod();
281 SG_LOG (SG_GENERAL, SG_BULK, "Traffic Manager: Flight is pending, departure in "
282 << flight->getDepartureTime() - now << " seconds ");
285 // departure / arrival coincident
286 remainingTimeEnroute = totalTimeEnroute = 0.0;
287 elapsedTimeEnroute = 0;
288 position = dep->geod();
291 // cartesian calculations are more numerically stable over the (potentially)
292 // large distances involved here: see bug #80
293 distanceToUser = dist(userCart, SGVec3d::fromGeod(position)) * SG_METER_TO_NM;
295 // If distance between user and simulated aircaft is less
296 // then 500nm, create this flight. At jet speeds 500 nm is roughly
297 // one hour flight time, so that would be a good approximate point
298 // to start a more detailed simulation of this aircraft.
299 SG_LOG (SG_GENERAL, SG_BULK, "Traffic manager: " << registration << " is scheduled for a flight from "
300 << dep->getId() << " to " << arr->getId() << ". Current distance to user: "
302 if (distanceToUser >= TRAFFICTOAIDISTTOSTART) {
303 return true; // out of visual range, for the moment.
305 return createAIAircraft(flight, speed, deptime);
308 bool FGAISchedule::createAIAircraft(FGScheduledFlight* flight, double speedKnots, time_t deptime)
310 FGAirport* dep = flight->getDepartureAirport();
311 FGAirport* arr = flight->getArrivalAirport();
312 string flightPlanName = dep->getId() + "-" + arr->getId() + ".xml";
313 SG_LOG(SG_GENERAL, SG_INFO, "Traffic manager: Creating AIModel from:" << flightPlanName);
315 // Only allow traffic to be created when the model path (or the AI version of mp) exists
316 SGPath mp(globals->get_fg_root());
319 mp.append(modelPath);
321 mp_ai.append(modelPath);
323 if (!mp.exists() && !mp_ai.exists()) {
324 SG_LOG(SG_GENERAL, SG_WARN, "TrafficManager: Could not load model " << mp.str());
328 FGAIAircraft *aircraft = new FGAIAircraft(this);
329 aircraft->setPerformance(m_class); //"jet_transport";
330 aircraft->setCompany(airline); //i->getAirline();
331 aircraft->setAcType(acType); //i->getAcType();
332 aircraft->setPath(modelPath.c_str());
333 //aircraft->setFlightPlan(flightPlanName);
334 aircraft->setLatitude(position.getLatitudeDeg());
335 aircraft->setLongitude(position.getLongitudeDeg());
336 aircraft->setAltitude(flight->getCruiseAlt()*100); // convert from FL to feet
337 aircraft->setSpeed(0);
338 aircraft->setBank(0);
340 courseToDest = SGGeodesy::courseDeg(position, arr->geod());
341 FGAIFlightPlan *fp = new FGAIFlightPlan(aircraft, flightPlanName, courseToDest, deptime,
342 dep, arr, true, radius,
343 flight->getCruiseAlt()*100,
344 position.getLatitudeDeg(),
345 position.getLongitudeDeg(),
346 speedKnots, flightType, acType,
348 if (fp->isValidPlan()) {
349 aircraft->SetFlightPlan(fp);
350 FGAIManager* aimgr = (FGAIManager *) globals-> get_subsystem("ai-model");
351 aimgr->attach(aircraft);
352 AIManagerRef = aircraft->getID();
357 //hand back the flights that had already been scheduled
358 while (!flights.empty()) {
359 flights.front()->release();
360 flights.erase(flights.begin());
366 // Create an initial heading for user controlled aircraft.
367 void FGAISchedule::setHeading()
369 courseToDest = SGGeodesy::courseDeg((*flights.begin())->getDepartureAirport()->geod(), (*flights.begin())->getArrivalAirport()->geod());
372 void FGAISchedule::scheduleFlights(time_t now)
374 if (!flights.empty()) {
377 //string startingPort;
378 string userPort = fgGetString("/sim/presets/airport-id");
379 SG_LOG(SG_GENERAL, SG_BULK, "Scheduling Flights for : " << modelPath << " " << registration << " " << homePort);
380 FGScheduledFlight *flight = NULL;
382 if (currentDestination.empty()) {
383 flight = findAvailableFlight(userPort, flightIdentifier, now, (now+6400));
385 flight = findAvailableFlight(currentDestination, flightIdentifier);
387 flight = findAvailableFlight(currentDestination, flightIdentifier);
392 currentDestination = flight->getArrivalAirport()->getId();
393 //cerr << "Current destination " << currentDestination << endl;
395 string departurePort = flight->getDepartureAirport()->getId();
396 if (userPort == departurePort) {
407 dep = flight->getDepartureTime();
408 arr = flight->getArrivalTime();
409 string depT = asctime(gmtime(&dep));
410 string arrT = asctime(gmtime(&arr));
412 depT = depT.substr(0,24);
413 arrT = arrT.substr(0,24);
414 SG_LOG(SG_GENERAL, SG_BULK, " Flight " << flight->getCallSign() << ":"
415 << " " << flight->getDepartureAirport()->getId() << ":"
416 << " " << depT << ":"
417 << " \"" << flight->getArrivalAirport()->getId() << "\"" << ":"
418 << " " << arrT << ":");
420 flights.push_back(flight);
421 } while (currentDestination != homePort);
422 SG_LOG(SG_GENERAL, SG_BULK, " Done ");
425 bool FGAISchedule::next()
427 if (!flights.empty()) {
428 flights.front()->release();
429 flights.erase(flights.begin());
432 FGScheduledFlight *flight = findAvailableFlight(currentDestination, flightIdentifier);
437 currentDestination = flight->getArrivalAirport()->getId();
440 dep = flight->getDepartureTime();
441 arr = flight->getArrivalTime();
442 string depT = asctime(gmtime(&dep));
443 string arrT = asctime(gmtime(&arr));
445 depT = depT.substr(0,24);
446 arrT = arrT.substr(0,24);
447 //cerr << " " << flight->getCallSign() << ":"
448 // << " " << flight->getDepartureAirport()->getId() << ":"
449 // << " " << depT << ":"
450 // << " \"" << flight->getArrivalAirport()->getId() << "\"" << ":"
451 // << " " << arrT << ":" << endl;
453 flights.push_back(flight);
457 FGScheduledFlight* FGAISchedule::findAvailableFlight (const string ¤tDestination,
459 time_t min, time_t max)
461 time_t now = time(NULL) + fgGetLong("/sim/time/warp");
463 FGTrafficManager *tmgr = (FGTrafficManager *) globals->get_subsystem("traffic-manager");
464 FGScheduledFlightVecIterator fltBegin, fltEnd;
465 fltBegin = tmgr->getFirstFlight(req);
466 fltEnd = tmgr->getLastFlight(req);
469 //cerr << "Finding available flight " << endl;
471 // Traverse every registered flight
472 if (fltBegin == fltEnd) {
473 //cerr << "No Flights Scheduled for " << req << endl;
476 for (FGScheduledFlightVecIterator i = fltBegin; i != fltEnd; i++) {
477 (*i)->adjustTime(now);
478 //sort(fltBegin, fltEnd, compareScheduledFlights);
479 //cerr << counter++ << endl;
481 std::sort(fltBegin, fltEnd, compareScheduledFlights);
482 for (FGScheduledFlightVecIterator i = fltBegin; i != fltEnd; i++) {
485 if (!(*i)->isAvailable()) {
486 //cerr << (*i)->getCallSign() << "is no longer available" << endl;
489 if (!((*i)->getRequirement() == req)) {
492 if (!(((*i)->getArrivalAirport()) && ((*i)->getDepartureAirport()))) {
495 if (!(currentDestination.empty())) {
496 if (currentDestination != (*i)->getDepartureAirport()->getId()) {
497 //cerr << (*i)->getCallSign() << "Doesn't match destination" << endl;
498 //cerr << "Current Destination " << currentDestination << "Doesnt match flight's " <<
499 // (*i)->getArrivalAirport()->getId() << endl;
503 if (flights.size()) {
504 time_t arrival = flights.back()->getArrivalTime();
505 int groundTime = groundTimeFromRadius();
506 if ((*i)->getDepartureTime() < (arrival+(groundTime)))
510 time_t dep = (*i)->getDepartureTime();
511 if ((dep < min) || (dep > max))
515 // So, if we actually get here, we have a winner
516 //cerr << "found flight: " << req << " : " << currentDestination << " : " <<
517 // (*i)->getArrivalAirport()->getId() << endl;
522 // if currentDestination has a value, does it match departure of next flight?
523 // is departure time later than planned arrival?
524 // is departure port valid?
525 // is arrival port valid?
526 //cerr << "Ack no flight found: " << endl;
530 int FGAISchedule::groundTimeFromRadius()
534 else if (radius < 15)
536 else if (radius < 20)
538 else if (radius < 25)
540 else if (radius < 30)
547 double FGAISchedule::getSpeed()
549 FGScheduledFlightVecIterator i = flights.begin();
551 FGAirport* dep = (*i)->getDepartureAirport(),
552 *arr = (*i)->getArrivalAirport();
553 double dist = SGGeodesy::distanceNm(dep->geod(), arr->geod());
554 double remainingTimeEnroute = (*i)->getArrivalTime() - (*i)->getDepartureTime();
556 double speed = dist / (remainingTimeEnroute/3600.0);
557 SG_CLAMP_RANGE(speed, 300.0, 500.0);
561 void FGAISchedule::setScore ()
564 score = ((double) hits / (double) runCount);
566 if (homePort == fgGetString("/sim/presets/airport-id")) {
575 bool compareSchedules(FGAISchedule*a, FGAISchedule*b)
580 bool FGAISchedule::operator< (const FGAISchedule &other) const
582 //cerr << "Sorting " << registration << " and " << other.registration << endl;
583 double currentScore = score * (1.5 - lastRun);
584 double otherScore = other.score * (1.5 - other.lastRun);
585 return currentScore > otherScore;