pardiso.hh

Go to the documentation of this file.

 
#ifndef PARDISO_HH
#define PARDISO_HH
 
#include <cstdlib>
#include "basics/typedefs.hh"
#include "basics/fortranException.hh"
#include "compositions.hh"
 
namespace concepts {
 
  // forwards
  template<typename F>
  class SparseMatrix;
 
  // *************************************************************** Pardiso **
 
  class Pardiso : public Operator<Real> {
  public:
    enum MATRIX_TYPE {STRUCT_SYMM = 1, SPD = 2, SYMM_INDEF = -2, UNSYMM = 11};
    Pardiso(const SparseMatrix<Real>& matrix, enum MATRIX_TYPE type);
    virtual ~Pardiso();
    virtual void operator()(const Function<Real>& fncY, Function<Real>& fncX);
    void operator()(const Vector<Real>& fncY, Vector<Real>& fncX);
    virtual const Space<Real>& spaceX() const { return spcX_; }
    virtual const Space<Real>& spaceY() const { return spcY_; }
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    const uint n_;
    const Space<Real>& spcX_;
    const Space<Real>& spcY_;
    bool symmetric_;
    bool factored_;
 
    uint nnz_;
    void* A_;
    Array<Real> a_;
    Array<uint> ia_, ja_;
    Array<size_t> handle;
    Array<int> iparam;
    int max_factorizations, matrix_number, matrix_type, msglvl, error,
      nb, max_threads;
 
    int solver_pardiso_factor_();
    int solver_pardiso_solve_(double*, double*);
    void errorcheck_(int error) const throw(FortranException);
  };
 
} // namespace concepts
 
#endif // PARDISO_HH