dbpp_LDeltaOperator.h

Go to the documentation of this file.

#pragma once
 
// C++ include (numeric)
#include <valarray>
// App include
#include "Sfx/Sfx_StateVectorField.h"
 
namespace dbpp { class SweRhsAlgorithm; }
 
namespace dbpp 
{
  using float64 = double; // temp fix
  class SweLDeltaOperator 
  {
  public: 
    SweLDeltaOperator(/*GblDiscr??*/SweRhsAlgorithm* aSweRhsAlgo /*predictor step*/) : m_rhsAlgo{aSweRhsAlgo} {}
    
    // should be virtual pure (interface) and overriden by subclass
    SweRhsEquations applyTo( const std::valarray<float64>& aU1, 
      const std::valarray<float64>& aU2, float64 aDt)
    {
      switch (switch_on)
      {
      default:
        break;
      }
 
      // -------------------------PREDICTOR
      Sfx::StateVectorField w_stateVecField{aName,aU1,aU2};
      SweRhsEquationsData w_rhsEqTerms;
 
      auto w_rhsTerms = m_rhsAlgo->calculate(/*GblDiscr??*/w_stateVecField);
      w_rhsEqTerms.setMassEq(aU1 - lambda * w_rhsTerms.dFF1); // mass equations
      w_rhsEqTerms.setMomentumEq(aU2 - lambda * w_rhsTerms.dFF1 - dt * w_rhsTerms.S); // momentum equation
      
      return w_rhsEqTerms;
      
      // -------------------------CORRECTOR
      auto w_rhsTerms = m_rhsAlgo->calculate(w_stateVecField.asStateVector());
      w_rhsEqTerms.setMassEq(aU1 - aU2 - lambda* w_rhsTerms.dFF1); // mass equation
      w_rhsEqTerms.setMomentumEq(0.5*(aU2 - lambda * w_rhsTerms.dFF1 - dt * w_rhsTerms.S)); // momentum equation
      
      return w_rhsEqTerms;
    }
 
    bool isTimeDependent() const { return true; }
    unsigned getNumberEquations() const { return 2; }
    void setStep(/*predictor/corrector*/) {/*two-time-steps*/ }
  protected:
    // ??? all contains in SweRhsAlgorithm!! isn't?
    virtual void setCalculFF(/*FluxAlgorithm*/)
    {
      aFluxAlgo->isGodunovType();
      m_rhsAlgo->setFaceFluxAlgo(aFFalgo);
    }
    virtual void setSourceTerms(/*SrcTermDiscr*/) {/**/ }
    virtual void setPressureTerm(/**/) {/**/}
  private:
    /*GblDiscr?? m_gblDiscr*/
    SweRhsAlgorithm* m_rhsAlgo;
  };
} // End of namespace