dbpp_VS19WaveSimulator.h

Go to the documentation of this file.

#pragma once
 
// C++ include
#include <fstream>
// SfxBase19 library includes
#include "include/Sfx_DefineTypes.h"
#include "include/Sfx_WaveFunc.h"
#include "include/Sfx_Simulator.h"
#include "include/Sfx_Observer.h"
// App include
#include "dbpp_ListSectionsFlow.h"
 
// forward declaration
namespace Sfx{ class SfxTimePrm; }
namespace dbpp { class Uh; }
 
namespace dbpp 
{
  class VS19WaveSimulator : public Sfx::Simulator,
                            virtual public Sfx::Observer,
    public std::enable_shared_from_this<VS19WaveSimulator>
  {
  public:
    using ListofSectionsFlow = std::shared_ptr<dbpp::ListSectionsFlow>;

    enum class eSurfaceWaveForm
    {
      SINE = 0, 
      STEP = 1  
    };

    enum class eRiverBedType
    {
      FLAT = 0,   
      STEP = 1,   
      JUMP = 2,   
      RANDOM = 3  
    };
  public:
    VS19WaveSimulator( int32 aNbIterationsMax = 50, float64 aCFL = 0.600000000)
    : m_nbIterationsMax{ aNbIterationsMax },
      m_cfl{aCFL},
      m_saveResult2File{}
    {}

    void run() override final;
    bool init() override final;
    void unInitialize() override final {/*nothing to do*/ }
    VS19WaveSimulator(const VS19WaveSimulator&) = delete;
    VS19WaveSimulator& operator= (const VS19WaveSimulator&) = delete;
    VS19WaveSimulator( VS19WaveSimulator&&) = delete;
    VS19WaveSimulator& operator= (const VS19WaveSimulator&&) = delete;
    void update() override final;
    std::valarray<float64> getH() const noexcept
    {
      return std::valarray<float64>(m_waveDepthProfile.data(), m_waveDepthProfile.size());
    }

    std::shared_ptr<VS19WaveSimulator> sharedPtrTo() { return shared_from_this(); }
  private:
    int32 m_nbIterationsMax;                  
    float64 m_cfl;                           
    std::string m_activeAlgo;                
    std::string m_validationDir;             
    bool m_saveResult2File;                  
    std::ofstream m_fileRes;                 
    std::vector<float64> m_waveDepthProfile; 
    ListofSectionsFlow m_ListSectFlow;       
 
    void setH();  // load H into lambda for efficiency
    void saveResult( const Uh& aUh, const Sfx::SfxTimePrm& aTimeSimEnd);
#if 0 // may be re-use some of them if it make sense
    fieldptr m_lambda;                   
    Sfx::SfxTimePrm m_tip;               
    wavfuncptr m_H;                      
    wavfuncptr m_I;                      
    emcil::ListSectFlow* m_ListSectFlow; 
    double m_CFL;                        
    void setH();                         
    std::string m_activeAlgo;            
    std::string m_validationDir;         
    bool m_saveResult2File;              
    unsigned int m_NumberIterationsMax;  
    double m_Phi1;                       
    double m_Phi0;                       
    double m_finalTime;                  
    eIC m_IC;                            
#endif
  };
} // End of namespace