Eigenvector_Printer_Eigen.h

#ifndef QNANO_NEW_EIGENVECTOR_PRINTER_EIGEN_DEFINED_H
#define QNANO_NEW_EIGENVECTOR_PRINTER_EIGEN_DEFINED_H

#include "tools/Eigenvector_Printer.h"
#include "eigen/Eigen_Diagonalize_Sparse_Matrix.h"


class Eigenvector_Printer_Eigen: public Eigenvector_Printer{
public:
  Eigen::SelfAdjointEigenSolver<Eigen::MatrixXcd> *solver;
  bool use_target;
  double target;
  int eps_nev;

  virtual void setup(Parameter_Map &param_map);


  virtual int get_n_converged() const;
//  virtual int get_actual_nr_print()const;
  virtual double get_eigenvalue_original_order(int i) const;
  virtual double get_eigenvalue(int i) const;
  
  int get_nth_closest_to_target(int n)const;
  virtual void print_eigenvector(int i)const;
  virtual void print_eigenvector_Kramers(int i)const;
//  virtual void print_eigenvalues()const;


  //Kramers part:
  virtual bool has_Kramers_partner(int i)const;
  virtual void calculate_alpha_beta_Kramers_sz( \
    std::complex<double> &alpha, std::complex<double> &beta, int i) const;
  virtual std::complex<double> calculate_max_phase_factor( \
    std::complex<double> alpha, std::complex<double> beta, int i) const;

  //Constructors
  Eigenvector_Printer_Eigen(  \
    Eigen::SelfAdjointEigenSolver<Eigen::MatrixXcd> *esolv, \
    Organize_Blocks *org_bl, bool spinless=true)  \
    : Eigenvector_Printer(org_bl, spinless), solver(esolv){
  }

  Eigenvector_Printer_Eigen(Parameter_Map &param_map, \
    Eigen::SelfAdjointEigenSolver<Eigen::MatrixXcd> *esolv, \
    Organize_Blocks *org_bl, bool spinless=true) \
    : Eigenvector_Printer(org_bl, spinless), solver(esolv){

    setup(param_map);
  }
  
  Eigenvector_Printer_Eigen(Eigen_Diagonalize_Sparse_Matrix *solv, Organize_Blocks *org_bl, bool spinless=true) \
   : Eigenvector_Printer(org_bl, spinless), solver(&solv->solver){
  }
};

#endif