diff options
| author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-05-07 23:07:23 -0300 |
|---|---|---|
| committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-05-07 23:07:23 -0300 |
| commit | e3f311a91b9684924262108a4e0c8e934f5d1d70 (patch) | |
| tree | 01d989abdc354ab28e9a4a551052039972d2cc42 /log-fourier_integrator.cpp | |
| parent | 1062dfa43c4206760fa02bf0c4a971fd782cb75c (diff) | |
| download | code-e3f311a91b9684924262108a4e0c8e934f5d1d70.tar.gz code-e3f311a91b9684924262108a4e0c8e934f5d1d70.tar.bz2 code-e3f311a91b9684924262108a4e0c8e934f5d1d70.zip | |
Revert "Average first two values of C to fix μ"
This reverts commit b9e5a1f2448aeaca8aa807f3baccbe66f5dc49d2.
Diffstat (limited to 'log-fourier_integrator.cpp')
| -rw-r--r-- | log-fourier_integrator.cpp | 3 |
1 files changed, 1 insertions, 2 deletions
diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp index 017def7..55f358a 100644 --- a/log-fourier_integrator.cpp +++ b/log-fourier_integrator.cpp @@ -140,8 +140,7 @@ int main(int argc, char* argv[]) { std::vector<Real> Rₜ₊₁ = fft.inverse(Ȓₜ₊₁); std::vector<Real> Cₜ₊₁ = fft.inverse(Ĉₜ₊₁); - Real C₀ = (Cₜ₊₁[0] + Cₜ₊₁[1]) / 2; - μₜ *= pow(tanh(C₀-1)+1, x); + μₜ *= pow(tanh(Cₜ₊₁[0]-1)+1, x); ΔCₜ = 0; for (unsigned i = 0; i < N; i++) { |