Coulomb_NN_Storage.h

#ifndef QNANO_NEW_COULOMB_NN_STORAGE_DEFINED_H
#define QNANO_NEW_COULOMB_NN_STORAGE_DEFINED_H

#include "tools/Tensor4.h"
#include "tools/Reader.h"
#include "structure/Bond_Type.h"
#include "orbitals/Select_Orbitals_List.h"

class Coulomb_NN_Storage_Entry: public Tensor4, public Bond_Type{
public: 
  
  //Constructors and copy
  void copy(const Coulomb_NN_Storage_Entry &other){
    Tensor4::copy(other);
    Bond_Type::copy(other);
    //orbs=other.orbs;
  }
  Coulomb_NN_Storage_Entry &operator=(const Coulomb_NN_Storage_Entry &other){
    copy(other);
    return *this;
  }
  Coulomb_NN_Storage_Entry(const Coulomb_NN_Storage_Entry &other){
    copy(other);
  }
  Coulomb_NN_Storage_Entry(const Lattice *l, int el0, int el1, const Vec3d &v, const Select_Orbitals_List &orblist):Bond_Type(l,el0,el1,v), Tensor4(){
    int orb0=orblist.find_check(el0);
    int orb1=orblist.find_check(el1);
    Tensor4::resize(orb0,orb1,orb1,orb0);
  }
  Coulomb_NN_Storage_Entry():Tensor4(){
  }
};
class Coulomb_NN_Storage: public List_Class<Coulomb_NN_Storage_Entry>{
public:

  void read(const Lattice &lat, const Select_Orbitals_List &orblist, const std::string &filename);  
  int find(int element0, int element1, const Vec3d &vec, double tolerance=0.05)const;

  Coulomb_NN_Storage();
  Coulomb_NN_Storage(const Lattice &lat, const Select_Orbitals_List &orblist, const std::string &filename){
    read(lat, orblist, filename);
  }
};

#endif