RadiationDiffusionFDOpSplitPrec.h

Go to the documentation of this file.

#ifndef RAD_DIF_OP_SPLIT_PREC
#define RAD_DIF_OP_SPLIT_PREC
 
#include "AMP/IO/AsciiWriter.h"
#include "AMP/IO/PIO.h"
#include "AMP/operators/LinearOperator.h"
#include "AMP/operators/Operator.h"
#include "AMP/operators/OperatorParameters.h"
#include "AMP/operators/radiationDiffusionFD/RadiationDiffusionFDBDFWrappers.h"
#include "AMP/operators/radiationDiffusionFD/RadiationDiffusionFDDiscretization.h"
#include "AMP/solvers/SolverFactory.h"
#include "AMP/solvers/SolverStrategy.h"
#include "AMP/solvers/SolverStrategyParameters.h"
#include "AMP/utils/AMPManager.h"
#include "AMP/vectors/MultiVector.h"
#include "AMP/vectors/Vector.h"
 
#include <iomanip>
#include <iostream>
 
namespace AMP::Solver {
 
 
class BDFRadDifOpPJacOpSplitPrec : public AMP::Solver::SolverStrategy
{
 
    //
public:
    BDFRadDifOpPJacOpSplitPrec( std::shared_ptr<AMP::Solver::SolverStrategyParameters> params );
 
    static std::unique_ptr<AMP::Solver::SolverStrategy>
    create( std::shared_ptr<AMP::Solver::SolverStrategyParameters> params )
    {
        return std::make_unique<BDFRadDifOpPJacOpSplitPrec>( params );
    };
 
    std::string type() const override { return "BDFRadDifOpPJacOpSplitPrec"; };
 
    void apply( std::shared_ptr<const AMP::LinearAlgebra::Vector> bET_,
                std::shared_ptr<AMP::LinearAlgebra::Vector> ET_ ) override;
 
    void reset( std::shared_ptr<AMP::Solver::SolverStrategyParameters> params ) override;
 
 
    //
private:
    static constexpr double SINGULAR_TOL = 1e-15;
 
    std::unique_ptr<AMP::Solver::SolverStrategy> d_difSolverE = nullptr;
    std::unique_ptr<AMP::Solver::SolverStrategy> d_difSolverT = nullptr;
 
 
    void diffusionSolve( std::shared_ptr<const AMP::LinearAlgebra::Vector> bE,
                         std::shared_ptr<const AMP::LinearAlgebra::Vector> bT,
                         std::shared_ptr<AMP::LinearAlgebra::Vector> E,
                         std::shared_ptr<AMP::LinearAlgebra::Vector> T ) const;
 
    void reactionSolve( std::shared_ptr<const AMP::LinearAlgebra::Vector> bE,
                        std::shared_ptr<const AMP::LinearAlgebra::Vector> bT,
                        std::shared_ptr<AMP::LinearAlgebra::Vector> E,
                        std::shared_ptr<AMP::LinearAlgebra::Vector> T );
 
    void scalar2x2Solve(
        double a, double b, double c, double d, double e, double f, double &x, double &y ) const;
 
    void setDiffusionSolvers();
};
 
 
} // namespace AMP::Solver
 
 
#endif