Testvs19_ProtoAlgoPOM.h

Go to the documentation of this file.

 
#pragma once
 
// Stl include
#include <vector>
// SfxBase package include
#include "include/Sfx_IPhysicalAlgorithm.h"
// SweNumeric package includes
#include "../../SfxBase/Sfx_ImposeBnd.h"
#include "../../SfxBase/Sfx_LDeltaOperator.h"
 
namespace Sfx {class PhysicalSystem; }
namespace Testvs19 { class RetRhsType; }
 
namespace Testvs19 
{
    // Actually, numerical integration of the ODE U_t=L(U).
    // We use a Runge-Kutta solver family to solve this problem.
    class EmcilAlgorithm : public Sfx::IPhysicalAlgorithm
    {
    public:
        using scalarField = std::shared_ptr<Sfx::FieldLattice>;
    public:
        EmcilAlgorithm( std::string aName);
        void calculate( Sfx::PhysicalSystem* aPhySys) override;
 
  protected: // can be overwritten by subclass, don't see why? 
    void traitementTermeS( const Sfx::PhysicalSystem* aPsys, const Sfx::EMcNeilBndCnd& aBC, RetRhsType* aRetRHS);    
    void traitementTermeP( const Sfx::PhysicalSystem* aPsys, const Sfx::EMcNeilBndCnd& aBC, RetRhsType* aRetRHS);
    void calculFF( const Sfx::PhysicalSystem* aPsys, const Sfx::EMcNeilBndCnd& aBC, RetRhsType* aRetRHS);
 
    private:
    std::string m_Name;                    
    Sfx::LDeltaOperator::SWE_RHS m_sweRhs; 
    };
} // End of namespace