[a9b86d] | 1 | /*
|
---|
| 2 | * MatrixContainer.hpp
|
---|
| 3 | *
|
---|
| 4 | * Created on: Sep 15, 2011
|
---|
| 5 | * Author: heber
|
---|
| 6 | */
|
---|
| 7 |
|
---|
| 8 | #ifndef MATRIXCONTAINER_HPP_
|
---|
| 9 | #define MATRIXCONTAINER_HPP_
|
---|
| 10 |
|
---|
| 11 | // include config.h
|
---|
| 12 | #ifdef HAVE_CONFIG_H
|
---|
| 13 | #include <config.h>
|
---|
| 14 | #endif
|
---|
| 15 |
|
---|
[9758f7] | 16 | #include <vector>
|
---|
[a9b86d] | 17 |
|
---|
[942906] | 18 | #include "boost/serialization/vector.hpp"
|
---|
| 19 |
|
---|
[a9b86d] | 20 | /** Contains a sequence of matrices parsed from file.
|
---|
| 21 | *
|
---|
| 22 | */
|
---|
| 23 | class MatrixContainer {
|
---|
| 24 | public:
|
---|
[9758f7] | 25 | typedef std::vector<std::vector<double> > MatrixArray;
|
---|
| 26 | typedef std::vector<MatrixArray> ArrayOfMatrices;
|
---|
| 27 | typedef std::vector<std::string> StringVector;
|
---|
| 28 | typedef std::vector<int> IntVector;
|
---|
| 29 | typedef std::vector<IntVector> VectorOfIntVectors;
|
---|
| 30 |
|
---|
| 31 | ArrayOfMatrices Matrix;
|
---|
| 32 | VectorOfIntVectors Indices;
|
---|
| 33 | StringVector Header;
|
---|
[a9b86d] | 34 | int MatrixCounter;
|
---|
[9758f7] | 35 | IntVector RowCounter;
|
---|
| 36 | IntVector ColumnCounter;
|
---|
[a9b86d] | 37 |
|
---|
| 38 | MatrixContainer();
|
---|
| 39 | virtual ~MatrixContainer();
|
---|
| 40 |
|
---|
| 41 | bool InitialiseIndices(class MatrixContainer *Matrix = NULL);
|
---|
[9758f7] | 42 | bool ParseMatrix(std::istream &input, int skiplines, int skipcolumns, size_t MatrixNr);
|
---|
[a2d653] | 43 | bool AddMatrix(const std::string &header, const MatrixArray &matrix, size_t MatrixNr);
|
---|
[9758f7] | 44 | virtual bool ParseFragmentMatrix(const std::string name, const std::string prefix, std::string suffix, int skiplines, int skipcolumns);
|
---|
| 45 | bool AllocateMatrix(StringVector GivenHeader, int MCounter, IntVector RCounter, IntVector CCounter);
|
---|
[a9b86d] | 46 | bool ResetMatrix();
|
---|
| 47 | double FindMinValue();
|
---|
| 48 | double FindMaxValue();
|
---|
| 49 | bool SetLastMatrix(double value, int skipcolumns);
|
---|
[9758f7] | 50 | bool SetLastMatrix(const MatrixArray &values, int skipcolumns);
|
---|
[a9b86d] | 51 | //bool ParseIndices();
|
---|
| 52 | //bool SumSubValues();
|
---|
| 53 | bool SumSubManyBodyTerms(class MatrixContainer &Matrix, class KeySetsContainer &KeySets, int Order);
|
---|
[9758f7] | 54 | bool WriteTotalFragments(const std::string name, const std::string prefix);
|
---|
| 55 | bool WriteLastMatrix(const std::string name, const std::string prefix, const std::string suffix);
|
---|
[1b145f] | 56 |
|
---|
| 57 | bool operator==(const MatrixContainer &other) const;
|
---|
| 58 |
|
---|
| 59 | bool operator!=(const MatrixContainer &other) const {
|
---|
| 60 | return !(*this == other);
|
---|
| 61 | }
|
---|
[942906] | 62 |
|
---|
| 63 | private:
|
---|
| 64 | friend class boost::serialization::access;
|
---|
| 65 | // serialization
|
---|
| 66 | template<class Archive>
|
---|
| 67 | void serialize(Archive & ar, const unsigned int version)
|
---|
| 68 | {
|
---|
| 69 | ar & Matrix;
|
---|
| 70 | ar & Indices;
|
---|
| 71 | ar & Header;
|
---|
| 72 | ar & MatrixCounter;
|
---|
| 73 | ar & RowCounter;
|
---|
| 74 | ar & ColumnCounter;
|
---|
| 75 | }
|
---|
[a9b86d] | 76 | };
|
---|
| 77 |
|
---|
| 78 | std::ostream & operator << (std::ostream &ost, const MatrixContainer &m);
|
---|
| 79 |
|
---|
| 80 | #endif /* MATRIXCONTAINER_HPP_ */
|
---|