From 9775f7ae7edb06dda347aca501bc47548d87020f Mon Sep 17 00:00:00 2001 From: gereon Date: Mon, 24 Dec 2012 13:37:35 +0100 Subject: [PATCH] Added CTMC model, doing c&p from dtmc. I have no idea, if this is what we actually want... :-) --- src/models/Ctmc.h | 189 ++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 189 insertions(+) create mode 100644 src/models/Ctmc.h diff --git a/src/models/Ctmc.h b/src/models/Ctmc.h new file mode 100644 index 000000000..5101cd556 --- /dev/null +++ b/src/models/Ctmc.h @@ -0,0 +1,189 @@ +/* + * Ctmc.h + * + * Created on: 14.11.2012 + * Author: Christian Dehnert + */ + +#ifndef MRMC_MODELS_CTMC_H_ +#define MRMC_MODELS_CTMC_H_ + +#include +#include +#include +#include + +#include "AtomicPropositionsLabeling.h" +#include "GraphTransitions.h" +#include "src/storage/SquareSparseMatrix.h" +#include "src/exceptions/InvalidArgumentException.h" + +namespace mrmc { + +namespace models { + +/*! + * This class represents a discrete-time Markov chain (DTMC) whose states are + * labeled with atomic propositions. + */ +template +class Ctmc { + +public: + //! Constructor + /*! + * Constructs a CTMC object from the given transition rate matrix and + * the given labeling of the states. + * @param rateMatrix The transition rate function of the + * CTMC given by a matrix. + * @param stateLabeling The labeling that assigns a set of atomic + * propositions to each state. + */ + Ctmc(std::shared_ptr> rateMatrix, + std::shared_ptr stateLabeling, + std::shared_ptr> stateRewards = nullptr, + std::shared_ptr> transitionRewardMatrix = nullptr) + : rateMatrix(rateMatrix), stateLabeling(stateLabeling), + stateRewards(stateRewards), transitionRewardMatrix(transitionRewardMatrix), + backwardTransitions(nullptr) { + } + + //! Copy Constructor + /*! + * Copy Constructor. Performs a deep copy of the given CTMC. + * @param ctmc A reference to the CTMC that is to be copied. + */ + Ctmc(const Ctmc &ctmc) : rateMatrix(ctmc.rateMatrix), + stateLabeling(ctmc.stateLabeling), stateRewards(ctmc.stateRewards), + transitionRewardMatrix(ctmc.transitionRewardMatrix) { + if (ctmc.backardTransitions != nullptr) { + this->backwardTransitions = new mrmc::models::GraphTransitions(*ctmc.backwardTransitions); + } + } + + //! Destructor + /*! + * Destructor. Frees the matrix and labeling associated with this CTMC. + */ + ~Ctmc() { + if (this->backwardTransitions != nullptr) { + delete this->backwardTransitions; + } + } + + /*! + * Returns the state space size of the CTMC. + * @return The size of the state space of the CTMC. + */ + uint_fast64_t getNumberOfStates() const { + return this->rateMatrix->getRowCount(); + } + + /*! + * Returns the number of (non-zero) transitions of the CTMC. + * @return The number of (non-zero) transitions of the CTMC. + */ + uint_fast64_t getNumberOfTransitions() const { + return this->rateMatrix->getNonZeroEntryCount(); + } + + /*! + * Returns a bit vector in which exactly those bits are set to true that + * correspond to a state labeled with the given atomic proposition. + * @param ap The atomic proposition for which to get the bit vector. + * @return A bit vector in which exactly those bits are set to true that + * correspond to a state labeled with the given atomic proposition. + */ + mrmc::storage::BitVector* getLabeledStates(std::string ap) const { + return this->stateLabeling->getAtomicProposition(ap); + } + + /*! + * Returns a pointer to the matrix representing the transition probability + * function. + * @return A pointer to the matrix representing the transition probability + * function. + */ + std::shared_ptr> getTransitionRateMatrix() const { + return this->rateMatrix; + } + + /*! + * Returns a pointer to the matrix representing the transition rewards. + * @return A pointer to the matrix representing the transition rewards. + */ + std::shared_ptr> getTransitionRewardMatrix() const { + return this->transitionRewardMatrix; + } + + /*! + * Returns a pointer to the vector representing the state rewards. + * @return A pointer to the vector representing the state rewards. + */ + std::shared_ptr> getStateRewards() const { + return this->stateRewards; + } + + /*! + * + */ + std::set const getPropositionsForState(uint_fast64_t const &state) const { + return stateLabeling->getPropositionsForState(state); + } + + /*! + * Retrieves a reference to the backwards transition relation. + * @return A reference to the backwards transition relation. + */ + mrmc::models::GraphTransitions& getBackwardTransitions() { + if (this->backwardTransitions == nullptr) { + this->backwardTransitions = new mrmc::models::GraphTransitions(this->probabilityMatrix, false); + } + return *this->backwardTransitions; + } + + /*! + * Prints information about the model to the specified stream. + * @param out The stream the information is to be printed to. + */ + void printModelInformationToStream(std::ostream& out) const { + out << "-------------------------------------------------------------- " + << std::endl; + out << "Model type: \t\tCTMC" << std::endl; + out << "States: \t\t" << this->getNumberOfStates() << std::endl; + out << "Transitions: \t\t" << this->getNumberOfTransitions() << std::endl; + this->stateLabeling->printAtomicPropositionsInformationToStream(out); + out << "Size in memory: \t" + << (this->rateMatrix->getSizeInMemory() + + this->stateLabeling->getSizeInMemory() + + sizeof(*this))/1024 << " kbytes" << std::endl; + out << "-------------------------------------------------------------- " + << std::endl; + } + +private: + + /*! A matrix representing the transition rate function of the CTMC. */ + std::shared_ptr> rateMatrix; + + /*! The labeling of the states of the CTMC. */ + std::shared_ptr stateLabeling; + + /*! The state-based rewards of the CTMC. */ + std::shared_ptr> stateRewards; + + /*! The transition-based rewards of the CTMC. */ + std::shared_ptr> transitionRewardMatrix; + + /*! + * A data structure that stores the predecessors for all states. This is + * needed for backwards directed searches. + */ + mrmc::models::GraphTransitions* backwardTransitions; +}; + +} // namespace models + +} // namespace mrmc + +#endif /* MRMC_MODELS_DTMC_H_ */