From 10b10574c65be1267b321ab83694113d80c2f881 Mon Sep 17 00:00:00 2001 From: "J.-S. Caux" Date: Fri, 6 Dec 2019 16:45:10 +0100 Subject: [PATCH] Flip order iK - it in output of S(q,t) --- src/EXECS/Heis_Fourier_to_Sqt.cc | 12 +++++++----- 1 file changed, 7 insertions(+), 5 deletions(-) diff --git a/src/EXECS/Heis_Fourier_to_Sqt.cc b/src/EXECS/Heis_Fourier_to_Sqt.cc index bf09ece..5d70543 100644 --- a/src/EXECS/Heis_Fourier_to_Sqt.cc +++ b/src/EXECS/Heis_Fourier_to_Sqt.cc @@ -92,16 +92,18 @@ int main(int argc, char* argv[]) if (FTim_outfile.fail()) ABACUSerror("Could not open FTim_outfile... "); - for (int iK = iKmin; iK <= iKmax; ++iK) { + for (int it = 0; it <= Npts_t; ++it) { FTre_outfile << setprecision(16); FTim_outfile << setprecision(16); - if (iK > iKmin) { + + if (it > 0) { FTre_outfile << endl; FTim_outfile << endl; } - FTre_outfile << real(FT[0][iK]); - FTim_outfile << real(FT[0][iK]); - for (int it = 1; it <= Npts_t; ++it) { + FTre_outfile << real(FT[it][iKmin]); + FTim_outfile << real(FT[it][iKmin]); + + for (int iK = iKmin + 1; iK <= iKmax; ++iK) { FTre_outfile << "\t" << real(FT[it][iK]); FTim_outfile << "\t" << imag(FT[it][iK]); }