Coulomb_Matrix_Elements.h

#ifndef QNANO_NEW_COULOMB_MATRIX_ELEMENTS_DEFINED_H
#define QNANO_NEW_COULOMB_MATRIX_ELEMENTS_DEFINED_H

#include "tools/Tensor4.h"
#include "tools/QNANO_Matrix.h"
#include "tools/Reader.h"
#include "manybody/Manybody_Hilbert_Space.h"

class Coulomb_Matrix_Elements: public Tensor4{
public: 
  enum Type{ none, direct, exchange, ee, hh} type;
  static std::string get_type_string(enum Type type_);

  void resize(int nr_e, int nr_h);
  void resize(const Manybody_Hilbert_Space &HS);
  void setup(Type type_, const Manybody_Hilbert_Space &HS, const std::string &filename=""){
    type=type_;
    resize(HS);
    if(filename!="")read(filename);
  }

  void check_size_consistency()const;
  void check_Hermitian(double accuracy=1e-6)const;
  void check_consistency(double accuracy=1e-6)const{check_size_consistency(); check_Hermitian(accuracy);}
  
  bool is_hole(size_t i)const;

  void read(const std::string &filename, bool read_size_from_file=false);
  void read_noheader(std::ifstream &ifs, const std::string &filename_for_error="");
  void print(const std::string &filename, int round_figure=8)const;
 

  void apply_trafo(const QNANO_Matrix &trafo_e, const QNANO_Matrix &trafo_h); 

  void make_Hermitian();

  //Copy and Constructors
  void copy_noresize(const Coulomb_Matrix_Elements& other);

  void copy(const Coulomb_Matrix_Elements& other){
    type=other.type;
    Tensor4::copy(other);
  }
  Coulomb_Matrix_Elements &operator=(const Coulomb_Matrix_Elements &other){
    copy(other);
    return *this;
  }
  Coulomb_Matrix_Elements(const Coulomb_Matrix_Elements &other){
    copy(other);
  }
  Coulomb_Matrix_Elements(Type type_=none) : type(type_) {}
  Coulomb_Matrix_Elements(Type type_, size_t i, size_t j, size_t k, size_t l)
    : Tensor4(i,j,k,l), type(type_) {
//    check_consistency();
  }
  Coulomb_Matrix_Elements(Type type_, int nr_e, int nr_h) : type(type_){
    resize(nr_e,nr_h);
    check_consistency(); 
  }
  Coulomb_Matrix_Elements(Type type_, const Manybody_Hilbert_Space &HS, const std::string &filename="") : type(type_){
    resize(HS);
    if(filename!="")read(filename);
  }
  Coulomb_Matrix_Elements(Type type_, const std::string &filename)
    : type(type_) {
    read(filename, true);
    check_consistency();
  }



};


#endif