linsolverksp.h

Go to the documentation of this file.

 
#pragma once
 
#include <petscksp.h>
 
#include <petibm/linsolver.h>
 
namespace petibm
{
namespace linsolver
{
class LinSolverKSP : public LinSolverBase
{
public:
    LinSolverKSP(const std::string &solverName, const std::string &file);
 
    virtual ~LinSolverKSP();
 
    virtual PetscErrorCode destroy();
 
    virtual PetscErrorCode setMatrix(const Mat &A);
 
    virtual PetscErrorCode solve(Vec &x, Vec &b);
 
    virtual PetscErrorCode getIters(PetscInt &iters);
 
    virtual PetscErrorCode getResidual(PetscReal &res);
 
protected:
    KSP ksp;
 
    virtual PetscErrorCode init();
 
};  // LinSolverKSP
 
}  // end of namespace linsolver
 
}  // end of namespace petibm