summaryrefslogtreecommitdiff
path: root/log-fourier_integrator.cpp
diff options
context:
space:
mode:
authorJaron Kent-Dobias <jaron@kent-dobias.com>2025-05-18 18:32:16 -0300
committerJaron Kent-Dobias <jaron@kent-dobias.com>2025-05-18 18:32:16 -0300
commite4a2618e4c25e528850565d5e75e5d016df64c4e (patch)
treed9b5b914c6f07185c02c6c8cef46e7522013ba4b /log-fourier_integrator.cpp
parent2653586bab4f9519e18be93caea2dea7f891a8e6 (diff)
downloadcode-e4a2618e4c25e528850565d5e75e5d016df64c4e.tar.gz
code-e4a2618e4c25e528850565d5e75e5d016df64c4e.tar.bz2
code-e4a2618e4c25e528850565d5e75e5d016df64c4e.zip
Small notation change
Diffstat (limited to 'log-fourier_integrator.cpp')
-rw-r--r--log-fourier_integrator.cpp6
1 files changed, 3 insertions, 3 deletions
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<Real> Cₜ₊₁ = fft.inverse(Ĉₜ₊₁);