PETScLinearSolver.h

Go to the documentation of this file.

 
#pragma once
 
#include <petscksp.h>
 
#include <string>
 
#include "PETScMatrix.h"
#include "PETScVector.h"
 
namespace MathLib
{
class PETScLinearSolver
{
public:
    PETScLinearSolver(std::string const& prefix,
                      std::string const& petsc_options);
 
    ~PETScLinearSolver() { KSPDestroy(&solver_); }
    // TODO check if some args in LinearSolver interface can be made const&.
    bool solve(PETScMatrix& A, PETScVector& b, PETScVector& x);
 
    PetscInt getNumberOfIterations() const
    {
        PetscInt its = 0;
        KSPGetIterationNumber(solver_, &its);
        return its;
    }
 
    double getElapsedTime() const { return elapsed_ctime_; }
 
private:
    KSP solver_;  
    PC pc_;       
 
    double elapsed_ctime_ = 0.0;  
};
 
}  // namespace MathLib