dbpp_SemiDiscreteMethod.h

Go to the documentation of this file.

#pragma once
 
// App include
#include "dbpp_NumericalMethod.h"
#include "../SfxBase/Sfx_LDeltaOperator.h"
//#include "dbpp_FiniteVolumeDiscretization.h"
 
// forward declaration
namespace dbpp {
  class ODESolver1D;
  class SweRhsAlgorithm; 
}
 
namespace dbpp
{
  class SemiDiscreteMethod final : public NumericalMethod 
  {
  public:
    SemiDiscreteMethod();
    SemiDiscreteMethod( /*ODESolver1D* aODEsolver,*/ SweRhsAlgorithm* aRhsAlgo);
    ~SemiDiscreteMethod();
    void setSweRhsAlgorithm( SweRhsAlgorithm* aRhsAlgo) override final 
    {
      // Design Note
      //  Shall implement discretization of the RHS (Right-Hand-Side)
      //  For example, use convective flux algorithm and source algo
      //  similar what we do in Physical algorithm.
      //  --- NujicENOFlux
      //  --- JBSourceTerms
      //  --- PressusreTerms
      //  m_rhsAlgo = create SweRhsAlgorithm(...)
      // or may use already defined algorithm available such as
      // RhsFluxSrc or whatever use decide to prototype
 
      m_rhsAlgo = aRhsAlgo; // this is wrong!!! don't make sense!!!
    }

    SweRhsAlgorithm* getSweRhsAlgorithm() const override final { return m_rhsAlgo; }
    SweRhsAlgorithm* getSweRhsAlgorithm() override final { return m_rhsAlgo; }
    bool isSemiDiscreteMethod() const override final { return true; }
    bool isTimeDependent() const override final      { return true; }
    ODESolver1D* getODESolver1D() const { return m_odeSolver1D; }
    void mainLoop(const std::shared_ptr<dbpp::FiniteVolumeDiscretization>& aGlbDiscr,
      Sfx::SfxTimePrm& aTime) override final;
    Uh updatedValues() override final;
  private:
    ODESolver1D* m_odeSolver1D;     
    Sfx::LDeltaOperator m_ldeltaOp; 
    SweRhsAlgorithm* m_rhsAlgo;     
    float64 timeStepUpdate( std::vector<float64>&& aVecA,
      std::vector<float64>&& aVecQ);
  };
} // End of namespace