PETScSchurPC.h

Go to the documentation of this file.

// $Id$
//==============================================================================
//==============================================================================
 
#include "PETScMatrix.h"
#include <memory>
 
 
class PETScSchurPC {
public:
  PETScSchurPC(PC& pc_init, const std::vector<Mat>& blocks,
               const LinSolParams::BlockParams& params, const ProcessAdm& adm);
 
  ~PETScSchurPC();
 
  static PetscErrorCode Apply_Schur(Mat A, Vec x, Vec y);
 
  static PetscErrorCode Apply_Outer(PC pc, Vec x, Vec y);
 
  static PetscErrorCode Destroy(PC pc);
 
protected:
  KSP inner_ksp; 
  KSP outer_ksp; 
  Mat outer_mat; 
  Vec tmp; 
  const std::vector<Mat>* m_blocks; 
};