NetCageWakeField.h

#pragma once
 
#include <array>
#include <atomic>
#include <functional>
#include <marenv/AtomicSharedPtr.h>
#include <marenv/marenv.h>
#include <marenv/wake/WakeField.h>
#include <memory>
 
namespace environment
{
 
class NetCageWakeField : public marenv::wake::WakeField
{
public:
   using AmbientCurrentFunction = std::function<marenv::Status(double time, const double pos[3], double vel[3])>;
 
   NetCageWakeField(const double pos[3], double fullWakeDepth, double fullWakeDiameter, double wakeAngle, double initialBoundaryLayer, double maxDeficit, double deficitCutoff, AmbientCurrentFunction ambientCurrent);
   ~NetCageWakeField() override = default;
 
   void SetPosition(const double pos[3]);
 
   marenv::Status Sample(double time, const double pos[3], marenv::wake::WakeSample& out) const override;
 
private:
   using Position = std::array<double, 3>;
 
   double MaxDeficitBehindB(double x) const;
 
   double DeficitInCageFrame(double x, double y, double z) const;
 
   double m_depthA;          
   double m_distAB;          
   double m_radiusFullWakeA; 
   double m_wakeTangentX;    
   double m_tanWakeAngle;    
   double m_radiusNoWakeA;   
   double m_radiusNoWakeB;   
   double m_maxDeficit;      
   double m_bounds_x;        
   double m_bounds_z;        
 
   marenv::AtomicSharedPtr<const Position> m_pos;            
   AmbientCurrentFunction                       m_ambientCurrent; 
};
 
} // namespace environment