diff options
| author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-04-19 15:38:14 -0300 |
|---|---|---|
| committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-04-19 15:38:14 -0300 |
| commit | 00ca38555f905506967180ede7494169102c60cd (patch) | |
| tree | d76a7e474a6e440f1c3b42e9cae19902ae2140e8 /log-fourier_integrator.cpp | |
| parent | ebd38828e38760954affea2c5c9073497c2926a7 (diff) | |
| download | code-00ca38555f905506967180ede7494169102c60cd.tar.gz code-00ca38555f905506967180ede7494169102c60cd.tar.bz2 code-00ca38555f905506967180ede7494169102c60cd.zip | |
More feedback
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 */ |