VirtualCurrentInDir.h

#pragma once
 
#include <SimObject.h>
#include <CEnvironment.h>
#include <memory>
 
class VirtualCurrentInDir : public SimObject, ICurrentField
{
public:
  VirtualCurrentInDir(const string& simObjectName, ISimObjectCreator* const creator);
 
  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;
  CEnvironment  *Environment;
 
  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
};