PETScSolParams.h

Go to the documentation of this file.

// $Id$
//==============================================================================
//==============================================================================
 
#ifndef _PETSC_SOL_PARAMS_H
#define _PETSC_SOL_PARAMS_H
 
#include "LinSolParams.h"
#include "PETScSupport.h"
 
#include <set>
#include <string>
#include <vector>
 
class ProcessAdm;
class SettingMap;
 
typedef std::vector<int>         IntVec;       
typedef std::vector<IntVec>      IntMat;       
typedef std::vector<std::string> StringVec;    
typedef std::vector<StringVec>   StringMat;    
typedef std::vector<IS>          ISVec;        
typedef std::vector<ISVec>       ISMat;        
 
enum SchurPrec { SIMPLE, MSIMPLER, PCD };
 
 
class PETScSolParams
{
public:
  PETScSolParams(const LinSolParams& spar, const ProcessAdm& padm) :
    params(spar), adm(padm)
  {}
 
  ~PETScSolParams()
  {
  }
 
  void setupPC(PC& pc, size_t block,
               const std::string& prefix,
               const std::set<int>& blockEqs,
               bool setup);
 
  size_t getNoBlocks() const { return params.getNoBlocks(); }
 
  const LinSolParams::BlockParams& getBlock(size_t i) const { return params.getBlock(i); }
 
  int getIntValue(const std::string& key) const { return params.getIntValue(key); }
 
  double getDoubleValue(const std::string& key) const { return params.getDoubleValue(key); }
 
  std::string getStringValue(const std::string& key) const { return params.getStringValue(key); }
 
  bool hasValue(const std::string& key) const { return params.hasValue(key); }
 
  LinAlg::LinearSystemType getLinSysType() const { return params.getLinSysType(); }
 
protected:
  bool addDirSmoother(PC pc, const Mat& P,
                      int iBlock, const ISMat& dirIndexSet);
 
  void setMLOptions(const std::string& prefix, const SettingMap& map);
 
  void setGAMGOptions(const std::string& prefix, const SettingMap& map);
 
  void setHypreOptions(const std::string& prefix, const SettingMap& map);
 
  void setupCoarseSolver(PC& pc, const std::string& prefix, const SettingMap& map);
 
  void setupSmoothers(PC& pc, size_t iBlock,
                      const ISMat& dirIndexSet,
                      const std::set<int>& blockEqs,
                      bool setup);
 
  void setupAdditiveSchwarz(PC& pc, size_t block,
                            bool asmlu, bool smoother,
                            const std::set<int>& blockEqs,
                            bool setup);
 
  const LinSolParams& params; 
  const ProcessAdm& adm;      
};
 
#endif