mumps.hh

Go to the documentation of this file.

 
#ifndef mumps_hh_
#define mumps_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]
#define CNTL(I) cntl[(I)-1]
 
 
namespace concepts {
 
  // **************************************************************** Mumps **
 
  template<class F>
  class Mumps : public VecOperator<F> {
  public:
#ifdef HAS_MPI
    Mumps(const Operator<F>& A, const int rank=0, const int comm=-987654);
    Mumps(const Operator<F>& A, const bool symmetric, const int rank=0, const int comm=-987654);
#else
    Mumps(const Operator<F>& A, const int rank=-1, const int comm=-987654);
    Mumps(const Operator<F>& A, const bool symmetric, const int rank=-1, const int comm=-987654);
#endif
    Mumps(const int rank, const int comm=-987654);
 
    virtual ~Mumps();
 
    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_;
    int* irn_;
    int* jcn_;
    uint n_;
    uint nz_;
    int rank_;
    int comm_;
    void* id_mumps_;
    bool trivial_matrix_;
    bool analysis_;
  };
 
}
 
#endif // mumps_hh