dbpp_HLLPhysicalFlux.h

Go to the documentation of this file.

 
#pragma once
 
// ...
#include "../SfxBase/StVenant1DEquations.h"
 
namespace dbpp 
{
  class HLLPhysicalFlux
  {             
  public:
    HLLPhysicalFlux( Sfx::SCLEquation* aMthEq) {}
    Sfx::SCLEquation* getSclEquations() { return nullptr; }
    void setSclEquations( const Sfx::SCLEquation* aSclEq) {}
    // Maybe use a factory Dsn::createSclEquations() 
    // File Dsn_SclEquations i think should be more appropriate
    cellFaceFlux calculFF( const Sfx::cellFaceVariables& aULR/*, Sfx::SCLEquation* aScleq*/) const
    {
      // DESIGN NOTE
      //  may be use a factory (why pass it as an arg?) unique_ptr<>  return by factory
      // see Scott Meyers book "Modern Effective C++"
      // new Dsn::StVenant1DEq( std::string("St-Venant1D FrictionLess with Convective flux"))
 
      Sfx::StVenant1D w_Scleq; // incomplete as default;
      
      Sfx::StateVector w_UL = aULR.UL(); // left state
      Sfx::StateVector w_UR = aULR.UR(); // right state
 
      const auto w_FL1 = w_UL.Q(); // flux left state
      const auto w_FR1 = w_UR.Q(); // flux right state
 
      // create a math equation from the factory?? avoid passing math eq as arg
      // "shr_ptr createSCLEquation(incomplete)" since we need it  
      // to evaluate (physical)flux of state variable (node value)  
 
      // compute flux according to math equation
      const auto w_FL2 = w_Scleq.flux( Sfx::StateVector(w_UL.A(), w_UL.Q())); //FL2
      const auto w_FR2 = w_Scleq.flux( Sfx::StateVector(w_UR.A(), w_UR.Q())); //FR2
 
      return { w_FL1, w_FL2, w_FR1, w_FR2 };
    }
  };
}//End of namespace