dbpp_ODESolver1D.h

Go to the documentation of this file.

#pragma once
 
// app includes
#include "dbpp_NujicIntegrator.h"
#include "../include/DamBreakProb.h"
 
// forward declaration
namespace Sfx { class LDeltaOperator; }
 
// forward declaration
namespace dbpp {
  class PhysicalSystem;
  class FiniteVolumeDiscretization;
}
 
namespace dbpp 
{
  class ODESolver1D 
  {
    using pairvarray = std::pair<std::valarray<float64>, std::valarray<float64>>;
  public:
    ODESolver1D() 
    : m_dbProb{1.,10.}, // default value 
      m_integrator{}, 
      m_physSys{ nullptr },
      m_LdeltaOp{nullptr}
    {}

    ODESolver1D( DamBreakProblem aDbProb) 
    : m_dbProb{ aDbProb }, m_integrator{}, m_physSys{ nullptr }, m_LdeltaOp{nullptr} {}

    ODESolver1D( Sfx::LDeltaOperator* aLdeltaOp) 
    : m_dbProb{0.,0.}, m_integrator{}, m_physSys{ nullptr }, m_LdeltaOp{ aLdeltaOp } {}

    void solve( const std::shared_ptr<dbpp::FiniteVolumeDiscretization>& aFVdiscr, float64 aTime);
    void solve( const Sfx::StateVectorField& aUfield, float64 aTime);
    pairvarray getSolution() const { return m_integrator.getSolutionAsVarray(); } // shall return what?
    void registerPhysicalSystem( PhysicalSystem* aPhysys) noexcept { m_physSys = aPhysys; }
    PhysicalSystem* getPhysicalSystem() const noexcept { return m_physSys;}
  private:
    DamBreakProblem m_dbProb;        
    NujicIntegrator m_integrator;    
    PhysicalSystem* m_physSys;       
    Sfx::LDeltaOperator* m_LdeltaOp; 
  };
} // End of namespace