fem/pardiso_solve.hh

Go to the documentation of this file.

/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
/*                                                                           */
/*  This file is part of the library KASKADE 7                               */
/*  https://www.zib.de/research/projects/kaskade7-finite-element-toolbox     */
/*                                                                           */
/*  Copyright (C) 2002-2011 Zuse Institute Berlin                            */
/*                                                                           */
/*  KASKADE 7 is distributed under the terms of the ZIB Academic License.    */
/*    see $KASKADE/academic.txt                                              */
/*                                                                           */
/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
 
#ifndef PARDISO_SOLVE_HH
#define PARDISO_SOLVE_HH
 
 
#include <vector>
#include <memory>
class PARDISOFactorization
{
public:
 
  enum MatrixStructure {
    STRUCTURALLY_SYMMETRIC = 1,
    SYMMETRIC_POSITIVE_DEFINITE = 2,
    SYMMETRIC_INDEFINITE = -2,
    NONSYMMETRIC = 11 
  };
  
 
 
  PARDISOFactorization(int n,
                       int mtype_,
                       std::vector<int> const& ridx,
                       std::vector<int> const& cidx,
                       std::vector<double> const& values);
                                     
  ~PARDISOFactorization();
 
  void solve(std::vector<double> const& b,
             std::vector<double>& x) const;
private:
  mutable int N; 
  mutable int n; 
  mutable std::vector<int> Ap, Ai;
  mutable std::vector<double> Az;
  mutable int phase;
 
  mutable int maxfct, mnum, mtype, nrhs, msglvl, error,  idum;
  mutable double ddum;
  static bool first;
  mutable std::vector<int> iparm;
  mutable std::vector<long int> pt;
};
 
 
#endif