mumpsoverlap.hh

Go to the documentation of this file.

 
#ifndef mumpsoverlap_hh_
#define mumpsoverlap_hh_
 
#include "compositions.hh"
#include "basics/typedefs.hh"
#ifdef HAS_MPI
#include <mpi.h>
#endif
 
#include "mumps_flags.hh"
 
// Define MUMPS internal constants and functions
#define JOB_INIT -1
#define JOB_END -2
#define ICNTL(I) icntl[(I)-1]
 
 
namespace concepts {
 
  // **************************************************************** MUMPS **
 
  template<class F>
  class MumpsOverlap : public VecOperator<F> {
  public:
    MumpsOverlap(const int NGlob, const Operator<F>& A, const int* LocToGLob, const int rank, const int comm=-987654, const bool CConvention = true);
 
    MumpsOverlap(const Operator<F>& A, const int rank, const int comm=-987654, const bool CConvention = true);
 
    virtual ~MumpsOverlap();
 
    virtual void show_messages();
 
    virtual void hide_messages();
 
    virtual void apply_(const Vector<F>& fncY, Vector<F>& fncX);
 
    virtual void apply_();
  protected:
    virtual std::ostream & info(std::ostream &) const;
  private:
    void* A_loc_;
    int* irn_loc_;
    int* jcn_loc_;
    uint n_;
    uint nz_loc_;
    int rank_;
    int comm_;
    void* id_mumps;
  };
 
}
 
#endif // mumpsoverlap_hh