Ambient.h

#include <complex>
#include <random>
#include <point3.h>
#include <iparamb2.h>
#include "Grid.h"

typedef std::complex<float> complex;

class Ambient : public Grid {
    std::tr1::mt19937 engine;
    std::tr1::normal_distribution<float> dist;

    const float         gravity;                    
    const float         omega_0;                    
    const float         T;                          
    const unsigned long seed;                       
    float               Lx;                         
    float               Ly;                         
    float               l;                          
    float               V;                          
    Point3              w_hat;                      
    float               t;                          
    // Values precached on initialization.
    float               P_h__L;                     
    float               P_h__l_2;                   
    complex*            h_tildes_in;
    complex*            h_tildes_out;

public:
    static const float  GRAVITY_METRIC;             
    static const float  GRAVITY_US;                 
    Ambient(float width, float length, int widthSegs, int lengthSegs, unsigned long rngSeed, float phaseDuration, float accelerationGravity);

    ~Ambient();

    void                Simulate(float time, float speed, float direction, float scale, float waveSizeLimit, float desiredMaxHeight);

private:
    float               omega(Point3 k);

    float               P_h(Point3 k);

    complex             h_tilde_0(Point3 k);

    complex             h_tilde(Point3 k);

    float               max_height();

    void                heights(float heightScale, float time);
};