LocalWaveVessel.h

#ifndef LocalWaveVessel_H
#define LocalWaveVessel_H
 
#include <array>
#include <fhsim/simobject/SimObject.h>
#include <fhsim_environment/EnvironmentProvider.h>
#include <string>
#include <vector>
 
class LocalWaveVessel : public SimObject
{
public:
   using Dof = std::array<double, 3>;
 
   LocalWaveVessel(std::string simObjectName, ISimObjectCreator* creator);
 
   void FinalSetup(double T, const double* X, ISimObjectCreator* creator) override;
 
   void OdeFcn(double T, const double* X, double* XDot) const override;
 
   bool HasJacobians() const override { return true; }
   void OdeJacobian(double T, const double* X, double* J, int nStates) override;
   int  GetJacobianSparsity(int nStates, int* rowPtr, int* colIdx) override { return -1; }
 
#ifdef FH_VISUALIZATION
   void RenderInit(Ogre::Root* ogreRoot, ISimObjectCreator* creator) override;
   void RenderUpdate(double T, const double* X) override;
#endif
 
   Dof RawWaveLoads(double T, const double* X) const;
 
   const double* MotionPort(double T, const double* X);
   const double* MotionRatePort(double T, const double* X);
   const double* ExcitationPort(double T, const double* X);
   const double* PositionPort(double T, const double* X);
   const double* RotationPort(double T, const double* X);
 
private:
   struct SamplePoint
   {
      double x;
      double y;
   };
 
   void   BuildFootprint(double length, double breadth, double lcf, int sampleLevel);
   double RotationalStiffness(ISimObjectCreator* creator, const std::string& gmName, const std::string& periodName,
         double mass, double inertia) const;
 
   Dof Accelerations(const double* motion, const double* rate, const double* load) const;
 
   ISignalPort*                      m_inPosition  = nullptr;
   ISignalPort*                      m_inRotation  = nullptr;
   environment::EnvironmentProvider* m_environment = nullptr;
 
   int m_indexMotion     = 0;
   int m_indexRate       = 0;
   int m_indexExcitation = 0;
 
   double m_rho                  = 0.0;
   double m_waterplaneArea       = 0.0;
   double m_filterTime           = 0.0;
   double m_rollQuadraticDamping = 0.0;
 
   Dof    m_inertia {};   
   Dof    m_stiffness {}; 
   Dof    m_damping {};   
   double m_a35 = 0.0;
   double m_b35 = 0.0;
   double m_c35 = 0.0;
 
   std::vector<SamplePoint> m_samples;
   double                   m_sampleWeight = 0.0; 
 
   std::array<double, 3> m_positionOut {}; 
   std::array<double, 3> m_rotationOut {}; 
 
   double m_length  = 0.0;
   double m_breadth = 0.0;
   double m_lcf     = 0.0;
   double m_draught = 0.0; 
 
#ifdef FH_VISUALIZATION
   Ogre::SceneNode* m_renderNode = nullptr;
#endif
};
 
#endif