vector_epa.h

#ifndef VECTOR_EPA_H
#define VECTOR_EPA_H
 
#include "Matrix.h"
#include "PMatrix.h"
#include "active_bandstructure.h"
#include "context.h"
#include "eigen.h"
 
class VectorEPA {
 public:
   VectorEPA(StatisticsSweep &statisticsSweep_, const int &numStates_,
                const int &dimensionality_ = 3);
 
   VectorEPA(const VectorEPA &that);
 
   VectorEPA &operator=(const VectorEPA &that);
 
  virtual Eigen::MatrixXd dot(const VectorEPA &that);
 
  VectorEPA operator*(VectorEPA &that);
 
  VectorEPA operator*(const double &scalar);
 
  VectorEPA operator*(const Eigen::VectorXd &vector);
 
  VectorEPA operator*(ParallelMatrix<double> &matrix);
 
  VectorEPA operator+(VectorEPA &that);
 
  VectorEPA operator-(VectorEPA &that);
 
  VectorEPA operator-();
 
  VectorEPA operator/(VectorEPA &that);
 
  VectorEPA sqrt();
  VectorEPA reciprocal();
 
  void setConst(const double &constant);
 
  double &operator()(const int &iCalc, const int &iDim, const int &iState);
 
  const double &operator()(const int &iCalc, const int &iDim,
                           const int &iState) const;
 
  Eigen::MatrixXd data;
 
  // we store auxiliary objects and parameters
  StatisticsSweep &statisticsSweep;
  int numCalculations;
  int numStates;
  int numChemPots;
  int numTemps;
  int dimensionality;
 
  int glob2Loc(const ChemPotIndex &imu, const TempIndex &it,
                const CartIndex &iDim) const;
  std::tuple<ChemPotIndex, TempIndex, CartIndex> loc2Glob(const int &i) const;
 
  std::vector<int> excludeIndices;
 
 protected:
   VectorEPA baseOperator(VectorEPA &that, const int &operatorType);
  const int operatorSums = 0;
  const int operatorDivs = 1;
  const int operatorProd = 2;
  const int operatorDiff = 3;
};
 
#endif