Testvs19_Nujic2StepsIntegrator.h

Go to the documentation of this file.

 
#pragma once
 
// BaseNumType includes
#include "Sfx/Sfx_NodalValues.h"
#include "Sfx/Sfx_FieldLattice.h"
//  SfxBase include
#include "numeric/Sfx_TwoStepIntegrator.h"
// Swe package includes
#include "../../SfxBase/Sfx_ImposeBnd.h"
#include "../../SfxBase/Sfx_Logger.h"
// ...
#include "Testvs19_RetRhsType.h"
 
namespace Testvs19 
{
  class Nujic2StepsIntegrator final : public Sfx::TwoStepIntegrator
  {
    using pairofield = std::pair<Sfx::FieldLattice, Sfx::FieldLattice>;
 
  public:
    Nujic2StepsIntegrator(const Sfx::FieldLattice& aU1, const Sfx::FieldLattice& aU2);
    Nujic2StepsIntegrator( const Sfx::FieldLattice& aU1,
      const Sfx::FieldLattice& aU2, Sfx::LDeltaOperator* aLdeltaOp);
 
    //
    // TwoStepIntegrator interface
    //

    void initialize() override final {};
    void predictor() override final;
    void corrector() override final;
    RetRhsType* previousState() const override { return nullptr; /*RetRhsType{ m_midvlArrdU1,  m_midvlArrdU1 }.asStdVector();*/ }
    RetRhsType* currentState() const override 
    { 
      if (getStepType() == TwoStepIntegrator::eStepType::predictorStep)
      {
        return new RetRhsType{ m_midU1.values().asStdVector(), m_midU2.values().asStdVector() };
      }
      else
      {
        return new RetRhsType{ m_currU1.values().asStdVector(), m_currU2.values().asStdVector() };
      }
    }

    pairofield getFinalSolution() const noexcept
    {
      // C++17 copy elison is mandatory??? use copy ctor to create field return type?? 
      // return by value!
      return std::make_pair(Sfx::FieldLattice{ m_currU1 }, Sfx::FieldLattice{ m_currU2 });
    }
    // C++17 mandatory copy Elision (materialize temporaries to xValue) no copy!!!!
//     RetRhsType::pairofstdvec midState() const { return RetRhsType { m_midvlArrdU1,  m_midvlArrdU2 }.asStdVector(); }
//     RetRhsType::pairofstdvec curState() const { return RetRhsType { m_currvlArrU1,  m_currvlArrU2 }.asStdVector(); }
 
    //RetRhsType rhsmidState() const { return RetRhsType{ m_midU1State, m_midU2State }; }
    //RetRhsType rhscurState() const { return RetRhsType{ m_curU1State,  m_curU2State }; }
 
    // this is simple, we want final state of the two-step algo
    // Note sure about this one. I added this implementation in
    // the first implementation of the NujicIntegrator.
//     virtual pairofield getCurrState() const override { return std::make_pair({ nullptr,nullptr }); }
//     virtual pairofield getPrevState() const override { return std::make_pair(nullptr, nullptr); }
    
    void setInitialState( const Sfx::FieldLattice& aU1, 
      const Sfx::FieldLattice& aU2) override final;
    void setBCAtBothEnds( const Sfx::EMcNeilBndCnd::pair_node& aBcNodes) 
    {
      m_upstreamNode = aBcNodes.first;
      m_downstreamNode = aBcNodes.second;
    }

    void setBCValuesAtBothEnds( const std::pair<Sfx::Object*, Sfx::Object*>& aBcNodes) override final
    {
      m_upstreamNode = *static_cast<Sfx::Nodal_Value*>(aBcNodes.first);
      m_downstreamNode = *static_cast<Sfx::Nodal_Value*>(aBcNodes.second);
    }

    void step();
 
   private:
     Sfx::LDeltaOperator* m_deltaOp;   
     // Design Note
     //  i don't see why declaring a scalarField when we need a numerical array 
     //  just make things more complicated (predictor/corrector) want to add 
     //  values and with numerical array that support math operations its easier
     Sfx::FieldLattice m_midU1;         
     Sfx::FieldLattice m_midU2;         
     Sfx::FieldLattice m_currU1;           
     Sfx::FieldLattice m_currU2;        
     Sfx::Nodal_Value m_upstreamNode;   
     Sfx::Nodal_Value m_downstreamNode; 
  };
} // End of namespace