dbpp_HLLegacyFluxAlgorithm.h

Go to the documentation of this file.

#pragma once
 
// Nov19 DamBreak includes
#include "dbpp_FluxTensor.h"
#include "dbpp_FluxAlgorithm.h"
#include "dbpp_HLLSolver1D.h"
#include "dbpp_RiemannProblem.h"
#include "SimulationConfig.h"     // pointer-2-function
// SfxBase19 include
#include "include/Sfx_cellFaceVariables.h"
 
// forward declarations 
namespace dbpp { 
  class cellFace;
  class PhysicalBoundaryCnd; 
}
 
namespace dbpp 
{
  using valarray64_t = std::valarray<float64>;            
  using ff12pair = std::pair<valarray64_t, valarray64_t>; 

  class Ptr2FLegacyFluxAlgorithm : public FluxAlgorithm
  {
  private:
    std::string m_ptr2fName; 
    int32 m_varOrder;        
    CalculFF m_ptr2func;     
  protected:
    FluxTensorMap calculFF( const std::valarray<float64>& aU1, const std::valarray<float64>& aU2) const;
  public:
    Ptr2FLegacyFluxAlgorithm(std::string aFluxAlgoName, CalculFF aPtr2Func);
    FluxTensor calculFF( const Sfx::cellFaceVariables& aCellFaceVar) override final;
    FluxTensorMap
      calculFF( const Sfx::scalarField1D& U1, const Sfx::scalarField1D& U2,
        const Omega& aDomain, const PhysicalBoundaryCnd& aPhysbc) override final;
    FluxTensor calculFF( const Sfx::StateVector& aUL, const Sfx::StateVector& aUR,
      const cellFace& aCellFace) override final;
    FluxTensorMap calculFF( const Sfx::scalarField1D& U1, const Sfx::scalarField1D& U2,
        const PhysicalBoundaryCnd& aPhysbc) override final;
    // MUSCL as default
    //  eReconstrType getReconstrType() const { return m_reconstrType; }
    int32 getReconstrVarOrder() const noexcept { return m_varOrder; }
    bool useReconstr() const noexcept override { return true; }
    void setReconstrType(eReconstrType aRecnstrType) noexcept override {/*no implemnentation*/ }  //ON/OFF??
    virtual eReconstrType getReconstrType() const noexcept override { return eReconstrType::None; }
    bool usePhysicalInCompleteFlux() const noexcept override { return true; }
   // void setNujicAlphaCoeff(float64 aAlphaCoeff) noexcept { m_alphaCoeff = aAlphaCoeff; }
   // float64 getNujicAlphaCoeff() const noexcept { return m_alphaCoeff; }
    std::string PtrToFuncName() const noexcept { return m_ptr2fName; }
  };
} // End of namespace