52 : myDistricts(dc), myNumLoaded(0), myNumWritten(0), myNumDiscarded(0), myBegin(-1), myEnd(-1) {}
57 for (
RORoute*
const r : cell->pathsVector) {
68 SUMOTime end,
const std::string& origin,
const std::string& destination,
69 const std::string& vehicleType,
const bool originIsEdge,
const bool destinationIsEdge) {
72 WRITE_WARNING(
"Missing origin '" + origin +
"' and destination '" + destination +
"' (" +
toString(vehicleNumber) +
" vehicles).");
77 }
else if (!originIsEdge &&
myDistricts.
get(origin) == 0 && vehicleNumber > 0) {
82 }
else if (!destinationIsEdge &&
myDistricts.
get(destination) == 0 && vehicleNumber > 0) {
83 WRITE_ERROR(
"Missing destination '" + destination +
"' (" +
toString(vehicleNumber) +
" vehicles).");
89 WRITE_ERROR(
"District '" + origin +
"' has no source.");
93 WRITE_ERROR(
"District '" + destination +
"' has no sink.");
113 const std::string& fromTaz,
const std::string& toTaz,
114 const std::string& vehicleType,
const bool originIsEdge,
const bool destinationIsEdge) {
121 std::vector<ODCell*>& odList =
myShortCut[std::make_pair(fromTaz, toTaz)];
123 for (std::vector<ODCell*>::const_reverse_iterator c = odList.rbegin(); c != odList.rend(); ++c) {
124 if ((*c)->begin <= depart && (*c)->end > depart && (*c)->vehicleType == vehicleType) {
129 if (cell ==
nullptr) {
131 const int intervalIdx = (int)(depart / interval);
132 if (
add(1., intervalIdx * interval, (intervalIdx + 1) * interval, fromTaz, toTaz, vehicleType, originIsEdge, destinationIsEdge)) {
134 odList.push_back(cell);
149 int& vehName, std::vector<ODVehicle>& into,
150 const bool uniform,
const bool differSourceSink,
151 const std::string& prefix) {
157 if (vehicles2insert == 0) {
161 const double offset = (double)(cell->
end - cell->
begin) / (double) vehicles2insert / (
double) 2.;
162 for (
int i = 0; i < vehicles2insert; ++i) {
175 }
while (canDiffer && differSourceSink && (veh.
to == veh.
from));
176 if (!canDiffer && differSourceSink && (veh.
to == veh.
from)) {
188 const ODCell*
const cell) {
194 if (oc.
isSet(
"departlane") && oc.
getString(
"departlane") !=
"default") {
197 if (oc.
isSet(
"departpos")) {
200 if (oc.
isSet(
"departspeed") && oc.
getString(
"departspeed") !=
"default") {
203 if (oc.
isSet(
"arrivallane")) {
206 if (oc.
isSet(
"arrivalpos")) {
209 if (oc.
isSet(
"arrivalspeed")) {
218 const bool differSourceSink,
const bool noVtype,
219 const std::string& prefix,
const bool stepLog,
220 bool pedestrians,
bool persontrips,
221 const std::string& modes) {
225 std::map<std::pair<std::string, std::string>,
double> fractionLeft;
230 std::vector<ODCell*>::iterator next =
myContainer.begin();
231 std::vector<ODVehicle> vehicles;
234 for (
SUMOTime t = begin; t < end;) {
235 if (stepLog && t - lastOut >=
DELTA_T) {
236 std::cout <<
"Parsing time " +
time2string(t) <<
'\r';
240 bool changed =
false;
241 while (next !=
myContainer.end() && (*next)->begin <= t && (*next)->end > t) {
242 std::pair<std::string, std::string> odID = std::make_pair((*next)->origin, (*next)->destination);
244 if (fractionLeft.find(odID) != fractionLeft.end()) {
245 (*next)->vehicleNumber += fractionLeft[odID];
246 fractionLeft[odID] = 0;
249 const int oldSize = (int)vehicles.size();
250 const double fraction =
computeDeparts(*next, vehName, vehicles, uniform, differSourceSink, prefix);
251 if (oldSize != (
int)vehicles.size()) {
255 fractionLeft[odID] = fraction;
264 std::string personDepartPos = oc.
isSet(
"departpos") ? oc.
getString(
"departpos") :
"random";
265 std::string personArrivalPos = oc.
isSet(
"arrivalpos") ? oc.
getString(
"arrivalpos") :
"random";
266 for (std::vector<ODVehicle>::reverse_iterator i = vehicles.rbegin(); i != vehicles.rend() && (*i).depart == t; ++i) {
278 }
else if (persontrips) {
298 while (vehicles.size() != 0 && vehicles.back().depart == t) {
301 if (!vehicles.empty()) {
302 t = vehicles.back().depart;
304 if (next !=
myContainer.end() && (t > (*next)->begin || vehicles.empty())) {
307 if (next ==
myContainer.end() && vehicles.empty()) {
317 const std::string& prefix,
318 bool asProbability,
bool pedestrians,
bool persontrips,
319 const std::string& modes) {
327 const ODCell*
const c = *i;
328 if (c->
end > begin && c->
begin < end) {
330 if (probability <= 0) {
337 if (!asProbability) {
340 if (probability > 1) {
341 WRITE_WARNING(
"Flow density of " +
toString(probability) +
" vehicles per second, cannot be represented with a simple probability. Falling back to even spacing.");
354 }
else if (persontrips) {
357 if (!asProbability) {
360 if (probability > 1) {
361 WRITE_WARNING(
"Flow density of " +
toString(probability) +
" vehicles per second, cannot be represented with a simple probability. Falling back to even spacing.");
383 if (!asProbability) {
386 if (probability > 1) {
387 WRITE_WARNING(
"Flow density of " +
toString(probability) +
" vehicles per second, cannot be represented with a simple probability. Falling back to even spacing.");
406 const std::string line = lr.
readLine();
407 if (line[0] !=
'*') {
417 if (time.find(
'.') == std::string::npos) {
420 std::string hours = time.substr(0, time.find(
'.'));
421 std::string minutes = time.substr(time.find(
'.') + 1);
426 std::pair<SUMOTime, SUMOTime>
438 throw ProcessError(
"Broken period definition '" + line +
"'.");
440 throw ProcessError(
"Broken period definition '" + line +
"'.");
459 std::string vehType,
bool matrixHasVehType) {
463 if (matrixHasVehType) {
471 std::pair<SUMOTime, SUMOTime> times =
readTime(lr);
482 std::vector<std::string> names;
483 while ((
int)names.size() != numDistricts && lr.
hasMore()) {
487 names.push_back(st2.
next());
495 for (std::vector<std::string>::iterator si = names.begin(); si != names.end(); ++si) {
496 std::vector<std::string>::iterator di = names.begin();
501 throw ProcessError(
"Missing line for district " + (*si) +
".");
503 if (line.length() == 0) {
509 assert(di != names.end());
511 if (vehNumber != 0) {
512 add(vehNumber, begin, end, *si, *di, vehType);
514 if (di == names.end()) {
515 throw ProcessError(
"More entries than districts found.");
520 throw ProcessError(
"Not numeric vehicle number in line '" + line +
"'.");
525 }
while (di != names.end());
533 std::string vehType,
bool matrixHasVehType) {
537 if (matrixHasVehType) {
546 std::pair<SUMOTime, SUMOTime> times =
readTime(lr);
556 if (line.length() == 0) {
560 if (st2.
size() == 0) {
564 std::string sourceD = st2.
next();
565 std::string destD = st2.
next();
567 if (vehNumber != 0) {
568 add(vehNumber, begin, end, sourceD, destD, vehType);
571 throw ProcessError(
"Missing at least one information in line '" + line +
"'.");
573 throw ProcessError(
"Not numeric vehicle number in line '" + line +
"'.");
601 const std::vector<double>& times = ps.
getVals();
602 for (
int i = 0; i < (int)times.size() - 1; ++i) {
610 newCells.push_back(ncell);
619 for (std::vector<ODCell*>::iterator i = oldCells.begin(); i != oldCells.end(); ++i) {
620 std::vector<ODCell*> newCells;
622 copy(newCells.begin(), newCells.end(), back_inserter(
myContainer));
630 std::vector<std::string> files = oc.
getStringVector(
"od-matrix-files");
631 for (std::vector<std::string>::iterator i = files.begin(); i != files.end(); ++i) {
638 if (type.find(
';') != std::string::npos) {
639 type = type.substr(0, type.find(
';'));
642 if (type.length() > 1 && type[1] ==
'V') {
644 if (type.find(
'N') != std::string::npos) {
645 throw ProcessError(
"'" + *i +
"' does not contain the needed information about the time described.");
648 }
else if (type.length() > 1 && type[1] ==
'O') {
650 if (type.find(
'N') != std::string::npos) {
651 throw ProcessError(
"'" + *i +
"' does not contain the needed information about the time described.");
655 throw ProcessError(
"'" + *i +
"' uses an unknown matrix type '" + type +
"'.");
658 std::vector<std::string> amitranFiles = oc.
getStringVector(
"od-amitran-files");
659 for (std::vector<std::string>::iterator i = amitranFiles.begin(); i != amitranFiles.end(); ++i) {
661 throw ProcessError(
"Could not access matrix file '" + *i +
"' to load.");
676 std::vector<std::string> routeFiles = oc.
getStringVector(
"route-files");
677 for (std::vector<std::string>::iterator i = routeFiles.begin(); i != routeFiles.end(); ++i) {
679 throw ProcessError(
"Could not access route file '" + *i +
"' to load.");
694 if (timelineDayInHours) {
695 if (def.size() != 24) {
696 throw ProcessError(
"Assuming 24 entries for a day timeline, but got " +
toString(def.size()) +
".");
698 for (
int chour = 0; chour < 24; ++chour) {
701 result.
add(24 * 3600., 0.);
703 for (
int i = 0; i < (int)def.size(); i++) {
705 if (st2.
size() != 2) {
706 throw ProcessError(
"Broken time line definition: missing a value in '" + def[i] +
"'.");
#define WRITE_WARNING(msg)
#define PROGRESS_DONE_MESSAGE()
#define PROGRESS_FAILED_MESSAGE()
#define PROGRESS_BEGIN_MESSAGE(msg)
std::string time2string(SUMOTime t)
convert SUMOTime to string
SUMOTime string2time(const std::string &r)
convert string to SUMOTime
@ SUMO_TAG_FLOW
a flow definitio nusing a from-to edges instead of a route (used by router)
@ SUMO_TAG_TRIP
a single trip definition (used by router)
@ SUMO_ATTR_BEGIN
weights: time range begin
@ SUMO_ATTR_END
weights: time range end
std::string toString(const T &t, std::streamsize accuracy=gPrecision)
static bool isReadable(std::string path)
Checks whether the given file is readable.
Retrieves a file linewise and reports the lines to a handler.
bool good() const
Returns the information whether the stream is readable.
bool readLine(LineHandler &lh)
Reads a single (the next) line from the file and reports it to the given LineHandler.
bool hasMore() const
Returns whether another line may be read (the file was not read completely)
std::string getFileName() const
Returns the name of the used file.
T get(const std::string &id) const
Retrieves an item.
An XML-Handler for districts.
A container for districts.
std::string getRandomSourceFromDistrict(const std::string &name) const
Returns the id of a random source from the named district.
std::string getRandomSinkFromDistrict(const std::string &name) const
Returns the id of a random sink from the named district.
int sourceNumber() const
Returns the number of sources.
int sinkNumber() const
Returns the number of sinks.
Used for sorting the cells by the begin time they describe.
Used for sorting vehicles by their departure (latest first)
double getNumLoaded() const
Returns the number of loaded vehicles.
void readV(LineReader &lr, double scale, std::string vehType, bool matrixHasVehType)
read a VISUM-matrix with the V Format
bool add(double vehicleNumber, SUMOTime begin, SUMOTime end, const std::string &origin, const std::string &destination, const std::string &vehicleType, const bool originIsEdge=false, const bool destinationIsEdge=false)
Builds a single cell from the given values, verifying them.
double readFactor(LineReader &lr, double scale)
ODMatrix(const ODDistrictCont &dc)
Constructor.
double computeDeparts(ODCell *cell, int &vehName, std::vector< ODVehicle > &into, const bool uniform, const bool differSourceSink, const std::string &prefix)
Computes the vehicle departs stored in the given cell and saves them in "into".
SUMOTime myBegin
parsed time bounds
SUMOTime parseSingleTime(const std::string &time)
void writeFlows(const SUMOTime begin, const SUMOTime end, OutputDevice &dev, const bool noVtype, const std::string &prefix, bool asProbability=false, bool pedestrians=false, bool persontrips=false, const std::string &modes="")
Writes the flows stored in the matrix.
void applyCurve(const Distribution_Points &ps)
Splits the stored cells dividing them on the given time line.
std::map< const std::pair< const std::string, const std::string >, std::vector< ODCell * > > myShortCut
The loaded cells indexed by origin and destination.
std::pair< SUMOTime, SUMOTime > readTime(LineReader &lr)
void readO(LineReader &lr, double scale, std::string vehType, bool matrixHasVehType)
read a VISUM-matrix with the O Format
std::set< std::string > myMissingDistricts
The missing districts already warned about.
void write(SUMOTime begin, const SUMOTime end, OutputDevice &dev, const bool uniform, const bool differSourceSink, const bool noVtype, const std::string &prefix, const bool stepLog, bool pedestrians, bool persontrips, const std::string &modes)
Writes the vehicles stored in the matrix assigning the sources and sinks.
const ODDistrictCont & myDistricts
The districts to retrieve sources/sinks from.
Distribution_Points parseTimeLine(const std::vector< std::string > &def, bool timelineDayInHours)
split the given timeline
void writeDefaultAttrs(OutputDevice &dev, const bool noVtype, const ODCell *const cell)
Helper function for flow and trip output writing the depart and arrival attributes.
double myNumLoaded
Number of loaded vehicles.
double myNumWritten
Number of written vehicles.
void loadMatrix(OptionsCont &oc)
read a matrix in one of several formats
void loadRoutes(OptionsCont &oc, SUMOSAXHandler &handler)
read SUMO routes
double getNumWritten() const
Returns the number of written vehicles.
std::vector< ODCell * > myContainer
The loaded cells.
double myNumDiscarded
Number of discarded vehicles.
double getNumDiscarded() const
Returns the number of discarded vehicles.
std::string getNextNonCommentLine(LineReader &lr)
A storage for options typed value containers)
bool isSet(const std::string &name, bool failOnNonExistant=true) const
Returns the information whether the named option is set.
double getFloat(const std::string &name) const
Returns the double-value of the named option (only for Option_Float)
std::string getString(const std::string &name) const
Returns the string-value of the named option (only for Option_String)
const StringVector & getStringVector(const std::string &name) const
Returns the list of string-value of the named option (only for Option_StringVector)
static OptionsCont & getOptions()
Retrieves the options.
Static storage of an output device and its base (abstract) implementation.
OutputDevice & openTag(const std::string &xmlElement)
Opens an XML tag.
OutputDevice & writeAttr(const SumoXMLAttr attr, const T &val)
writes a named attribute
bool closeTag(const std::string &comment="")
Closes the most recently opened tag and optionally adds a comment.
void setPrecision(int precision=gPrecision)
Sets the precison or resets it to default.
A complete router's route.
static double rand(std::mt19937 *rng=nullptr)
Returns a random real number in [0, 1)
double getOverallProb() const
Return the sum of the probabilites assigned to the members.
const std::vector< double > & getProbs() const
Returns the probabilities assigned to the members of the distribution.
const std::vector< T > & getVals() const
Returns the members of the distribution.
bool add(T val, double prob, bool checkDuplicates=true)
Adds a value with an assigned probability to the distribution.
SAX-handler base for SUMO-files.
int size() const
returns the number of existing substrings
static const int WHITECHARS
identifier for splitting the given string at all whitespace characters
bool hasNext()
returns the information whether further substrings exist
std::string next()
returns the next substring when it exists. Otherwise the behaviour is undefined
static double toDouble(const std::string &sData)
converts a string into the double value described by it by calling the char-type converter
static std::string prune(const std::string &str)
Removes trailing and leading whitechars.
static int toInt(const std::string &sData)
converts a string into the integer value described by it by calling the char-type converter,...
static bool runParser(GenericSAXHandler &handler, const std::string &file, const bool isNet=false, const bool isRoute=false)
Runs the given handler on the given file; returns if everything's ok.
A single O/D-matrix cell.
std::string destination
Name of the destination district.
std::map< SUMOTime, std::vector< std::string > > departures
mapping of departure times to departing vehicles, if already fixed
std::string vehicleType
Name of the vehicle type.
std::string origin
Name of the origin district.
double vehicleNumber
The number of vehicles.
bool originIsEdge
the origin "district" is an edge id
SUMOTime end
The end time this cell describes.
SUMOTime begin
The begin time this cell describes.
bool destinationIsEdge
the destination "district" is an edge id
An internal representation of a single vehicle.
SUMOTime depart
The departure time of the vehicle.
std::string from
The edge the vehicles shall start at.
ODCell * cell
The cell of the ODMatrix which generated the vehicle.
std::string to
The edge the vehicles shall end at.
std::string id
The id of the vehicle.