BiCGSTABSolver.h

Go to the documentation of this file.

#ifndef included_AMP_BiCGSTABSolver
#define included_AMP_BiCGSTABSolver
 
#include "AMP/solvers/SolverStrategy.h"
#include "AMP/solvers/SolverStrategyParameters.h"
#include "AMP/utils/AMP_MPI.h"
 
 
namespace AMP::Solver {
 
 
template<typename T = double>
class BiCGSTABSolver : public SolverStrategy
{
public:
    BiCGSTABSolver();
 
    explicit BiCGSTABSolver( std::shared_ptr<SolverStrategyParameters> params );
 
    static std::unique_ptr<SolverStrategy>
    createSolver( std::shared_ptr<SolverStrategyParameters> params )
    {
        return std::make_unique<BiCGSTABSolver<T>>( params );
    }
 
    virtual ~BiCGSTABSolver();
 
    std::string type() const override { return "BiCGSTABSolver"; }
 
    void apply( std::shared_ptr<const AMP::LinearAlgebra::Vector> f,
                std::shared_ptr<AMP::LinearAlgebra::Vector> u ) override;
 
    void initialize( std::shared_ptr<const SolverStrategyParameters> params ) override;
 
    void registerOperator( std::shared_ptr<AMP::Operator::Operator> op ) override;
 
protected:
    void getFromInput( std::shared_ptr<AMP::Database> db );
    void allocateScratchVectors( std::shared_ptr<const AMP::LinearAlgebra::Vector> u );
 
private:
    int d_restarts = 0; 
 
    bool d_bUsesPreconditioner = false;
 
    std::shared_ptr<AMP::LinearAlgebra::Vector> d_r;
    std::shared_ptr<AMP::LinearAlgebra::Vector> d_r_tilde;
    std::shared_ptr<AMP::LinearAlgebra::Vector> d_p;
    std::shared_ptr<AMP::LinearAlgebra::Vector> d_p_hat;
    std::shared_ptr<AMP::LinearAlgebra::Vector> d_s;
    std::shared_ptr<AMP::LinearAlgebra::Vector> d_s_hat;
    std::shared_ptr<AMP::LinearAlgebra::Vector> d_t;
    std::shared_ptr<AMP::LinearAlgebra::Vector> d_v;
};
} // namespace AMP::Solver
 
#endif