SAMpatchPETSc.h

Go to the documentation of this file.

// $Id$
//==============================================================================
//==============================================================================
 
#ifndef _SAM_PATCH_PETSC_H
#define _SAM_PATCH_PETSC_H
 
#include "SAMpatch.h"
#include "PETScSupport.h"
#include <map>
 
class ProcessAdm;
 
 
class SAMpatchPETSc : public SAMpatch
{
public:
  SAMpatchPETSc(const std::map<int,int>& g2ln, const ProcessAdm& padm);
  virtual ~SAMpatchPETSc();
 
  virtual Real dot(const Vector& x, const Vector& y, char dofType = 'D') const;
 
  virtual Real normL2(const Vector& x, char dofType = 'D') const;
 
  virtual Real normInf(const Vector& x, size_t& comp, char dofType = 'D') const;
 
  virtual bool expandSolution(const SystemVector& solVec,
                              Vector& dofVec, Real scaleSD) const;
 
private:
  // Parameters for parallel computing
  int    nProc;      
  IntVec ghostNodes; 
  IntVec l2gn;       
 
  const ProcessAdm& adm; 
 
#ifdef HAVE_MPI
  mutable IS glob2LocEq = nullptr; 
 
  struct DofIS {
    IS local;        
    IS global;       
    int nDofs;       
    VecScatter ctx;  
    bool scatterCreated = false; 
  };
  mutable std::map<char,DofIS> dofIS; 
 
  DofIS& getIS(char dofType) const;
#endif
};
 
#endif