From c02f4e2adc2d8f8d09e6f8d930c6ea654a2215ba Mon Sep 17 00:00:00 2001 From: dehnert Date: Thu, 14 Nov 2013 16:17:46 +0100 Subject: [PATCH] Markov Automata transitions can now be parsed. Next up: a parser that combines transition and label parsing for Markov automata. Former-commit-id: 77db051f1fa1c2d171963313bf4a4c4aa2c80342 --- .../MarkovAutomataSparseTransitionParser.cpp | 207 +++++++++++++----- .../MarkovAutomataSparseTransitionParser.h | 20 +- src/storm.cpp | 3 - 3 files changed, 165 insertions(+), 65 deletions(-) diff --git a/src/parser/MarkovAutomataSparseTransitionParser.cpp b/src/parser/MarkovAutomataSparseTransitionParser.cpp index e8f7d7002..0293fb2e6 100644 --- a/src/parser/MarkovAutomataSparseTransitionParser.cpp +++ b/src/parser/MarkovAutomataSparseTransitionParser.cpp @@ -1,5 +1,7 @@ #include "MarkovAutomataSparseTransitionParser.h" +#include "src/settings/Settings.h" + namespace storm { namespace parser { @@ -11,16 +13,13 @@ namespace storm { // Skip the format hint. buf = storm::parser::forwardToNextLine(buf, lineEndings); - /* - * Read all transitions. - */ + // Now read the transitions. int_fast64_t source, target = -1; int_fast64_t lastsource = -1; bool encounteredEOF = false; + bool stateHasMarkovianChoice = false; while (buf[0] != '\0' && !encounteredEOF) { - /* - * Read source state and choice. - */ + // At the current point, the next thing to read is the source state of the next choice to come. source = checked_strtol(buf, &buf); // Check if we encountered a state index that is bigger than all previously seen ones and record it if necessary. @@ -36,6 +35,11 @@ namespace storm { ++result.numberOfChoices; } + // If we have moved to the next state, we need to clear the flag that stores whether or not the source has a Markovian choice. + if (source != lastsource) { + stateHasMarkovianChoice = false; + } + // Record that the current source was the last source. lastsource = source; @@ -63,24 +67,12 @@ namespace storm { isMarkovianChoice = false; } - // Proceed to beginning of next line. - switch (lineEndings) { - case SupportedLineEndingsEnum::SlashN: - buf += strcspn(buf, " \t\n"); - break; - case SupportedLineEndingsEnum::SlashR: - buf += strcspn(buf, " \t\r"); - break; - case SupportedLineEndingsEnum::SlashRN: - buf += strcspn(buf, " \t\r\n"); - break; - default: - case storm::parser::SupportedLineEndingsEnum::Unsupported: - // This line will never be reached as the parser would have thrown already. - throw; - break; + if (isMarkovianChoice && stateHasMarkovianChoice) { + LOG4CPLUS_ERROR(logger, "The state " << source << " has multiple Markovian choices."); + throw storm::exceptions::WrongFormatException() << "The state " << source << " has multiple Markovian choices."; } - buf = trimWhitespaces(buf); + + buf = forwardToNextLine(buf, lineEndings); // Now that we have the source state and the information whether or not the current choice is probabilistic or Markovian, we need to read the list of successors and the probabilities/rates. bool hasSuccessorState = false; @@ -132,24 +124,7 @@ namespace storm { // As we found a new successor, we need to increase the number of nonzero entries. ++result.numberOfNonzeroEntries; - // Proceed to beginning of next line. - switch (lineEndings) { - case SupportedLineEndingsEnum::SlashN: - buf += strcspn(buf, " \t\n"); - break; - case SupportedLineEndingsEnum::SlashR: - buf += strcspn(buf, " \t\r"); - break; - case SupportedLineEndingsEnum::SlashRN: - buf += strcspn(buf, " \t\r\n"); - break; - default: - case storm::parser::SupportedLineEndingsEnum::Unsupported: - // This line will never be reached as the parser would have thrown already. - throw; - break; - } - buf = trimWhitespaces(buf); + buf = forwardToNextLine(buf, lineEndings); } else { // If it was not a "*", we have to assume that we encountered the beginning of a new choice definition. In this case, we don't move the pointer // to the buffer, because we still need to read the new source state. @@ -157,37 +132,159 @@ namespace storm { } } while (!encounteredEOF && !encounteredNewDistribution); } - - exit(-1); - + return result; } MarkovAutomataSparseTransitionParser::ResultType MarkovAutomataSparseTransitionParser::secondPass(char* buf, SupportedLineEndingsEnum lineEndings, FirstPassResult const& firstPassResult, RewardMatrixInformationStruct* rewardMatrixInformation) { - return ResultType(); + ResultType result(firstPassResult); + + bool fixDeadlocks = storm::settings::Settings::getInstance()->isSet("fixDeadlocks"); + + // Skip the format hint. + buf = storm::parser::forwardToNextLine(buf, lineEndings); + + // Now read the transitions. + int_fast64_t source, target = -1; + int_fast64_t lastsource = -1; + bool encounteredEOF = false; + uint_fast64_t currentChoice = 0; + while (buf[0] != '\0' && !encounteredEOF) { + // At the current point, the next thing to read is the source state of the next choice to come. + source = checked_strtol(buf, &buf); + + // If we have skipped some states, we need to insert self-loops if requested. + if (source > lastsource + 1) { + if (fixDeadlocks) { + for (uint_fast64_t index = lastsource + 1; index < source; ++source) { + result.nondeterministicChoiceIndices[index] = currentChoice; + result.transitionMatrix.addNextValue(currentChoice, index, 1); + ++currentChoice; + } + } else { + LOG4CPLUS_ERROR(logger, "Found deadlock states (e.g. " << lastsource + 1 << ") during parsing. Please fix them or set the appropriate flag."); + throw storm::exceptions::WrongFormatException() << "Found deadlock states (e.g. " << lastsource + 1 << ") during parsing. Please fix them or set the appropriate flag."; + } + } + + if (source != lastsource) { + // If we skipped to a new state we need to record the beginning of the choices in the nondeterministic choice indices. + result.nondeterministicChoiceIndices[source] = currentChoice; + } + + // Record that the current source was the last source. + lastsource = source; + + char actionNameBuffer[20]; +#ifdef WINDOWS + int length = sscanf_s(buf, "%20s\n", actionNameBuffer, 20); +#else + int length = sscanf(buf, "%20s\n", actionNameBuffer); +#endif + + // If the number of arguments filled is not one, there was an error. + if (length != 1) { + LOG4CPLUS_ERROR(logger, "Parsing error."); + throw storm::exceptions::WrongFormatException() << "Parsing error."; + } else { + // If the action name was parsed successfully, we need to move by the corresponding number of characters. + buf += strlen(actionNameBuffer); + } + + // Depending on the action name, the choice is either a probabilitic one or a markovian one. + bool isMarkovianChoice = false; + if (strcmp(actionNameBuffer, "!") == 0) { + isMarkovianChoice = true; + + // Mark the current choice as a Markovian one. + result.markovianChoices.set(currentChoice, true); + } else { + isMarkovianChoice = false; + } + + buf = forwardToNextLine(buf, lineEndings); + + // Now that we have the source state and the information whether or not the current choice is probabilistic or Markovian, we need to read the list of successors and the probabilities/rates. + bool hasSuccessorState = false; + bool encounteredNewDistribution = false; + + // At this point, we need to check whether there is an additional successor or we have reached the next choice for the same or a different state. + do { + // Now parse the next symbol to see whether there is another target state for the current choice + // or not. + char star[1]; +#ifdef WINDOWS + length = sscanf_s(buf, "%1s\n", star, 1); +#else + length = sscanf(buf, "%1s\n", star); +#endif + + // If the number of arguments filled is not one, there was an error. + if (length == EOF) { + if (!hasSuccessorState) { + LOG4CPLUS_ERROR(logger, "Premature end-of-file. Expected at least one successor state for state " << source << " under action " << actionNameBuffer << "."); + throw storm::exceptions::WrongFormatException() << "Premature end-of-file. Expected at least one successor state for state " << source << " under action " << actionNameBuffer << "."; + } else { + // If there was at least one successor for the current choice, this is legal and we need to move on. + encounteredEOF = true; + } + } else if (length != 1) { + LOG4CPLUS_ERROR(logger, "Parsing error."); + throw storm::exceptions::WrongFormatException() << "Parsing error."; + } else if (strcmp(star, "*") == 0) { + // We need to record that we found at least one successor state for the current choice. + hasSuccessorState = true; + + // As we have encountered a "*", we know that there is an additional successor state for the current choice. + buf += strlen(star); + + // Now we need to read the successor state and check if we already saw a higher state index. + target = checked_strtol(buf, &buf); + + // And the corresponding probability/rate. + double val = checked_strtod(buf, &buf); + + // Record the value as well as the exit rate in case of a Markovian choice. + result.transitionMatrix.addNextValue(currentChoice, target, val); + if (isMarkovianChoice) { + result.exitRates[currentChoice] += val; + } + + buf = forwardToNextLine(buf, lineEndings); + } else { + // If it was not a "*", we have to assume that we encountered the beginning of a new choice definition. In this case, we don't move the pointer + // to the buffer, because we still need to read the new source state. + encounteredNewDistribution = true; + } + } while (!encounteredEOF && !encounteredNewDistribution); + + ++currentChoice; + } + + // As we have added all entries at this point, we need to finalize the matrix. + result.transitionMatrix.finalize(); + + // Put a sentinel element at the end. + result.nondeterministicChoiceIndices[firstPassResult.highestStateIndex + 1] = currentChoice; + + return result; } MarkovAutomataSparseTransitionParser::ResultType MarkovAutomataSparseTransitionParser::parseMarkovAutomataTransitions(std::string const& filename, RewardMatrixInformationStruct* rewardMatrixInformation) { - /* - * Enforce locale where decimal point is '.'. - */ + // Set the locale to correctly recognize floating point numbers. setlocale(LC_NUMERIC, "C"); if (!fileExistsAndIsReadable(filename.c_str())) { LOG4CPLUS_ERROR(logger, "Error while parsing " << filename << ": File does not exist or is not readable."); - throw storm::exceptions::WrongFormatException(); + throw storm::exceptions::WrongFormatException() << "Error while parsing " << filename << ": File does not exist or is not readable."; } bool isRewardFile = rewardMatrixInformation != nullptr; - /* - * Find out about the used line endings. - */ + // Determine used line endings. SupportedLineEndingsEnum lineEndings = findUsedLineEndings(filename, true); - /* - * Open file. - */ + // Open file and prepare pointer to buffer. MappedFile file(filename.c_str()); char* buf = file.data; diff --git a/src/parser/MarkovAutomataSparseTransitionParser.h b/src/parser/MarkovAutomataSparseTransitionParser.h index ed7f06ae2..933208b19 100644 --- a/src/parser/MarkovAutomataSparseTransitionParser.h +++ b/src/parser/MarkovAutomataSparseTransitionParser.h @@ -10,13 +10,6 @@ namespace storm { class MarkovAutomataSparseTransitionParser { public: - struct ResultType { - storm::storage::SparseMatrix transitionMatrix; - std::vector nondeterministicChoiceIndices; - storm::storage::BitVector markovianChoices; - std::vector exitRates; - }; - struct FirstPassResult { FirstPassResult() : numberOfNonzeroEntries(0), highestStateIndex(0), numberOfChoices(0) { @@ -28,6 +21,19 @@ namespace storm { uint_fast64_t numberOfChoices; }; + struct ResultType { + + ResultType(FirstPassResult const& firstPassResult) : transitionMatrix(firstPassResult.numberOfChoices, firstPassResult.highestStateIndex + 1), nondeterministicChoiceIndices(firstPassResult.highestStateIndex + 2), markovianChoices(firstPassResult.numberOfChoices), exitRates(firstPassResult.numberOfChoices) { + transitionMatrix.initialize(firstPassResult.numberOfNonzeroEntries); + // Intentionally left empty. + } + + storm::storage::SparseMatrix transitionMatrix; + std::vector nondeterministicChoiceIndices; + storm::storage::BitVector markovianChoices; + std::vector exitRates; + }; + static ResultType parseMarkovAutomataTransitions(std::string const& filename, RewardMatrixInformationStruct* rewardMatrixInformation = nullptr); private: diff --git a/src/storm.cpp b/src/storm.cpp index 52c157d4d..19db92fca 100644 --- a/src/storm.cpp +++ b/src/storm.cpp @@ -415,9 +415,6 @@ int main(const int argc, const char* argv[]) { std::string const chosenTransitionSystemFile = s->getOptionByLongName("explicit").getArgument(0).getValueAsString(); std::string const chosenLabelingFile = s->getOptionByLongName("explicit").getArgument(1).getValueAsString(); - // FIXME: Remove after testing. - storm::parser::MarkovAutomataSparseTransitionParser::ResultType result = storm::parser::MarkovAutomataSparseTransitionParser::parseMarkovAutomataTransitions(chosenTransitionSystemFile, nullptr); - std::string chosenStateRewardsFile = ""; if (s->isSet("stateRewards")) { chosenStateRewardsFile = s->getOptionByLongName("stateRewards").getArgument(0).getValueAsString();