dbpp_FiniteDifferenceModel.h

Go to the documentation of this file.

 
#pragma once
 
// C++ include
#include <string>
// SfxBase19 includes
#include "include/Sfx_TimePrm.h"
#include "include/Sfx_DefineTypes.h"
// SfxBase19 include
#include "include/Sfx_Directory.h"
// App includes
#include "dbpp_NujicIntegrator.h"
#include "dbpp_PhysicalSystem.h"
#include "dbpp_FiniteVolumeDiscretization.h"
 
// FiniteDifference Model hold a model with finite difference discretization
 
namespace dbpp 
{
    class SweRhsAlgorithm;
  class NujicIntegrator;
//  class GodunovTypeScheme;
    class FiniteVolumeDiscretization;
}
 
namespace dbpp 
{
    class FiniteDifferenceModel 
    {
        using Time = double;
 
      public:
          FiniteDifferenceModel( std::string aName, SweRhsAlgorithm* aRhsAlgo)
            :// m_twoStepsIntegrator{},
        //      m_rhsAlgo{aRhsAlgo},
                m_U12{nullptr,nullptr}
            {}

          void advance( const std::shared_ptr<FiniteVolumeDiscretization>& aGblDiscr, Time aFrom,
                          Time aTo, float64 aNts /*aStepCondition*/) const 
            {           
                // Time parameters (start, stop and time step) read from simulation property file
                Sfx::SfxTimePrm w_tip; // { 0., 1., 0.01);
                using enum Sfx::SfxTimePrm::eTimeStepMode;
                w_tip.setStartTime(aFrom);
                w_tip.setStopTime(aTo); //hudson (0.01) or mcneil (22.5)
                w_tip.setTimeStepMode(VARIABLE_TIME_STEP);
                w_tip.initTimeLoop();
 
                // ... 
                auto w_numrep = aGblDiscr->method();
                w_numrep->isGodunovType();
 
                w_numrep->mainLoop(aGblDiscr, w_tip);
 
            //  const auto& w_Uvalues = aGblDiscr->Uh(); // values()
            //  Sfx::scalarField w_U1; // shall be shared pointer!!
            //  Sfx::scalarField w_U2;
            //      w_tip.increaseTime();
            }

        //  void setInitialState( const PhysicalSystem aSystm) {}
//  private:
        // Concept of an evolver to step through time
    //  NujicIntegrator m_twoStepsIntegrator; /**< explicit integrator*/
    //  StateVector m_U12;                    /**< scalar field over a fd grid*/
    };
}  // End of namespace