LinearWaveTheory.h

 
 
#ifndef LINEAR_WAVE_THEORY
#define LINEAR_WAVE_THEORY
 
#include <vector>
using std::vector;
 
 
 
class LinearWaveTheory
{
public:
    
    enum EWaveTheory_t
    {
        EWaveTheory_Airy=0,
        EWaveTheory_Gerstner,
    };
 
    LinearWaveTheory(bool bIrregularWaves, vector<double> & mdWaveAmplitude, vector<double> & mdPhaseAngle, vector<double> & adWaveFrequency, vector<double> & adWaveDirection);
    virtual ~LinearWaveTheory(){}
 
    virtual double  ComputeSurfaceElevation(const double dT, const double adPos[2]) = 0;
    virtual void    ComputeSurfaceUnitNormal(const double dT, const double adPos[2], double *adUnitNormalOut) = 0;
    virtual void    ComputeParticleVelocity(const double dT, const double adPos[3], double *adVelOut) = 0;
    virtual void    ComputeParticleAcceleration(const double dT, const double adPos[3], double *adAccOut) = 0;
    virtual double  ComputeDynamicPressure(const double dT, const double adPos[3]) = 0;
    virtual double  ComputeTotalPressure(const double dT, const double adPos[3]) = 0;
 
    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) = 0;
    
    std::vector<double>     m_mdZetaA;                  
    std::vector<double>     m_mdEpsilon;                
    std::vector<double>     m_adOmega;                  
    std::vector<double>     m_adWaveNumber;             
    std::vector<double>     m_adTheta;                  
 
    void   SetWaveHarmonicCuttoffPercentage(double percent);    
    double GetWaveHarmonicCuttoffPercentage();          
 
protected:
    double                  m_ln_percent;               
    std::vector<double>     m_adCosTheta;               
    std::vector<double>     m_adSinTheta;               
    double                  m_dGravity;                 
    double                  m_dDensity;                 
    int                     m_NumWaves;                 
//  bool                    m_bIrregularWaves;          /// Irregular wave identifier
 
};
 
#endif