diff options
Diffstat (limited to 'log-fourier_integrator.cpp')
| -rw-r--r-- | log-fourier_integrator.cpp | 8 |
1 files changed, 6 insertions, 2 deletions
diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp index eb72b49..5dc64ec 100644 --- a/log-fourier_integrator.cpp +++ b/log-fourier_integrator.cpp @@ -4,15 +4,18 @@ #include <iostream> int main(int argc, char* argv[]) { + /* Model parameters */ unsigned p = 2; unsigned s = 2; Real λ = 0.5; Real τ₀ = 0; + /* Log-Fourier parameters */ unsigned log2n = 8; Real Δτ = 0.1; Real k = 0.1; + /* Iteration parameters */ Real ε = 1e-13; Real γ = 1; Real βₘₐₓ = 0.7; @@ -69,7 +72,7 @@ int main(int argc, char* argv[]) { std::vector<Complex> Ĉₜ(N); std::vector<Complex> Ȓₜ(N); - // start from the exact solution for β = 0 + /* Start from the exact solution for β = 0 */ for (unsigned n = 0; n < N; n++) { if (τ₀ > 0) { Cₜ[n] = Γ₀ * (exp(-μ * fft.t(n)) - μ * τ₀ * exp(-fft.t(n) / τ₀)) / (μ - pow(μ, 3) * pow(τ₀, 2)); @@ -135,7 +138,8 @@ int main(int argc, char* argv[]) { it = 0; } -// std::cerr << ΔC << std::endl; + std::cerr << β << " " << μ << " " << ΔC << " " << γ; + std::cerr << "\r"; } /* Integrate the energy using Simpson's rule */ |