Eigen_Diagonalize_Sparse_Matrix.h

#ifndef QNANO_NEW_EIGEN_DIAGONALIZE_SPARSE_MATRIX_DEFINED_H
#define QNANO_NEW_EIGEN_DIAGONALIZE_SPARSE_MATRIX_DEFINED_H

#include "eigen/include_Eigen.h"
#include "tools/Sparse_Matrix_Generator.h"
#include "tools/Abstract_Solver.h"

class Eigen_Diagonalize_Sparse_Matrix: public Abstract_Solver{
public:
  Sparse_Matrix_Generator *generator;

  Eigen::MatrixXcd A;
  Eigen::SelfAdjointEigenSolver<Eigen::MatrixXcd> solver;

  
  void generate_Matrix();
  virtual void setup();
  virtual void solve();
  virtual bool should_print()const{
    return true;
  }
  //virtual size_t get_local_EV_start()const;
  //virtual size_t get_local_EV_blocksize()const;
  virtual double get_eigenvalue(size_t i);
  virtual void get_EV_local_part(size_t i, std::complex<double> *ptr);


  void solve(Sparse_Matrix_Generator *generator);

  void print_eigenvalues(std::ostream &ofs=std::cout)const;
  void print_eigenvalues(const std::string &filename)const;

  Eigen_Diagonalize_Sparse_Matrix() : generator(NULL) { } 
  Eigen_Diagonalize_Sparse_Matrix(Sparse_Matrix_Generator *generator_)
    : generator(generator_) {
  }

  virtual ~Eigen_Diagonalize_Sparse_Matrix(){
  }
};

#endif