dbpp_PhysicalBoundaryCnd.h

Go to the documentation of this file.

#pragma once
 
// C++ include
#include <tuple>
//..
#include "dbpp_PhysicalSystem.h"
#include "dbpp_EMcNeilPhysBC.h"
// ...
//#include "../May_VS2015/Test/NodalVariables.h"
 
// Note
//   this is a first implementaion (still under development)
namespace dbpp 
{
  using float64 = double;

  class PhysicalBoundaryCnd
  {
  public:
    //  template<typename T>
   // using PhyBCNodeConstraint = EMcNeilPhysBC::PhyBCNodeConstraint;
  
  public:
    PhysicalBoundaryCnd() 
    : m_phsys{ nullptr },
      m_leftEnd{},
      m_rightEnd{}
    {}

    PhysicalBoundaryCnd(PhysicalSystem* aPhySys);
    PhysicalBoundaryCnd(PhyBCNdlConstraint aLeftBC,
      PhyBCNdlConstraint aRightBC);
    PhyBCNdlConstraint getLeftEnd() const
    { 
      // characteristic s
      //PhysicalBCNodeConstraint w_ndlVal;  created a tuple by default
      return m_leftEnd;
    }

    PhyBCNdlConstraint getRightEnd() const
    { 
      // characteristic s
      //PhysicalBCNodeConstraint w_ndlVal;  created a tuple by default
      return m_rightEnd;
    }

    void update();

    PhysicalSystem* getPhysicalSystem() const noexcept { return m_phsys; }
    virtual void setLeftBC( const std::tuple<unsigned, float64, float64, float64>& aLeftPhyBnd)
    {
      // add code here
      // compute from E. McNeil algorithm (characteristic)
      m_leftEnd = std::make_from_tuple<PhyBCNdlConstraint>(aLeftPhyBnd);
    }

    virtual void setRightBC(std::tuple<unsigned, float64, float64, float64> aRightPhyBnd)
    {
      // add code here
      // compute from E. McNeil algorithm (characteristic)
      m_rightEnd = std::make_from_tuple<PhyBCNdlConstraint>(aRightPhyBnd);
    }
  private:
    PhysicalSystem* m_phsys;        
    // Characteristic solver method   
    PhyBCNdlConstraint m_leftEnd;  
    PhyBCNdlConstraint m_rightEnd; 
  };
} // End of namespace