BandedSolver.h

Go to the documentation of this file.

#ifndef included_AMP_BandedSolver
#define included_AMP_BandedSolver
 
#include "AMP/operators/Operator.h"
#include "AMP/solvers/SolverStrategy.h"
#include "AMP/vectors/Vector.h"
 
 
namespace AMP::Solver {
 
class BandedSolver : public SolverStrategy
{
public:
    explicit BandedSolver( std::shared_ptr<SolverStrategyParameters> params );
 
    static std::unique_ptr<SolverStrategy>
    createSolver( std::shared_ptr<SolverStrategyParameters> params )
    {
        return std::make_unique<BandedSolver>( params );
    }
 
    virtual ~BandedSolver();
 
    std::string type() const override { return "BandedSolver"; }
 
    void apply( std::shared_ptr<const AMP::LinearAlgebra::Vector> f,
                std::shared_ptr<AMP::LinearAlgebra::Vector> u ) override;
 
    void resetOperator( std::shared_ptr<const AMP::Operator::OperatorParameters> params ) override;
 
    void reset( std::shared_ptr<SolverStrategyParameters> params ) override;
 
 
protected:
    BandedSolver();
 
 
private:
    int N, M, KL, KU;
    double *AB;
    int *IPIV;
    std::shared_ptr<AMP::Discretization::DOFManager> rightDOF;
    std::shared_ptr<AMP::Discretization::DOFManager> leftDOF;
};
} // namespace AMP::Solver
 
#endif