linsolver.h

Go to the documentation of this file.

 
#pragma once
 
#include <memory>
#include <string>
 
#include <petscmat.h>
#include <petscsys.h>
#include <petscvec.h>
#include <yaml-cpp/yaml.h>
 
#include <petibm/type.h>
 
namespace petibm
{
namespace linsolver
{
class LinSolverBase
{
public:
    LinSolverBase() = default;
 
    LinSolverBase(const std::string &solverName, const std::string &file)
        : name(solverName), config(file){};
 
    virtual ~LinSolverBase() = default;
 
    virtual PetscErrorCode destroy();
 
    PetscErrorCode printInfo() const;
 
    PetscErrorCode getType(std::string &_type) const;
 
    virtual PetscErrorCode setMatrix(const Mat &A) = 0;
 
    virtual PetscErrorCode solve(Vec &x, Vec &b) = 0;
 
    virtual PetscErrorCode getIters(PetscInt &iters) = 0;
 
    virtual PetscErrorCode getResidual(PetscReal &res) = 0;
 
protected:
    std::string name;
 
    std::string config;
 
    std::string type;
 
    virtual PetscErrorCode init() = 0;
 
};  // LinSolver
 
}  // end of namespace linsolver
 
namespace type
{
typedef std::shared_ptr<linsolver::LinSolverBase> LinSolver;
 
}  // end of namespace type
 
namespace linsolver
{
PetscErrorCode createLinSolver(const std::string &solverName,
                               const YAML::Node &node, type::LinSolver &solver);
 
}  // end of namespace linsolver
 
}  // end of namespace petibm