Eigenvector_Printer_SLEPc.h

#ifndef QNANO_NEW_EIGENVECTOR_PRINTER_SLEPC_DEFINED_H
#define QNANO_NEW_EIGENVECTOR_PRINTER_SLEPC_DEFINED_H

#include <slepceps.h>
#include "tools/Eigenvector_Printer.h"
#include "slepc/SLEPc_Solver.h"
//#include "slepc/SLEPc_Kramers_Eigenvector.h"


class Eigenvector_Printer_SLEPc: public Eigenvector_Printer{
public:
  int MPI_rank, MPI_size;

  Mat *A;
  EPS *eps;
  int nrorb;

  Organize_Blocks_By_List local_orb_blocks;
  void setup_local_orb_blocks();

  void setup(Parameter_Map &param_map);
  void setup(){ Parameter_Map param_map; setup(param_map); }

  Eigenvector_Printer_SLEPc(Mat *A_, EPS *eps_, Organize_Blocks *org_bl) : Eigenvector_Printer(org_bl), A(A_), eps(eps_) {
    setup();
  }

  Eigenvector_Printer_SLEPc(SLEPc_Solver *solver, Organize_Blocks *org_bl) : Eigenvector_Printer(org_bl), A(&solver->A), eps(&solver->eps) {
    setup();
  }
/*
  Eigenvector_Printer_SLEPc(Parameter_Map &param_map) {
    A=NULL;
    eps=NULL;
    setup(param_map);
  }
  Eigenvector_Printer_SLEPc(){
    A=NULL; 
    eps=NULL;
    setup();
  }
*/
  virtual ~Eigenvector_Printer_SLEPc(){}

  void check_fully_initialized() const;

  virtual int get_n_converged() const;
  virtual double get_eigenvalue(int i) const;
  
//  PetscErrorCode print_SLEPc(int n_print) const;
  
  virtual void print_eigenvalues() const;
  virtual void print_eigenvector(int i) const;//{ print_SLEPc(i); }
  virtual void print_eigenvector_Kramers(int i)const;


  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;


};



#endif