GerstnerWaves.h

 
 
#ifndef GERSTNER_WAVES
#define GERSTNER_WAVES
 
#include "LinearWaveTheory.h"
class ISimObjectCreator;
 
class GerstnerWaves : public LinearWaveTheory
{
public:
    GerstnerWaves( bool bIrregularWaves, vector<double> & mdWaveAmplitude, vector<double> & mdPhaseAngle, vector<double> & adWaveFrequency, vector<double> & adWaveDirection, ISimObjectCreator* creator);
    ~GerstnerWaves(void);
 
    double      ComputeSurfaceElevation(const double dT, const double adPos[2]);
    void        ComputeSurfaceUnitNormal(const double dT, const double adPos[2], double *adUnitNormalOut);
    void        ComputeParticleVelocity(const double dT, const double adPos[3], double *adVelOut);
    void        ComputeParticleAcceleration(const double dT, const double adPos[3], double *adAccOut);
    double      ComputeDynamicPressure(const double dT, const double adPos[3]);
    double      ComputeTotalPressure(const double dT, const double adPos[3]);
 
    virtual void ComputePointQuery(double T, const double Position[3], double& surface_elevation, double surface_unit_normal[3], double particle_velocity[3], double particle_acceleration[3], double& pressure);
 
 
protected:
 
    void FindLagrangePoint3D(const double dT, const double adEulerPointIn[3], double *adLagrangePointOut);
    void FindLagrangePoint2D(const double dT, const double adEulerPointIn[2], double *adLagrangePointOut);
 
};
 
#endif