From e4a2618e4c25e528850565d5e75e5d016df64c4e Mon Sep 17 00:00:00 2001 From: Jaron Kent-Dobias Date: Sun, 18 May 2025 18:32:16 -0300 Subject: Small notation change --- log-fourier_integrator.cpp | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) (limited to 'log-fourier_integrator.cpp') diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp index f2a88c2..f553569 100644 --- a/log-fourier_integrator.cpp +++ b/log-fourier_integrator.cpp @@ -122,7 +122,7 @@ int main(int argc, char* argv[]) { Rₜ₋₁[n] = std::exp(-μ₀ * fft.t(n)); Ĉₜ₋₁[n] = 2 * Γ₀ / (pow(μ₀, 2) + pow(fft.ν(n), 2)) / (1 + pow(τ₀ * fft.ν(n), 2)); - Ȓₜ₋₁[n] = (Real)1.0 / (μ₀ + II * fft.ν(n)); + Ȓₜ₋₁[n] = (Real)1 / (μ₀ + II * fft.ν(n)); } } else { if (!logFourierLoad(Cₜ₋₁, Rₜ₋₁, Ĉₜ₋₁, Ȓₜ₋₁, p, s, λ, τ₀, β₀, log2n, Δτ, logShift)) { @@ -172,12 +172,12 @@ int main(int argc, char* argv[]) { ΔCₜ = 0; for (unsigned n = 0; n < N; n++) { ΔCₜ += std::norm(((2 * Γ[n] + D[n]) * std::conj(Ȓₜ[n]) + Σ[n] * Ĉₜ[n]) - Ĉₜ[n] * (μₜ + II * fft.ν(n))); - ΔCₜ += std::norm(((Real)1.0 + Σ[n] * Ȓₜ[n]) - Ȓₜ[n] * (μₜ + II * fft.ν(n))); + ΔCₜ += std::norm(((Real)1 + Σ[n] * Ȓₜ[n]) - Ȓₜ[n] * (μₜ + II * fft.ν(n))); } ΔCₜ = sqrt(ΔCₜ) / (2*N); for (unsigned n = 0; n < N; n++) { - Ȓₜ₊₁[n] = ((Real)1.0 + Σ[n] * Ȓₜ[n]) / (μₜ + II * fft.ν(n)); + Ȓₜ₊₁[n] = ((Real)1 + Σ[n] * Ȓₜ[n]) / (μₜ + II * fft.ν(n)); } std::vector Cₜ₊₁ = fft.inverse(Ĉₜ₊₁); -- cgit v1.2.3-70-g09d2