Line data Source code
1 : /****************************************************************************/
2 : // Eclipse SUMO, Simulation of Urban MObility; see https://eclipse.dev/sumo
3 : // Copyright (C) 2001-2026 German Aerospace Center (DLR) and others.
4 : // This program and the accompanying materials are made available under the
5 : // terms of the Eclipse Public License 2.0 which is available at
6 : // https://www.eclipse.org/legal/epl-2.0/
7 : // This Source Code may also be made available under the following Secondary
8 : // Licenses when the conditions for such availability set forth in the Eclipse
9 : // Public License 2.0 are satisfied: GNU General Public License, version 2
10 : // or later which is available at
11 : // https://www.gnu.org/licenses/old-licenses/gpl-2.0-standalone.html
12 : // SPDX-License-Identifier: EPL-2.0 OR GPL-2.0-or-later
13 : /****************************************************************************/
14 : /// @file RONet.cpp
15 : /// @author Daniel Krajzewicz
16 : /// @author Jakob Erdmann
17 : /// @author Michael Behrisch
18 : /// @date Sept 2002
19 : ///
20 : // The router's network representation
21 : /****************************************************************************/
22 : #include <config.h>
23 :
24 : #include <algorithm>
25 : #include <utils/router/RouteCostCalculator.h>
26 : #include <utils/vehicle/SUMOVTypeParameter.h>
27 : #include <utils/router/SUMOAbstractRouter.h>
28 : #include <utils/options/OptionsCont.h>
29 : #include <utils/common/UtilExceptions.h>
30 : #include <utils/common/ToString.h>
31 : #include <utils/common/StringUtils.h>
32 : #include <utils/common/RandHelper.h>
33 : #include <utils/common/SUMOVehicleClass.h>
34 : #include <utils/iodevices/OutputDevice.h>
35 : #include "ROEdge.h"
36 : #include "ROLane.h"
37 : #include "RONode.h"
38 : #include "ROPerson.h"
39 : #include "RORoute.h"
40 : #include "RORouteDef.h"
41 : #include "ROVehicle.h"
42 : #include "ROAbstractEdgeBuilder.h"
43 : #include "RONet.h"
44 :
45 :
46 : // ===========================================================================
47 : // static member definitions
48 : // ===========================================================================
49 : RONet* RONet::myInstance = nullptr;
50 :
51 :
52 : // ===========================================================================
53 : // method definitions
54 : // ===========================================================================
55 : RONet*
56 5902692 : RONet::getInstance(void) {
57 5902692 : if (myInstance != nullptr) {
58 5902692 : return myInstance;
59 : }
60 0 : throw ProcessError(TL("A network was not yet constructed."));
61 : }
62 :
63 :
64 4368 : RONet::RONet() :
65 4368 : myVehicleTypes(), myDefaultVTypeMayBeDeleted(true),
66 4368 : myDefaultPedTypeMayBeDeleted(true),
67 4368 : myDefaultBikeTypeMayBeDeleted(true),
68 4368 : myDefaultTaxiTypeMayBeDeleted(true),
69 4368 : myDefaultRailTypeMayBeDeleted(true),
70 4368 : myHaveActiveFlows(true),
71 4368 : myRoutesOutput(nullptr), myRouteAlternativesOutput(nullptr), myTypesOutput(nullptr),
72 4368 : myReadRouteNo(0), myDiscardedRouteNo(0), myWrittenRouteNo(0),
73 4368 : myHavePermissions(false),
74 4368 : myHaveParamRestrictions(false),
75 4368 : myNumInternalEdges(0),
76 4368 : myErrorHandler(OptionsCont::getOptions().exists("ignore-errors")
77 8289 : && OptionsCont::getOptions().getBool("ignore-errors") ? MsgHandler::getWarningInstance() : MsgHandler::getErrorInstance()),
78 4368 : myKeepVTypeDist(OptionsCont::getOptions().exists("keep-vtype-distributions")
79 8289 : && OptionsCont::getOptions().getBool("keep-vtype-distributions")),
80 4368 : myDoPTRouting(!OptionsCont::getOptions().exists("ptline-routing")
81 7823 : || OptionsCont::getOptions().getBool("ptline-routing")),
82 4368 : myKeepFlows(OptionsCont::getOptions().exists("keep-flows")
83 7823 : && OptionsCont::getOptions().getBool("keep-flows")),
84 4368 : myHasBidiEdges(false),
85 15646 : myMaxTraveltime(OptionsCont::getOptions().exists("max-traveltime") ? STEPS2TIME(string2time(OptionsCont::getOptions().getString("max-traveltime"))) : -1) {
86 4368 : if (myInstance != nullptr) {
87 0 : throw ProcessError(TL("A network was already constructed."));
88 : }
89 4368 : SUMOVTypeParameter* type = new SUMOVTypeParameter(DEFAULT_VTYPE_ID, SVC_PASSENGER);
90 4368 : type->onlyReferenced = true;
91 4368 : myVehicleTypes.add(type->id, type);
92 :
93 4368 : SUMOVTypeParameter* defPedType = new SUMOVTypeParameter(DEFAULT_PEDTYPE_ID, SVC_PEDESTRIAN);
94 4368 : defPedType->onlyReferenced = true;
95 4368 : defPedType->parametersSet |= VTYPEPARS_VEHICLECLASS_SET;
96 4368 : myVehicleTypes.add(defPedType->id, defPedType);
97 :
98 4368 : SUMOVTypeParameter* defBikeType = new SUMOVTypeParameter(DEFAULT_BIKETYPE_ID, SVC_BICYCLE);
99 4368 : defBikeType->onlyReferenced = true;
100 4368 : defBikeType->parametersSet |= VTYPEPARS_VEHICLECLASS_SET;
101 4368 : myVehicleTypes.add(defBikeType->id, defBikeType);
102 :
103 4368 : SUMOVTypeParameter* defTaxiType = new SUMOVTypeParameter(DEFAULT_TAXITYPE_ID, SVC_TAXI);
104 4368 : defTaxiType->onlyReferenced = true;
105 4368 : defTaxiType->parametersSet |= VTYPEPARS_VEHICLECLASS_SET;
106 4368 : myVehicleTypes.add(defTaxiType->id, defTaxiType);
107 :
108 4368 : SUMOVTypeParameter* defRailType = new SUMOVTypeParameter(DEFAULT_RAILTYPE_ID, SVC_RAIL);
109 4368 : defRailType->onlyReferenced = true;
110 4368 : defRailType->parametersSet |= VTYPEPARS_VEHICLECLASS_SET;
111 4368 : myVehicleTypes.add(defRailType->id, defRailType);
112 :
113 4368 : myInstance = this;
114 4368 : }
115 :
116 :
117 8279 : RONet::~RONet() {
118 4571 : for (const auto& routables : myRoutables) {
119 426 : for (RORoutable* const r : routables.second) {
120 213 : const ROVehicle* const veh = dynamic_cast<const ROVehicle*>(r);
121 : // delete routes and the vehicle
122 213 : if (veh != nullptr && veh->getRouteDefinition()->getID()[0] == '!') {
123 204 : if (!myRoutes.remove(veh->getRouteDefinition()->getID())) {
124 76 : delete veh->getRouteDefinition();
125 : }
126 : }
127 213 : delete r;
128 : }
129 : }
130 4394 : for (const RORoutable* const r : myPTVehicles) {
131 36 : const ROVehicle* const veh = dynamic_cast<const ROVehicle*>(r);
132 : // delete routes and the vehicle
133 36 : if (veh != nullptr && veh->getRouteDefinition()->getID()[0] == '!') {
134 36 : if (!myRoutes.remove(veh->getRouteDefinition()->getID())) {
135 8 : delete veh->getRouteDefinition();
136 : }
137 : }
138 36 : delete r;
139 : }
140 : myRoutables.clear();
141 4442 : for (const auto& vTypeDist : myVTypeDistDict) {
142 168 : delete vTypeDist.second;
143 : }
144 30069 : }
145 :
146 :
147 : void
148 17 : RONet::addSpeedRestriction(const std::string& id, const SUMOVehicleClass svc, const double speed) {
149 17 : mySpeedRestrictions[id][svc] = speed;
150 17 : }
151 :
152 :
153 : double
154 1326 : RONet::getPreference(const std::string& routingType, const SUMOVTypeParameter& pars) const {
155 1326 : if (gRoutingPreferences) {
156 1326 : auto it = myVTypePreferences.find(pars.id);
157 1326 : if (it != myVTypePreferences.end()) {
158 : auto it2 = it->second.find(routingType);
159 90 : if (it2 != it->second.end()) {
160 24 : return it2->second;
161 : }
162 : }
163 : auto it3 = myVClassPreferences.find(pars.vehicleClass);
164 1302 : if (it3 != myVClassPreferences.end()) {
165 : auto it4 = it3->second.find(routingType);
166 180 : if (it4 != it3->second.end()) {
167 48 : return it4->second;
168 : }
169 : }
170 : // fallback to generel preferences
171 1254 : it = myVTypePreferences.find("");
172 1254 : if (it != myVTypePreferences.end()) {
173 : auto it2 = it->second.find(routingType);
174 798 : if (it2 != it->second.end()) {
175 144 : return it2->second;
176 : }
177 : }
178 : }
179 : return 1;
180 : }
181 :
182 :
183 : void
184 24 : RONet::addPreference(const std::string& routingType, SUMOVehicleClass svc, double prio) {
185 24 : myVClassPreferences[svc][routingType] = prio;
186 24 : gRoutingPreferences = true;
187 24 : }
188 :
189 :
190 : void
191 32 : RONet::addPreference(const std::string& routingType, std::string vType, double prio) {
192 32 : myVTypePreferences[vType][routingType] = prio;
193 32 : gRoutingPreferences = true;
194 32 : }
195 :
196 :
197 :
198 : const std::map<SUMOVehicleClass, double>*
199 66 : RONet::getRestrictions(const std::string& id) const {
200 : std::map<std::string, std::map<SUMOVehicleClass, double> >::const_iterator i = mySpeedRestrictions.find(id);
201 66 : if (i == mySpeedRestrictions.end()) {
202 : return nullptr;
203 : }
204 26 : return &i->second;
205 : }
206 :
207 :
208 : bool
209 331278 : RONet::addEdge(ROEdge* edge) {
210 331278 : if (!myEdges.add(edge->getID(), edge)) {
211 36 : WRITE_ERRORF(TL("The edge '%' occurs at least twice."), edge->getID());
212 12 : delete edge;
213 12 : return false;
214 : }
215 331266 : if (edge->isInternal()) {
216 161488 : myNumInternalEdges += 1;
217 : }
218 : return true;
219 : }
220 :
221 :
222 : bool
223 1818 : RONet::addDistrict(const std::string id, ROEdge* source, ROEdge* sink) {
224 : if (myDistricts.count(id) > 0) {
225 0 : WRITE_ERRORF(TL("The TAZ '%' occurs at least twice."), id);
226 0 : delete source;
227 0 : delete sink;
228 0 : return false;
229 : }
230 : sink->setFunction(SumoXMLEdgeFunc::CONNECTOR);
231 1818 : if (!addEdge(sink)) {
232 : return false;
233 : }
234 : source->setFunction(SumoXMLEdgeFunc::CONNECTOR);
235 1818 : if (!addEdge(source)) {
236 : return false;
237 : }
238 : sink->setOtherTazConnector(source);
239 : source->setOtherTazConnector(sink);
240 3636 : myDistricts[id] = std::make_pair(std::vector<std::string>(), std::vector<std::string>());
241 1818 : return true;
242 : }
243 :
244 :
245 : bool
246 1383 : RONet::addDistrictEdge(const std::string tazID, const std::string edgeID, const bool isSource) {
247 : if (myDistricts.count(tazID) == 0) {
248 0 : WRITE_ERRORF(TL("The TAZ '%' is unknown."), tazID);
249 0 : return false;
250 : }
251 : ROEdge* edge = getEdge(edgeID);
252 1383 : if (edge == nullptr) {
253 0 : WRITE_ERRORF(TL("The edge '%' for TAZ '%' is unknown."), edgeID, tazID);
254 0 : return false;
255 : }
256 1383 : if (isSource) {
257 2073 : getEdge(tazID + "-source")->addSuccessor(edge);
258 691 : myDistricts[tazID].first.push_back(edgeID);
259 : } else {
260 2076 : edge->addSuccessor(getEdge(tazID + "-sink"));
261 692 : myDistricts[tazID].second.push_back(edgeID);
262 : }
263 : return true;
264 : }
265 :
266 :
267 : void
268 101 : RONet::addJunctionTaz(ROAbstractEdgeBuilder& eb) {
269 1544 : for (auto item : myNodes) {
270 : const std::string tazID = item.first;
271 8 : if (myDistricts.count(tazID) != 0) {
272 24 : WRITE_WARNINGF(TL("A TAZ with id '%' already exists. Not building junction TAZ."), tazID);
273 8 : continue;
274 : }
275 1435 : const std::string sourceID = tazID + "-source";
276 1435 : const std::string sinkID = tazID + "-sink";
277 : // sink must be added before source
278 2870 : ROEdge* sink = eb.buildEdge(sinkID, nullptr, nullptr, 0, "", "");
279 2870 : ROEdge* source = eb.buildEdge(sourceID, nullptr, nullptr, 0, "", "");
280 : sink->setOtherTazConnector(source);
281 : source->setOtherTazConnector(sink);
282 2870 : if (!addDistrict(tazID, source, sink)) {
283 : continue;
284 : }
285 1435 : auto& district = myDistricts[tazID];
286 1435 : const RONode* junction = item.second;
287 12452 : for (const ROEdge* edge : junction->getIncoming()) {
288 11017 : if (!edge->isInternal()) {
289 3129 : const_cast<ROEdge*>(edge)->addSuccessor(sink);
290 3129 : district.second.push_back(edge->getID());
291 : }
292 : }
293 12452 : for (const ROEdge* edge : junction->getOutgoing()) {
294 11017 : if (!edge->isInternal()) {
295 3129 : source->addSuccessor(const_cast<ROEdge*>(edge));
296 3129 : district.first.push_back(edge->getID());
297 : }
298 : }
299 : }
300 101 : }
301 :
302 :
303 : void
304 3637 : RONet::setBidiEdges(const std::map<ROEdge*, std::string>& bidiMap) {
305 16189 : for (const auto& item : bidiMap) {
306 12552 : ROEdge* bidi = myEdges.get(item.second);
307 12552 : if (bidi == nullptr) {
308 0 : WRITE_ERRORF(TL("The bidi edge '%' is not known."), item.second);
309 : }
310 12552 : item.first->setBidiEdge(bidi);
311 12552 : myHasBidiEdges = true;
312 : }
313 3637 : }
314 :
315 :
316 : void
317 74576 : RONet::addNode(RONode* node) {
318 74576 : if (!myNodes.add(node->getID(), node)) {
319 0 : WRITE_ERRORF(TL("The node '%' occurs at least twice."), node->getID());
320 0 : delete node;
321 : }
322 74576 : }
323 :
324 :
325 : void
326 1980 : RONet::addStoppingPlace(const std::string& id, const SumoXMLTag category, SUMOVehicleParameter::Stop* stop) {
327 3580 : if (!myStoppingPlaces[category == SUMO_TAG_TRAIN_STOP ? SUMO_TAG_BUS_STOP : category].add(id, stop)) {
328 12 : WRITE_ERRORF(TL("The % '%' occurs at least twice."), toString(category), id);
329 4 : delete stop;
330 : }
331 1980 : }
332 :
333 :
334 : bool
335 281276 : RONet::addRouteDef(RORouteDef* def) {
336 281276 : return myRoutes.add(def->getID(), def);
337 : }
338 :
339 :
340 : void
341 2934 : RONet::openOutput(const OptionsCont& options) {
342 8800 : if (options.isSet("output-file") && options.getString("output-file") != "") {
343 2933 : myRoutesOutput = &OutputDevice::getDevice(options.getString("output-file"));
344 2927 : if (myRoutesOutput->isNull()) {
345 1 : myRoutesOutput = nullptr;
346 : } else {
347 5852 : myRoutesOutput->writeXMLHeader("routes", "routes_file.xsd");
348 : }
349 : }
350 8598 : if (options.exists("alternatives-output") && options.isSet("alternatives-output")
351 13699 : && !(options.exists("write-trips") && options.getBool("write-trips"))) {
352 2547 : myRouteAlternativesOutput = &OutputDevice::getDevice(options.getString("alternatives-output"));
353 2547 : if (myRouteAlternativesOutput->isNull()) {
354 165 : myRouteAlternativesOutput = nullptr;
355 : } else {
356 4764 : myRouteAlternativesOutput->writeXMLHeader("routes", "routes_file.xsd");
357 : }
358 : }
359 5856 : if (options.isSet("vtype-output")) {
360 32 : myTypesOutput = &OutputDevice::getDevice(options.getString("vtype-output"));
361 64 : myTypesOutput->writeXMLHeader("routes", "routes_file.xsd");
362 : }
363 2928 : }
364 :
365 :
366 : void
367 2710 : RONet::writeIntermodal(const OptionsCont& options, ROIntermodalRouter& router) const {
368 5420 : if (options.exists("intermodal-network-output") && options.isSet("intermodal-network-output")) {
369 32 : OutputDevice::createDeviceByOption("intermodal-network-output", "intermodal");
370 32 : router.writeNetwork(OutputDevice::getDevice(options.getString("intermodal-network-output")));
371 : }
372 5420 : if (options.exists("intermodal-weight-output") && options.isSet("intermodal-weight-output")) {
373 16 : OutputDevice::createDeviceByOption("intermodal-weight-output", "weights", "meandata_file.xsd");
374 8 : OutputDevice& dev = OutputDevice::getDeviceByOption("intermodal-weight-output");
375 8 : dev.openTag(SUMO_TAG_INTERVAL);
376 8 : dev.writeAttr(SUMO_ATTR_ID, "intermodalweights");
377 8 : dev.writeAttr(SUMO_ATTR_BEGIN, 0);
378 8 : dev.writeAttr(SUMO_ATTR_END, SUMOTime_MAX);
379 8 : router.writeWeights(dev);
380 16 : dev.closeTag();
381 : }
382 2710 : }
383 :
384 :
385 : void
386 2933 : RONet::cleanup() {
387 : // end writing
388 2933 : if (myRoutesOutput != nullptr) {
389 2926 : myRoutesOutput->close();
390 : }
391 : // only if opened
392 2933 : if (myRouteAlternativesOutput != nullptr) {
393 2382 : myRouteAlternativesOutput->close();
394 : }
395 : // only if opened
396 2933 : if (myTypesOutput != nullptr) {
397 32 : myTypesOutput->close();
398 : }
399 : RouteCostCalculator<RORoute, ROEdge, ROVehicle>::cleanup();
400 : #ifdef HAVE_FOX
401 2933 : if (myThreadPool.size() > 0) {
402 236 : myThreadPool.clear();
403 : }
404 : #endif
405 2933 : }
406 :
407 :
408 :
409 : SUMOVTypeParameter*
410 498035 : RONet::getVehicleTypeSecure(const std::string& id) {
411 : // check whether the type was already known
412 : SUMOVTypeParameter* type = myVehicleTypes.get(id);
413 498035 : if (id == DEFAULT_VTYPE_ID) {
414 455574 : myDefaultVTypeMayBeDeleted = false;
415 42461 : } else if (id == DEFAULT_PEDTYPE_ID) {
416 5748 : myDefaultPedTypeMayBeDeleted = false;
417 36713 : } else if (id == DEFAULT_BIKETYPE_ID) {
418 54 : myDefaultBikeTypeMayBeDeleted = false;
419 36659 : } else if (id == DEFAULT_TAXITYPE_ID) {
420 93 : myDefaultTaxiTypeMayBeDeleted = false;
421 36566 : } else if (id == DEFAULT_RAILTYPE_ID) {
422 0 : myDefaultRailTypeMayBeDeleted = false;
423 : }
424 498035 : if (type != nullptr) {
425 : return type;
426 : }
427 : VTypeDistDictType::iterator it2 = myVTypeDistDict.find(id);
428 2462 : if (it2 != myVTypeDistDict.end()) {
429 2132 : return it2->second->get();
430 : }
431 330 : if (id == "") {
432 : // ok, no vehicle type or an unknown type was given within the user input
433 : // return the default type
434 0 : myDefaultVTypeMayBeDeleted = false;
435 : return myVehicleTypes.get(DEFAULT_VTYPE_ID);
436 : }
437 : return type;
438 : }
439 :
440 :
441 : bool
442 2234 : RONet::checkVType(const std::string& id) {
443 2234 : if (id == DEFAULT_VTYPE_ID) {
444 68 : if (myDefaultVTypeMayBeDeleted) {
445 68 : myVehicleTypes.remove(id);
446 68 : myDefaultVTypeMayBeDeleted = false;
447 : } else {
448 : return false;
449 : }
450 2166 : } else if (id == DEFAULT_PEDTYPE_ID) {
451 0 : if (myDefaultPedTypeMayBeDeleted) {
452 0 : myVehicleTypes.remove(id);
453 0 : myDefaultPedTypeMayBeDeleted = false;
454 : } else {
455 : return false;
456 : }
457 2166 : } else if (id == DEFAULT_BIKETYPE_ID) {
458 0 : if (myDefaultBikeTypeMayBeDeleted) {
459 0 : myVehicleTypes.remove(id);
460 0 : myDefaultBikeTypeMayBeDeleted = false;
461 : } else {
462 : return false;
463 : }
464 2166 : } else if (id == DEFAULT_TAXITYPE_ID) {
465 0 : if (myDefaultTaxiTypeMayBeDeleted) {
466 0 : myVehicleTypes.remove(id);
467 0 : myDefaultTaxiTypeMayBeDeleted = false;
468 : } else {
469 : return false;
470 : }
471 2166 : } else if (id == DEFAULT_RAILTYPE_ID) {
472 0 : if (myDefaultRailTypeMayBeDeleted) {
473 0 : myVehicleTypes.remove(id);
474 0 : myDefaultRailTypeMayBeDeleted = false;
475 : } else {
476 : return false;
477 : }
478 : } else {
479 2166 : if (myVehicleTypes.get(id) != 0 || myVTypeDistDict.find(id) != myVTypeDistDict.end()) {
480 0 : return false;
481 : }
482 : }
483 : return true;
484 : }
485 :
486 :
487 : bool
488 2150 : RONet::addVehicleType(SUMOVTypeParameter* type) {
489 2150 : if (checkVType(type->id)) {
490 2150 : myVehicleTypes.add(type->id, type);
491 : } else {
492 0 : WRITE_ERRORF(TL("The vehicle type '%' occurs at least twice."), type->id);
493 0 : delete type;
494 0 : return false;
495 : }
496 2150 : return true;
497 : }
498 :
499 :
500 : bool
501 84 : RONet::addVTypeDistribution(const std::string& id, RandomDistributor<SUMOVTypeParameter*>* vehTypeDistribution) {
502 84 : if (checkVType(id)) {
503 84 : myVTypeDistDict[id] = vehTypeDistribution;
504 84 : return true;
505 : }
506 0 : delete vehTypeDistribution;
507 : return false;
508 : }
509 :
510 :
511 : bool
512 298171 : RONet::addVehicle(const std::string& id, ROVehicle* veh) {
513 298171 : if (myVehIDs.find(id) == myVehIDs.end()) {
514 596302 : myVehIDs[id] = veh->getParameter().departProcedure == DepartDefinition::TRIGGERED ? -1 : veh->getDepartureTime();
515 :
516 298167 : if (veh->isPublicTransport()) {
517 44 : if (!veh->isPartOfFlow()) {
518 36 : myPTVehicles.push_back(veh);
519 : }
520 44 : if (!myDoPTRouting) {
521 : return true;
522 : }
523 : }
524 298131 : myRoutables[veh->getDepart()].push_back(veh);
525 298131 : return true;
526 : }
527 12 : WRITE_ERRORF(TL("Another vehicle with the id '%' exists."), id);
528 4 : delete veh;
529 : return false;
530 : }
531 :
532 :
533 : bool
534 256 : RONet::knowsVehicle(const std::string& id) const {
535 256 : return myVehIDs.find(id) != myVehIDs.end();
536 : }
537 :
538 :
539 : bool
540 16 : RONet::knowsVType(const std::string& id) const {
541 16 : return myVehicleTypes.get(id) != nullptr || myVTypeDistDict.find(id) != myVTypeDistDict.end();
542 : }
543 :
544 :
545 : SUMOTime
546 28 : RONet::getDeparture(const std::string& vehID) const {
547 : auto it = myVehIDs.find(vehID);
548 28 : if (it != myVehIDs.end()) {
549 28 : return it->second;
550 : } else {
551 0 : throw ProcessError(TLF("Requesting departure time for unknown vehicle '%'", vehID));
552 : }
553 : }
554 :
555 :
556 : bool
557 2586 : RONet::addFlow(SUMOVehicleParameter* flow, const bool randomize) {
558 2586 : if (randomize && flow->repetitionOffset >= 0) {
559 10 : myDepartures[flow->id].reserve(flow->repetitionNumber);
560 110 : for (int i = 0; i < flow->repetitionNumber; ++i) {
561 100 : myDepartures[flow->id].push_back(flow->depart + RandHelper::rand(flow->repetitionNumber * flow->repetitionOffset));
562 : }
563 10 : std::sort(myDepartures[flow->id].begin(), myDepartures[flow->id].end());
564 10 : std::reverse(myDepartures[flow->id].begin(), myDepartures[flow->id].end());
565 2576 : } else if (flow->repetitionOffset < 0) {
566 : // init poisson flow (but only the timing)
567 349 : flow->incrementFlow(1);
568 349 : flow->repetitionsDone--;
569 : }
570 2586 : const bool added = myFlows.add(flow->id, flow);
571 2586 : if (added) {
572 2582 : myHaveActiveFlows = true;
573 : }
574 2586 : return added;
575 : }
576 :
577 :
578 : bool
579 3551 : RONet::addPerson(ROPerson* person) {
580 : if (myPersonIDs.count(person->getID()) == 0) {
581 : myPersonIDs.insert(person->getID());
582 3551 : myRoutables[person->getDepart()].push_back(person);
583 3551 : return true;
584 : }
585 0 : WRITE_ERRORF(TL("Another person with the id '%' exists."), person->getID());
586 0 : return false;
587 : }
588 :
589 :
590 : void
591 56 : RONet::addContainer(const SUMOTime depart, const std::string desc) {
592 56 : myContainers.insert(std::pair<const SUMOTime, const std::string>(depart, desc));
593 56 : }
594 :
595 :
596 : void
597 13382 : RONet::checkFlows(SUMOTime time, MsgHandler* errorHandler) {
598 13382 : myHaveActiveFlows = false;
599 51426 : for (const auto& i : myFlows) {
600 38044 : SUMOVehicleParameter* const pars = i.second;
601 38044 : if (pars->line != "" && !myDoPTRouting) {
602 745 : continue;
603 : }
604 37299 : if (myKeepFlows) {
605 796 : if (pars->repetitionsDone < pars->repetitionNumber) {
606 : // each each flow only once
607 790 : pars->repetitionsDone = pars->repetitionNumber;
608 790 : const SUMOVTypeParameter* type = getVehicleTypeSecure(pars->vtypeid);
609 790 : if (type == nullptr) {
610 0 : type = getVehicleTypeSecure(DEFAULT_VTYPE_ID);
611 : } else {
612 : auto dist = getVTypeDistribution(pars->vtypeid);
613 0 : if (dist != nullptr) {
614 0 : WRITE_WARNINGF("Keeping flow '%' with a vTypeDistribution can lead to invalid routes if the distribution contains different vClasses", pars->id);
615 : }
616 : }
617 1580 : RORouteDef* route = getRouteDef(pars->routeid)->copy(pars->routeid, pars->depart);
618 790 : ROVehicle* veh = new ROVehicle(*pars, route, type, this, errorHandler);
619 790 : addVehicle(pars->id, veh);
620 : }
621 36503 : } else if (pars->repetitionProbability > 0) {
622 108 : if (pars->repetitionEnd > pars->depart && pars->repetitionsDone < pars->repetitionNumber) {
623 96 : myHaveActiveFlows = true;
624 : }
625 : const SUMOTime origDepart = pars->depart;
626 1020 : while (pars->depart < time && pars->repetitionsDone < pars->repetitionNumber) {
627 928 : if (pars->repetitionEnd <= pars->depart) {
628 : break;
629 : }
630 : // only call rand if all other conditions are met
631 912 : if (RandHelper::rand() < (pars->repetitionProbability * TS)) {
632 168 : SUMOVehicleParameter* newPars = new SUMOVehicleParameter(*pars);
633 336 : newPars->id = pars->id + "." + toString(pars->repetitionsDone);
634 168 : newPars->depart = pars->depart;
635 168 : for (StopParVector::iterator stop = newPars->stops.begin(); stop != newPars->stops.end(); ++stop) {
636 0 : if (stop->until >= 0) {
637 0 : stop->until += pars->depart - origDepart;
638 : }
639 : }
640 168 : pars->repetitionsDone++;
641 : // try to build the vehicle
642 168 : const SUMOVTypeParameter* type = getVehicleTypeSecure(pars->vtypeid);
643 168 : if (type == nullptr) {
644 32 : type = getVehicleTypeSecure(DEFAULT_VTYPE_ID);
645 136 : } else if (!myKeepVTypeDist) {
646 : // fix the type id in case we used a distribution
647 136 : newPars->vtypeid = type->id;
648 : }
649 168 : const SUMOTime stopOffset = pars->routeid[0] == '!' ? pars->depart - origDepart : pars->depart;
650 336 : RORouteDef* route = getRouteDef(pars->routeid)->copy("!" + newPars->id, stopOffset);
651 168 : ROVehicle* veh = new ROVehicle(*newPars, route, type, this, errorHandler);
652 168 : addVehicle(newPars->id, veh);
653 168 : delete newPars;
654 : }
655 912 : pars->depart += DELTA_T;
656 : }
657 : } else {
658 36395 : SUMOTime depart = static_cast<SUMOTime>(pars->depart + pars->repetitionTotalOffset);
659 55118 : while (pars->repetitionsDone < pars->repetitionNumber && pars->repetitionEnd >= depart) {
660 29477 : myHaveActiveFlows = true;
661 29477 : depart = static_cast<SUMOTime>(pars->depart + pars->repetitionTotalOffset);
662 29477 : if (myDepartures.find(pars->id) != myDepartures.end()) {
663 110 : depart = myDepartures[pars->id].back();
664 : }
665 29477 : if (depart >= time + DELTA_T) {
666 : break;
667 : }
668 18723 : if (myDepartures.find(pars->id) != myDepartures.end()) {
669 100 : myDepartures[pars->id].pop_back();
670 : }
671 18723 : SUMOVehicleParameter* newPars = new SUMOVehicleParameter(*pars);
672 37446 : newPars->id = pars->id + "." + toString(pars->repetitionsDone);
673 18723 : newPars->depart = depart;
674 18883 : for (StopParVector::iterator stop = newPars->stops.begin(); stop != newPars->stops.end(); ++stop) {
675 160 : if (stop->until >= 0) {
676 40 : stop->until += depart - pars->depart;
677 : }
678 160 : if (stop->arrival >= 0) {
679 40 : stop->arrival += depart - pars->depart;
680 : }
681 : }
682 18723 : pars->incrementFlow(1);
683 : // try to build the vehicle
684 18723 : const SUMOVTypeParameter* type = getVehicleTypeSecure(pars->vtypeid);
685 18723 : if (type == nullptr) {
686 62 : type = getVehicleTypeSecure(DEFAULT_VTYPE_ID);
687 18661 : } else if (!myKeepVTypeDist) {
688 : // fix the type id in case we used a distribution
689 18621 : newPars->vtypeid = type->id;
690 : }
691 18723 : const SUMOTime stopOffset = pars->routeid[0] == '!' ? depart - pars->depart : depart;
692 37446 : RORouteDef* route = getRouteDef(pars->routeid)->copy("!" + newPars->id, stopOffset);
693 18723 : ROVehicle* veh = new ROVehicle(*newPars, route, type, this, errorHandler);
694 18723 : addVehicle(newPars->id, veh);
695 18723 : delete newPars;
696 : }
697 : }
698 : }
699 13382 : }
700 :
701 :
702 : void
703 74 : RONet::createBulkRouteRequests(const RORouterProvider& provider, const SUMOTime time, const bool removeLoops) {
704 : std::map<const int, std::vector<RORoutable*> > bulkVehs;
705 : int numBulked = 0;
706 111 : for (RoutablesMap::const_iterator i = myRoutables.begin(); i != myRoutables.end(); ++i) {
707 74 : if (i->first >= time) {
708 : break;
709 : }
710 694 : for (RORoutable* const routable : i->second) {
711 657 : const ROEdge* const depEdge = routable->getDepartEdge();
712 657 : bulkVehs[depEdge->getNumericalID()].push_back(routable);
713 657 : RORoutable* const first = bulkVehs[depEdge->getNumericalID()].front();
714 657 : numBulked++;
715 657 : if (first->getMaxSpeed() != routable->getMaxSpeed()) {
716 6 : WRITE_WARNINGF(TL("Bulking different maximum speeds ('%' and '%') may lead to suboptimal routes."), first->getID(), routable->getID());
717 : }
718 657 : if (first->getVClass() != routable->getVClass()) {
719 6 : WRITE_WARNINGF(TL("Bulking different vehicle classes ('%' and '%') may lead to invalid routes."), first->getID(), routable->getID());
720 : }
721 : }
722 : }
723 : #ifdef HAVE_FOX
724 : int workerIndex = 0;
725 : #endif
726 74 : if ((int)bulkVehs.size() < numBulked) {
727 72 : WRITE_MESSAGE(TLF("Using bulk-mode for % entities from % origins", numBulked, bulkVehs.size()));
728 : }
729 235 : for (std::map<const int, std::vector<RORoutable*> >::const_iterator i = bulkVehs.begin(); i != bulkVehs.end(); ++i) {
730 : #ifdef HAVE_FOX
731 161 : if (myThreadPool.size() > 0) {
732 : bool bulk = true;
733 132 : for (RORoutable* const r : i->second) {
734 120 : myThreadPool.add(new RoutingTask(r, removeLoops, myErrorHandler), workerIndex);
735 120 : if (bulk) {
736 12 : myThreadPool.add(new BulkmodeTask(true), workerIndex);
737 : bulk = false;
738 : }
739 : }
740 12 : myThreadPool.add(new BulkmodeTask(false), workerIndex);
741 12 : workerIndex++;
742 12 : if (workerIndex == (int)myThreadPool.size()) {
743 : workerIndex = 0;
744 : }
745 12 : continue;
746 12 : }
747 : #endif
748 686 : for (RORoutable* const r : i->second) {
749 537 : r->computeRoute(provider, removeLoops, myErrorHandler);
750 537 : provider.setBulkMode(true);
751 : }
752 149 : provider.setBulkMode(false);
753 : }
754 74 : }
755 :
756 :
757 : SUMOTime
758 19548 : RONet::saveAndRemoveRoutesUntil(OptionsCont& options, const RORouterProvider& provider,
759 : SUMOTime time) {
760 19548 : MsgHandler* mh = (options.getBool("ignore-errors") ?
761 19548 : MsgHandler::getWarningInstance() : MsgHandler::getErrorInstance());
762 19548 : if (myHaveActiveFlows) {
763 13382 : checkFlows(time, mh);
764 : }
765 : SUMOTime lastTime = -1;
766 19548 : const bool removeLoops = options.getBool("remove-loops");
767 : #ifdef HAVE_FOX
768 39096 : const int maxNumThreads = options.getInt("routing-threads");
769 : #endif
770 19548 : if (myRoutables.size() != 0) {
771 19040 : if (options.getBool("bulk-routing")) {
772 : #ifdef HAVE_FOX
773 82 : while ((int)myThreadPool.size() < maxNumThreads) {
774 8 : new WorkerThread(myThreadPool, provider);
775 : }
776 : #endif
777 74 : createBulkRouteRequests(provider, time, removeLoops);
778 : } else {
779 252489 : for (RoutablesMap::const_iterator i = myRoutables.begin(); i != myRoutables.end(); ++i) {
780 249645 : if (i->first >= time) {
781 : break;
782 : }
783 543871 : for (RORoutable* const routable : i->second) {
784 : #ifdef HAVE_FOX
785 : // add task
786 300828 : if (maxNumThreads > 0) {
787 : const int numThreads = (int)myThreadPool.size();
788 179640 : if (numThreads == 0) {
789 : // This is the very first routing. Since at least the CHRouter needs initialization
790 : // before it gets cloned, we do not do this in parallel
791 228 : routable->computeRoute(provider, removeLoops, myErrorHandler);
792 228 : new WorkerThread(myThreadPool, provider);
793 : } else {
794 : // add thread if necessary
795 179412 : if (numThreads < maxNumThreads && myThreadPool.isFull()) {
796 1153 : new WorkerThread(myThreadPool, provider);
797 : }
798 179412 : myThreadPool.add(new RoutingTask(routable, removeLoops, myErrorHandler));
799 : }
800 179640 : continue;
801 179640 : }
802 : #endif
803 121188 : routable->computeRoute(provider, removeLoops, myErrorHandler);
804 : }
805 : }
806 : }
807 : #ifdef HAVE_FOX
808 9504 : myThreadPool.waitAll();
809 : #endif
810 : }
811 38792 : const double scale = options.exists("scale-suffix") ? options.getFloat("scale") : 1;
812 : // write all vehicles (and additional structures)
813 262660 : while (myRoutables.size() != 0 || myContainers.size() != 0) {
814 : // get the next vehicle, person or container
815 : RoutablesMap::iterator routables = myRoutables.begin();
816 249759 : const SUMOTime routableTime = routables == myRoutables.end() ? SUMOTime_MAX : routables->first;
817 : ContainerMap::iterator container = myContainers.begin();
818 249759 : const SUMOTime containerTime = container == myContainers.end() ? SUMOTime_MAX : container->first;
819 : // check whether it shall not yet be computed
820 249759 : if (routableTime >= time && containerTime >= time) {
821 : lastTime = MIN2(routableTime, containerTime);
822 : break;
823 : }
824 : const SUMOTime minTime = MIN2(routableTime, containerTime);
825 243128 : if (routableTime == minTime) {
826 : // check whether to print the output
827 243080 : if (lastTime != routableTime && lastTime != -1) {
828 : // report writing progress
829 236764 : if (options.getInt("stats-period") >= 0 && ((int)routableTime % options.getInt("stats-period")) == 0) {
830 0 : WRITE_MESSAGE("Read: " + toString(myVehIDs.size()) + ", Discarded: " + toString(myDiscardedRouteNo) + ", Written: " + toString(myWrittenRouteNo));
831 : }
832 : }
833 : lastTime = routableTime;
834 544549 : for (const RORoutable* const r : routables->second) {
835 : // ok, check whether it has been routed
836 301469 : if (r->getRoutingSuccess()) {
837 : // write the route
838 282172 : int quota = getScalingQuota(scale, myWrittenRouteNo);
839 282172 : r->write(myRoutesOutput, myRouteAlternativesOutput, myTypesOutput, options, quota);
840 282172 : myWrittenRouteNo++;
841 : } else {
842 19297 : myDiscardedRouteNo++;
843 : }
844 : // we need to keep individual public transport vehicles but not the flows
845 301469 : if (!r->isPublicTransport() || r->isPartOfFlow()) {
846 : // delete routes and the vehicle
847 301469 : const ROVehicle* const veh = dynamic_cast<const ROVehicle*>(r);
848 301469 : if (veh != nullptr && veh->getRouteDefinition()->getID()[0] == '!') {
849 297927 : if (r->isPartOfFlow() || !myRoutes.remove(veh->getRouteDefinition()->getID())) {
850 20880 : delete veh->getRouteDefinition();
851 : }
852 : }
853 301469 : delete r;
854 : }
855 : }
856 : myRoutables.erase(routables);
857 : }
858 243128 : if (containerTime == minTime) {
859 48 : myRoutesOutput->writePreformattedTag(container->second);
860 48 : if (myRouteAlternativesOutput != nullptr) {
861 : myRouteAlternativesOutput->writePreformattedTag(container->second);
862 : }
863 : myContainers.erase(container);
864 : }
865 : }
866 19532 : return lastTime;
867 : }
868 :
869 :
870 : bool
871 41933 : RONet::furtherStored() {
872 41933 : return myRoutables.size() > 0 || (myFlows.size() > 0 && myHaveActiveFlows) || myContainers.size() > 0;
873 : }
874 :
875 :
876 : int
877 123 : RONet::getEdgeNumber() const {
878 123 : return myEdges.size();
879 : }
880 :
881 :
882 : int
883 1 : RONet::getInternalEdgeNumber() const {
884 1 : return myNumInternalEdges;
885 : }
886 :
887 :
888 : ROEdge*
889 9059 : RONet::getEdgeForLaneID(const std::string& laneID) const {
890 18118 : return getEdge(SUMOXMLDefinitions::getEdgeIDFromLane(laneID));
891 : }
892 :
893 :
894 : ROLane*
895 1140 : RONet::getLane(const std::string& laneID) const {
896 1140 : int laneIndex = SUMOXMLDefinitions::getIndexFromLane(laneID);
897 1140 : return getEdgeForLaneID(laneID)->getLanes()[laneIndex];
898 : }
899 :
900 :
901 : void
902 757 : RONet::adaptIntermodalRouter(ROIntermodalRouter& router) {
903 757 : double taxiWait = STEPS2TIME(string2time(OptionsCont::getOptions().getString("persontrip.taxi.waiting-time")));
904 1048 : for (const auto& stopType : myInstance->myStoppingPlaces) {
905 : // add access to all stopping places
906 291 : const SumoXMLTag element = stopType.first;
907 2078 : for (const auto& stop : stopType.second) {
908 1787 : router.getNetwork()->addAccess(stop.first, myInstance->getEdgeForLaneID(stop.second->lane),
909 1787 : stop.second->startPos, stop.second->endPos, 0., element, false, taxiWait);
910 : // add access to all public transport stops
911 1787 : if (element == SUMO_TAG_BUS_STOP) {
912 3843 : for (const auto& a : stop.second->accessPos) {
913 2076 : router.getNetwork()->addAccess(stop.first, myInstance->getEdgeForLaneID(std::get<0>(a)),
914 : std::get<1>(a), std::get<1>(a), std::get<2>(a), SUMO_TAG_BUS_STOP, true, taxiWait);
915 : }
916 : }
917 : }
918 : }
919 : // fill the public transport router with pre-parsed public transport lines
920 2220 : for (const auto& i : myInstance->myFlows) {
921 1463 : if (i.second->line != "") {
922 749 : const RORouteDef* const route = myInstance->getRouteDef(i.second->routeid);
923 : const StopParVector* addStops = nullptr;
924 2247 : if (route != nullptr && route->getFirstRoute() != nullptr) {
925 749 : addStops = &route->getFirstRoute()->getStops();
926 : }
927 749 : router.getNetwork()->addSchedule(*i.second, addStops);
928 : }
929 : }
930 793 : for (const RORoutable* const veh : myInstance->myPTVehicles) {
931 : // add single vehicles with line attribute which are not part of a flow
932 : // no need to add route stops here, they have been added to the vehicle before
933 36 : router.getNetwork()->addSchedule(veh->getParameter());
934 : }
935 : // add access to transfer from walking to taxi-use
936 757 : if ((router.getCarWalkTransfer() & ModeChangeOptions::TAXI_PICKUP_ANYWHERE) != 0) {
937 224775 : for (const ROEdge* edge : ROEdge::getAllEdges()) {
938 224030 : if (!edge->isTazConnector() && (edge->getPermissions() & SVC_PEDESTRIAN) != 0 && (edge->getPermissions() & SVC_TAXI) != 0) {
939 36652 : router.getNetwork()->addCarAccess(edge, SVC_TAXI, taxiWait);
940 : }
941 : }
942 : }
943 757 : }
944 :
945 :
946 : bool
947 4172852 : RONet::hasPermissions() const {
948 4172852 : return myHavePermissions;
949 : }
950 :
951 :
952 : void
953 254700 : RONet::setPermissionsFound() {
954 254700 : myHavePermissions = true;
955 254700 : }
956 :
957 : bool
958 12 : RONet::hasLoadedEffort() const {
959 38 : for (const auto& item : myEdges) {
960 35 : if (item.second->hasStoredEffort()) {
961 : return true;
962 : }
963 : }
964 : return false;
965 : }
966 :
967 : const std::string
968 1629 : RONet::getStoppingPlaceName(const std::string& id) const {
969 1821 : for (const auto& mapItem : myStoppingPlaces) {
970 : SUMOVehicleParameter::Stop* stop = mapItem.second.get(id);
971 1629 : if (stop != nullptr) {
972 : // see RONetHandler::parseStoppingPlace
973 : return stop->busstop;
974 : }
975 : }
976 0 : return "";
977 : }
978 :
979 : const std::string
980 1509 : RONet::getStoppingPlaceElement(const std::string& id) const {
981 1653 : for (const auto& mapItem : myStoppingPlaces) {
982 : SUMOVehicleParameter::Stop* stop = mapItem.second.get(id);
983 1509 : if (stop != nullptr) {
984 : // see RONetHandler::parseStoppingPlace
985 : return stop->actType;
986 : }
987 : }
988 0 : return toString(SUMO_TAG_BUS_STOP);
989 : }
990 :
991 :
992 : void
993 4 : RONet::addProhibition(const ROEdge* edge, const RouterProhibition& prohibition) {
994 : if (myProhibitions.count(edge) != 0) {
995 0 : throw ProcessError(TLF("Already loaded prohibition for edge '%'. (Only one prohibition per edge is supported)", edge->getID()));
996 : }
997 4 : myProhibitions[edge] = prohibition;
998 4 : myHavePermissions = true;
999 4 : }
1000 :
1001 :
1002 : void
1003 4 : RONet::addLaneProhibition(const ROLane* lane, const RouterProhibition& prohibition) {
1004 : if (myLaneProhibitions.count(lane) != 0) {
1005 0 : throw ProcessError(TLF("Already loaded prohibition for lane '%'. (Only one prohibition per lane is supported)", lane->getID()));
1006 : }
1007 : assert(prohibition.end > prohibition.begin);
1008 4 : myLaneProhibitions[lane] = prohibition;
1009 4 : myLaneProhibitionTimes[prohibition.begin].insert(lane);
1010 4 : myHavePermissions = true;
1011 4 : }
1012 :
1013 :
1014 : void
1015 20 : RONet::updateLaneProhibitions(SUMOTime begin) {
1016 20 : const double beginS = STEPS2TIME(begin);
1017 24 : while (myLaneProhibitionTimes.size() > 0 && myLaneProhibitionTimes.begin()->first <= beginS) {
1018 : const double t = myLaneProhibitionTimes.begin()->first;
1019 8 : for (const ROLane* const lane : myLaneProhibitionTimes.begin()->second) {
1020 : const SVCPermissions orig = lane->getPermissions();
1021 : assert(myLaneProhibitions.count(lane) != 0);
1022 4 : RouterProhibition& rp = myLaneProhibitions[lane];
1023 4 : const_cast<ROLane*>(lane)->setPermissions(rp.permissions);
1024 4 : lane->getEdge().resetSuccessors();
1025 4 : for (ROEdge* pred : lane->getEdge().getPredecessors()) {
1026 0 : pred->resetSuccessors();
1027 : }
1028 4 : if (t == rp.begin) {
1029 : // schedule restoration of original permissions. This works
1030 : // without a stack because there is at most one prohibition per lane
1031 2 : myLaneProhibitionTimes[rp.end].insert(lane);
1032 2 : rp.permissions = orig;
1033 : }
1034 : }
1035 : myLaneProhibitionTimes.erase(myLaneProhibitionTimes.begin());
1036 : }
1037 20 : }
1038 :
1039 :
1040 : #ifdef HAVE_FOX
1041 : // ---------------------------------------------------------------------------
1042 : // RONet::RoutingTask-methods
1043 : // ---------------------------------------------------------------------------
1044 : void
1045 179532 : RONet::RoutingTask::run(MFXWorkerThread* context) {
1046 179532 : myRoutable->computeRoute(*static_cast<WorkerThread*>(context), myRemoveLoops, myErrorHandler);
1047 179532 : }
1048 : #endif
1049 :
1050 :
1051 : /****************************************************************************/
|
