| [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_ */
 | 
|---|