VirtualCurrentInDir.h

#pragma once
 
#include <SimObject.h>
#include <EnvironmentProvider.h>
#include <memory>
 
class VirtualCurrentInDir;
 
class VirtualCurrentFieldInDir : public marenv::current::CurrentField {
public:
   VirtualCurrentFieldInDir() = default;
   virtual ~VirtualCurrentFieldInDir() = default;
   void inline setVirtualCurrentInDir(VirtualCurrentInDir* virtCurr){ m_virtCurr = virtCurr; }
   [[deprecated("This may not use up-to-date state info")]]
   marenv::Status GetCurrent(const double T, const double position[3], double currentVelocity[3]);
private:
   VirtualCurrentInDir* m_virtCurr;
};
 
class VirtualCurrentInDir : public SimObject
{
friend class VirtualCurrentFieldInDir;
public:
  VirtualCurrentInDir(const string& simObjectName, ISimObjectCreator* const creator);
  virtual ~VirtualCurrentInDir();
  void FinalSetup(const double T, const double *const X, ISimObjectCreator* const creator);
  void InitialConditionSetup(const double T, const double* const currentIC, double* const updatedIC, ISimObjectCreator* const creator);
  void GetCurrent(const double T, const double* const X, const double position[3], double* currentVelocity);
  void OdeFcn(const double T, const double* const X, double* const XDot, const bool bIsMajorTimeStep);
 
  void PreOdeFcn(const double T, const double *const X, IStateUpdater* updater);    
 
#ifdef FH_VISUALIZATION
  void RenderInit(Ogre::Root* const ogreRoot, ISimObjectCreator* const creator);
  void RenderUpdate( const double T, const double* const X);
#endif
protected:
  void SetOutputPortValues(const double T, const double* const X);
  const double * PlanarCurrent(const double T, const double * const X, int layer);
private:
  
  ICommonComputation* m_SetOutputPortValues;
  double m_Pp, m_Pd, m_Pi, m_IntegralSaturation, m_DerivativeTInv;
  int m_CurrentIndex, m_IntegralIndex, m_DerivativeIndex, m_layers;
 
  std::unique_ptr<double[]> m_PlanarCurrents;
  environment::EnvironmentProvider  *Environment;
  std::shared_ptr<VirtualCurrentFieldInDir> m_CurrentField = std::make_shared<VirtualCurrentFieldInDir>();
 
  std::vector<int> m_CurrentIndices, m_IntegralIndices, m_DerivativeIndices;
  std::unique_ptr<double[]> m_Depths;
 
  double m_MaximumCurrent;      
  std::unique_ptr<double[]> m_Directions;       
  
  double m_MaximumCurrentRel;       
  
  double m_MinimumCurrentTot;       
 
  bool  m_AdaMax;
  std::unique_ptr<double[]> m_AdaMaxRatio;
  std::unique_ptr<double[]> m_AdaErr;
  std::unique_ptr<double[]> m_AdaErrTol;
 
protected:
  std::unique_ptr< ISignalPort*[] > m_Input_e;
  std::unique_ptr< ISignalPort*[] > m_Input_d;
  std::unique_ptr< ISignalPort*[] > m_Input_v;
  std::unique_ptr< ISignalPort*[] > m_Input_vt;
 
#ifdef FH_VISUALIZATION
  Ogre::SceneNode*    m_RenderNode;  
#endif
};