WaveEnergySpecter.h

 
 
#ifndef WAVE_ENERGY_SPECTER
#define WAVE_ENERGY_SPECTER
 
#include <random>
#include <vector>
 
 
class WaveEnergySpecter
{
public:
    enum EWaveSpecter_t
    {
        EWaveSpecter_Jonswap=0,
        EWaveSpecter_Issc,
    };
    enum EDirectionSpecter_t
    {
        EDirectionSpecter_None=0,
        EDirectionSpecter_CosineSquared,
    };
 
    struct SSpecterParameters {
            double  dHs;        
            double  dT1;        
            double  dTheta; 
            double  dCosinePower; 
            bool    dIsLongCrested; 
            double  dOmegaNDimLowLim; 
            double  dOmegaNDimUppLim; 
    };
 
    WaveEnergySpecter(EWaveSpecter_t EWaveSpecter, EDirectionSpecter_t EDirectionSpecter, SSpecterParameters* pSSpecterParameters, int iNumWaves, unsigned int RandomSeed);
    ~WaveEnergySpecter(void);
 
 
    //double Cosinesquared      (double theta, double theta1, double theta2);
    //double RandomPhaseAngle       ();
    //void ComputeWaveAmplitudes    ();
 
 
    const std::vector<double>*  GetWaveFrequencyArray();
    const std::vector<double>*  GetWaveDirectionArray();
    const std::vector<double>*  GetSpectralValueArray();
    const std::vector<double>*  GetRandomPhaseArray();
    const std::vector<double>*  GetWaveNumberArray();
 
 
 
protected:
    typedef double (WaveEnergySpecter::*NormalizedProbabilityDistribution)(const double); 
    void fillRandomVector(NormalizedProbabilityDistribution f, double limit_low, double limit_high, std::vector<double>& result);
    void makeWaveFrequencyVectorWithRandomIntervals();
    double randomNumber(double limit_low, double limit_high);
 
    double cosPower(double theta);
 
 
    double SpectralEnergy       (double omega);
    double Jonswap              (double omega);
    double NormalizedConstantSteepnessJonswap   (double omega);
    double Issc                 (double omega);
    double SpectralDirectionDistribution    (double theta);
    double NormalizedConstantSteepnessIssc(double omega);
 
    // Member variables
    double                  m_dHs;      
    double                  m_dT1;      
    double                  m_cosPower;
    bool                    m_isLongCrested; 
    double                  m_mainWaveDirection;
    //double                    m_dTheta1;  /// Lower sector angle of direction spectrum
    //double                    m_dTheta2;  /// Upper sector angle of direction spectrum
    double                  m_dRho;     
    double                  m_dGravity; 
    double                  m_omgLowLim_nDim;   
    double                  m_omgUppLim_nDim;   
    
 
private:
    EWaveSpecter_t          m_EWaveSpecter;
    EDirectionSpecter_t     m_EDirectionSpecter;
    std::vector<double>     m_ZetaA; 
    std::vector<double>     m_Epsilon; 
    std::vector<double>     m_Omega; 
    std::vector<double>     m_deltaOmega; 
    std::vector<double>     m_WaveNumber; 
    std::vector<double>     m_Theta; 
    int                     m_iNumWaves;
 
    std::mt19937 m_rng;
    std::uniform_real_distribution<double> m_dist;
};
 
#endif