Marmote Core
The project aims at realizing the prototype of a software environment dedicated to modeling with Markov chains.
transitionStructure.h
1 /* Marmote is free software: you can redistribute it and/or modify
2 it under the terms of the GNU General Public License as published by
3 the Free Software Foundation, either version 3 of the License, or
4 (at your option) any later version.
5 
6 Marmote is distributed in the hope that it will be useful,
7 but WITHOUT ANY WARRANTY; without even the implied warranty of
8 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
9 GNU General Public License for more details.
10 
11 You should have received a copy of the GNU General Public License
12 along with Marmote. If not, see <http://www.gnu.org/licenses/>.
13 
14 Copyright 2015 Alain Jean-Marie, Jean-Michel Fourneau, Jean-Marc Vincent, Issam Rabhi */
15 
16 #ifndef TRANSITIONSTRUCTURE_H
17 #define TRANSITIONSTRUCTURE_H
18 
19 #include "../marmoteConstants.h"
20 #include "../Distribution/discreteDistribution.h"
21 #include <stdio.h>
22 #include <string>
23 
33 {
34  protected:
35  timeType _type;
37  long int _size;
38  double _uniformizationRate;
40  public:
41  // constructors
42  // transitionStructure(int);
47  virtual ~transitionStructure() {};
48 
49  public:
50  //accessors
56  int size() { return _size; };
62  timeType type() { return _type; };
75  void setType(timeType t) { _type = t; };
81  void setUniformizationRate(double rate) { _uniformizationRate = rate; };
90  virtual bool setEntry(int i, int j, double val) = 0;
99  virtual double getEntry(int i, int j) = 0;
107  virtual int getNbElts(int i) = 0;
116  virtual int getCol(int i, int k) = 0;
125  virtual double getEntryByCol(int i, int k) = 0;
132  virtual discreteDistribution* getTransDistrib(int i) = 0; // transitions from some state and their probas
133 
142  bool readEntry(FILE *input);
149  virtual double rowSum(int i) = 0;
150 
151  public:
157  virtual transitionStructure* copy() = 0;
166  virtual transitionStructure* uniformize() = 0;
176  virtual void evaluateMeasure(double* d,double* res) {
177  fprintf( stderr, "Warning in transitionStructure::evaluateDistribution: not implemented.\n" ); };
185  fprintf( stderr, "Warning in transitionStructure::evaluateDistribution: not implemented.\n" ); };
194  virtual void evaluateValue(double* v, double* res) = 0;
195 
196 protected:
197  // specific technical methods
209  int consolidate(int i, int* destinations, double* values);
210 
211 public:
218  virtual void write(FILE* out, std::string format) = 0;
219 
220 };
221 
222 #endif // TRANSITIONSTRUCTURE_H
virtual discreteDistribution * getTransDistrib(int i)=0
Method to get the transition from some state as a probability distribution.
virtual void evaluateMeasure(double *d, double *res)
Computing the action of the transition structure on some measure, the measure being represented as a ...
Definition: transitionStructure.h:176
int consolidate(int i, int *destinations, double *values)
Method to consolidate (aggregate) transition probabilities from a given state. The method returns the...
Definition: transitionStructure.cpp:36
timeType _type
the time type of the structure: discrete or continuous.
Definition: transitionStructure.h:20
virtual transitionStructure * uniformize()=0
Uniformizing a transition structure. The structure should be of continuous time type. The resulting one will be of discrete time type. If uniformization fails, a NULL pointer should be returned. If the origin structure is already of discrete-time type, a copy should be returned.
virtual int getCol(int i, int k)=0
Method to get the number of the state corresponding to transition number k in the list of possible tr...
virtual double getEntryByCol(int i, int k)=0
Method to get the value attached to transition number k in the list of possible transitions from some...
virtual int getNbElts(int i)=0
Method to get the number of non-zero entries in a transition from some state. This can be seen as the...
double _uniformizationRate
value related to the uniformization procedure.
Definition: transitionStructure.h:23
virtual void write(FILE *out, std::string format)=0
Output method for a transition structure.
long int _size
size of the state space.
Definition: transitionStructure.h:22
virtual void evaluateValue(double *v, double *res)=0
Computing the action of the transition structure on some vector of values. This corresponds to the mu...
bool readEntry(FILE *input)
Reading from a file, and adding an element to the matrix. The field must be in the form "row column v...
Definition: transitionStructure.cpp:19
virtual double rowSum(int i)=0
Calculating the sum of values corresponding to transitions from some state.
virtual bool setEntry(int i, int j, double val)=0
Method to set the value associated with some transition.
Abstract class for transition structures. These are structures which describe transitions to one stat...
Definition: transitionStructure.h:17
void setUniformizationRate(double rate)
Write accessor for the uniformization rate.
Definition: transitionStructure.h:81
virtual ~transitionStructure()
Destructor of the class.
Definition: transitionStructure.h:47
timeType type()
Read accessor for the time type.
Definition: transitionStructure.h:62
virtual double getEntry(int i, int j)=0
Method to get the value associated with some transition. When state parameters are out of bounds...
virtual void evaluateMeasure(discreteDistribution *d, discreteDistribution *res)
Computing the action of the transition structure on some probability distribution.
Definition: transitionStructure.h:184
void setType(timeType t)
Write accessor for the time type.
Definition: transitionStructure.h:75
int size()
Read accessor for the size of the state space.
Definition: transitionStructure.h:56
virtual transitionStructure * copy()=0
Copying a transition structure.
The general discrete distribution with finite support.
Definition: discreteDistribution.h:25
double uniformizationRate()
Read accessor for the uniformization rate. Relevant mostly for discrete-time structures created by un...
Definition: transitionStructure.h:69