Removed some more duplicate code.

1 files changed, 1 insertions, 3 deletions

diff --git a/get_energy.cpp b/get_energy.cpp
index a388e0c..b1527dc 100644
--- a/get_energy.cpp
+++ b/get_energy.cpp
@@ -75,9 +75,7 @@ int main(int argc, char* argv[]) {
       std::vector<Complex> Ct = fft.fourier(C);
       std::vector<Complex> Rt = fft.fourier(R);
 
-      auto [RddfCt, dfCt] = RddfCtdfCt(fft, C, R, p, s, λ);
-
-      Real z = ((std::conj(Rt[0]) + pow(y, 2) * (RddfCt[0] * Ct[0] + dfCt[0] * std::conj(Rt[0]))) / Ct[0]).real();
+      Real z = estimateZ(fft, C, Ct, R, Rt, p, s, λ, y);
 
       std::cout << y << " " << e << " " << Ct[0].real() << " " << z << std::endl;
     }