fdm_cg_solver2.h

Go to the documentation of this file.

// Copyright (c) 2018 Doyub Kim
//
// I am making my contributions/submissions to this project solely in my
// personal capacity and am not conveying any rights to any intellectual
// property of any third parties.
 
#ifndef INCLUDE_JET_FDM_CG_SOLVER2_H_
#define INCLUDE_JET_FDM_CG_SOLVER2_H_
 
#include <jet/fdm_linear_system_solver2.h>
 
namespace jet {
 
class FdmCgSolver2 final : public FdmLinearSystemSolver2 {
 public:
    FdmCgSolver2(unsigned int maxNumberOfIterations, double tolerance);
 
    bool solve(FdmLinearSystem2* system) override;
 
    bool solveCompressed(FdmCompressedLinearSystem2* system) override;
 
    unsigned int maxNumberOfIterations() const;
 
    unsigned int lastNumberOfIterations() const;
 
    double tolerance() const;
 
    double lastResidual() const;
 
 private:
    unsigned int _maxNumberOfIterations;
    unsigned int _lastNumberOfIterations;
    double _tolerance;
    double _lastResidual;
 
    // Uncompressed vectors
    FdmVector2 _r;
    FdmVector2 _d;
    FdmVector2 _q;
    FdmVector2 _s;
 
    // Compressed vectors
    VectorND _rComp;
    VectorND _dComp;
    VectorND _qComp;
    VectorND _sComp;
 
    void clearUncompressedVectors();
    void clearCompressedVectors();
};
 
typedef std::shared_ptr<FdmCgSolver2> FdmCgSolver2Ptr;
 
}  // namespace jet
 
#endif  // INCLUDE_JET_FDM_CG_SOLVER2_H_