SLEPc_Solver.h

#ifndef QNANO_NEW_SLEPC_SOLVER_DEFINED_H_
#define QNANO_NEW_SLEPC_SOLVER_DEFINED_H_


#include <sys/time.h>
#include <slepceps.h>
#include "MPI/MPI_Communicator.h"
#include "MPI/Distribute_Block_MPI.h"
#include "tools/Sparse_Matrix_Generator.h"
#include "tools/Abstract_Solver.h"

class SLEPc_Solver: public Abstract_Solver{
public:

  Sparse_Matrix_Generator *sparse_generator;

  Mat A;
  EPS eps;
  PetscInt DIM;
  Distribute_Block_MPI db;
  int mlocal, mstart;

  int MPI_rank, MPI_size;
  struct timeval timestart, timeend;
  bool was_set_up;


  PetscErrorCode solve_ierr();
  virtual void solve();


  PetscErrorCode assemble_SetValue();
  PetscErrorCode set_A_size_from_generator();
  PetscErrorCode set_A_from_generator();
  
  PetscErrorCode destroy();
  PetscErrorCode setup_ierr();
  virtual void setup();
  virtual bool should_print()const{
    return MPI_rank==0;
  }
  virtual Distribute_Block_Local* new_Distribute_Block_Local(size_t totalsize)const{
    return new Distribute_Block_MPI(MPI_COMM_WORLD,totalsize);
  }
  virtual Communicator* new_Communicator()const{
    return new MPI_Communicator();
  }


//  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);
//  virtual std::complex<double> get_local_EV_coeff(size_t i, size_t j)const;



  SLEPc_Solver(Sparse_Matrix_Generator *sparse_gen)   
    : sparse_generator(sparse_gen), DIM(0), was_set_up(false){
  }
  ~SLEPc_Solver(){
    if(was_set_up) destroy();
  }
 
};

#endif