dbpp_RK2HLLScheme.h

Go to the documentation of this file.

#pragma once
 
// stl include
#include <vector>
// project includes
#include "dbpp_GodunovTypeScheme.h"
#include "dbpp_NujicIntegrator.h"
 
namespace dbpp { 
    class PhysicalSystem; 
    class SweRhsAlgorithm; 
}
 
namespace dbpp
{
    class Rk2HllScheme : public dbpp::GodunovTypeScheme 
    {
    public:
        Rk2HllScheme( /*SweRhsAlgorithm const std::shared_ptr<ListSectFlow>& aListSectShrd*/);  //something like this!!
        Rk2HllScheme(dbpp::PhysicalSystem* aPhysys/*, dbpp::SweRhsAlgorithm* aSweRhsAlgo*/);
        ~Rk2HllScheme();
 
        // Validating EMcNeilPhysicalAlgorithm (need to be refactored!!!)
        // Design Note (next version)
        //  Main purpose of this method its let user implement its own algorithm
        //  similar to EMcNeilAlgorithm::calculate() and not to set pointer
        //  setting pointer is dangereous (who responsible ot delete the pointer?)
        //  result od dangling pointer. In this case pas the algorithm in ctor 
        //  Also, maybe pass a unique_ptr<SweRhsAlgorithm> as argument
        //  could provide another service witn no argument pass to method
        //  void setSweRhsAlgorithm() this method let user implement the algo
        void setSweRhsAlgorithm( SweRhsAlgorithm* aRhsAlgo) override
        { 
            // Design Note
            // NO!! see SemiDiscreteMethod note. this is wrong
            m_rhsAlgo = aRhsAlgo; 
        }

        SweRhsAlgorithm* getSweRhsAlgorithm() const override final 
        { 
            return m_rhsAlgo; 
        }

        SweRhsAlgorithm* getSweRhsAlgorithm() override final
        {
            return m_rhsAlgo;
        }

        void setSourceTerm( dbpp::SrcNumericalTreatment* aSrcTreatment) override final {}
    void setFaceFluxAlgorithm(/*FluxAlgorithm*/) override final {}
  //  void setApproximateRiemannSolver( emcil::ApprRiemannSolver* ARprob) override final {}
        /*
        *  @brief update nodal values
        *  @return global nodal values
        */
        Uh updatedValues() override final;
 
    protected:
        void advance(/*dbpp::PhysicalSystem* aPhysys,*/ 
            const std::shared_ptr<dbpp::FiniteVolumeDiscretization>& aGblDiscr, const double aDt) override final;
    
        //void setPhysicalBoundaryCnd( const PhysicalBoundaryCnd& aPhyBCond) 
        //{
        //  m_phyBC =  aBCond.getBCpair();
        //}
 
  private:
    NujicIntegrator m_twoStepsIntegrator; 
      dbpp::PhysicalSystem* m_physys;       
      dbpp::SweRhsAlgorithm* m_rhsAlgo;     
        FluxTensorMap computeFaceFlux( const dbpp::Omega& aOmega, 
            const std::vector<float64>& aVecA, const std::vector<float64>& aVecQ);
        FluxTensorMap computeFaceFlux( const dbpp::Omega& aOmega,
            std::vector<float64>&& aVecA, std::vector<float64>&& aVecQ);
    };
} // End of namespace