dbpp_RhsFluxSrc.h

Go to the documentation of this file.

#pragma once
 
// App includes
#include "dbpp_SweRhsAlgorithm.h"
//#include "dbpp_PhysicalBoundaryCnd.h"
#include "dbpp_HLLFluxAlgorithm.h"
#include "dbpp_SrcNumericalTreatment.h"
 
// forward declarations
namespace dbpp {
  class StateVectorField;
  class SrcNumericalTreatment;
  class FiniteVolumeDiscretization;
}
 
namespace dbpp 
{
    class RhsHLLFluxSrc : public SweRhsAlgorithm
    {
  private:
    std::string m_name{"RhsHLLFluxSrc"};          
    bool m_usePressure;                           
    PhyBCNdlConstraint m_leftEnd;                 
    PhyBCNdlConstraint m_rightEnd;                
    const PhysicalSystem* m_phySystem;            
    PhysicalBoundaryCnd m_phyBC;                  
    FluxAlgorithm* m_fluxAlgo{nullptr};           
    SrcNumericalTreatment* m_srcDiscr{ nullptr }; 
    public:
    RhsHLLFluxSrc(std::string aName, const PhysicalSystem* aPhysys, bool aUsePressure = false);
    RhsHLLFluxSrc( std::string aName, HLLFluxAlgorithm* aFluxAlgo, SrcNumericalTreatment* aSrcDiscr,
      const PhysicalBoundaryCnd& aPhybc);
    void setPhysicalBoundaryCnd( const PhysicalBoundaryCnd& aPhyBnd) override final
    {  
      m_phyBC = aPhyBnd;
      m_leftEnd = aPhyBnd.getLeftEnd();
      m_rightEnd = aPhyBnd.getRightEnd();
    }

    PhysicalBoundaryCnd getPhysicalBoundaryCnd() const override final 
    {
      return m_phyBC;
    }

    void setFluxAlgorithm( FluxAlgorithm* aFluxAlgo, 
      bool useIncompleteFlux = true) override 
    {
      m_fluxAlgo = aFluxAlgo;
    }

    void setSourceTermDiscr( SrcNumericalTreatment* aSrcDiscr,
      bool useManningFormula = true) override 
    {
      m_srcDiscr = aSrcDiscr;
    }

    bool useReconstruction() const noexcept override final { return m_fluxAlgo->useReconstr(); }
    bool isFrictionLess()    const noexcept override final { return m_srcDiscr->useFriction(); }
    bool useSourceTerms()    const noexcept override final { return m_srcDiscr!=nullptr; }
    SweRhsData<> calculate( const Sfx::StateVectorField& aU, // for now will do 
      const std::shared_ptr<FiniteVolumeDiscretization>& aGblDiscr) override; //override final
    SweRhsData<> calculate( const Sfx::StateVectorField& aU) override final;
    };
} // End of namespace