diff options
Diffstat (limited to 'src/ocean')
| -rw-r--r-- | src/ocean/ocean_alt.cpp | 2 | ||||
| -rw-r--r-- | src/ocean/ocean_alt.h | 4 |
2 files changed, 3 insertions, 3 deletions
diff --git a/src/ocean/ocean_alt.cpp b/src/ocean/ocean_alt.cpp index 0660a14..e0809a3 100644 --- a/src/ocean/ocean_alt.cpp +++ b/src/ocean/ocean_alt.cpp @@ -569,7 +569,7 @@ std::vector<Eigen::Vector2d> ifft_1d // calculate the twiddle factor int index = group * pow(2, stage) + butterfly; - float w = index * 2 * M_PI / pow(2, stage); + float w = -index * 2 * M_PI / pow(2, stage); Eigen::Vector2d twiddle_factor = {cos(w), sin(w)}; if (reading_buffer1) diff --git a/src/ocean/ocean_alt.h b/src/ocean/ocean_alt.h index 5df1eb5..88344d6 100644 --- a/src/ocean/ocean_alt.h +++ b/src/ocean/ocean_alt.h @@ -85,10 +85,10 @@ private: const double vertex_displacement = Lx / 2; const int N = num_rows*num_cols; // total number of grid points - const double lambda =0.0; // how much displacement matters + const double lambda = .5; // how much displacement matters const double spacing = 1.0; // spacing between grid points - const double A = 6.5; // numeric constant for the Phillips spectrum + const double A = 10; // numeric constant for the Phillips spectrum const double V = 500; // wind speed const double gravity = 9.81; const double L = V*V/gravity; |