ColumnSolver.h

Go to the documentation of this file.

#ifndef included_AMP_ColumnSolver
#define included_AMP_ColumnSolver
 
#include "AMP/solvers/SolverStrategy.h"
#include <string>
#include <vector>
 
namespace AMP::Solver {
 
typedef SolverStrategyParameters ColumnSolverParameters;
 
class ColumnSolver : public SolverStrategy
{
public:
    ColumnSolver() { d_resetColumnOperator = false; }
 
    explicit ColumnSolver( std::shared_ptr<const SolverStrategyParameters> parameters );
 
    static std::unique_ptr<SolverStrategy>
    createSolver( std::shared_ptr<SolverStrategyParameters> solverStrategyParameters )
    {
        return std::make_unique<ColumnSolver>( solverStrategyParameters );
    }
 
    virtual ~ColumnSolver() {}
 
    std::string type() const override { return "ColumnSolver"; }
 
    void apply( std::shared_ptr<const AMP::LinearAlgebra::Vector> f,
                std::shared_ptr<AMP::LinearAlgebra::Vector> u ) override;
 
    void setInitialGuess( std::shared_ptr<AMP::LinearAlgebra::Vector> initialGuess ) override;
 
    void append( std::shared_ptr<AMP::Solver::SolverStrategy> solver );
 
    std::shared_ptr<AMP::Solver::SolverStrategy> getSolver( const int i ) { return d_Solvers[i]; }
 
    int getNumberOfSolvers() { return d_Solvers.size(); }
 
    void resetOperator( std::shared_ptr<const AMP::Operator::OperatorParameters> params ) override;
 
protected:
    void GaussSeidel( std::shared_ptr<const AMP::LinearAlgebra::Vector> &f,
                      std::shared_ptr<AMP::LinearAlgebra::Vector> &u );
 
    void SymmetricGaussSeidel( std::shared_ptr<const AMP::LinearAlgebra::Vector> &f,
                               std::shared_ptr<AMP::LinearAlgebra::Vector> &u );
 
    std::vector<std::shared_ptr<AMP::Solver::SolverStrategy>> d_Solvers;
    std::string d_IterationType;
 
    bool d_resetColumnOperator;
 
private:
};
} // namespace AMP::Solver
 
#endif