governingequationsolver.hpp

#pragma once
 
#include <armadillo>
#include <memory>
 
class Pipeline;
class BoundaryConditions;
class InternalEnergyDiscretizer;
class MatrixEquation;
class Discretizer;
 
class GoverningEquationSolverBase
{
public:
    virtual ~GoverningEquationSolverBase();
 
    virtual arma::mat solve(
            const arma::uword dt,
            const Pipeline& currentState,
            const Pipeline& newState,
            const BoundaryConditions& boundaryConditions) = 0;
 
    bool isOverDetermined(const BoundaryConditions& boundaryConditions);
};
 
template<typename T = InternalEnergyDiscretizer>
class GoverningEquationSolver : public GoverningEquationSolverBase
{
public:
    virtual ~GoverningEquationSolver();
 
    explicit GoverningEquationSolver(const arma::uword nGridPoints);
 
    explicit GoverningEquationSolver(std::unique_ptr<Discretizer> discretizer);
 
    virtual arma::mat solve(
            const arma::uword dt,
            const Pipeline& currentState,
            const Pipeline& newState,
            const BoundaryConditions& boundaryConditions) override;
 
private:
    const arma::uword m_nVariables = 3; 
 
    std::unique_ptr<Discretizer> m_discretizer;
 
    std::unique_ptr<MatrixEquation> m_matrixEquation;
};