Commit 0e6457f9 authored by Reza Ansari's avatar Reza Ansari
Browse files

correction bug sauvegarde ds cartes Temps-Frequence

parent eab06fcd
......@@ -327,27 +327,27 @@ int main(int narg, const char* arg[])
POutPersist po(outfile);
po<<dt;
if (FgTFMAC) { // --- renormalizing and saving AutoCorr time-frequency maps
cout<<" visi2dtacx/Info: Saving 8 AutoCorr time-frequency maps to PPF file "<<TFMoutfile<<endl;
if (FgTFMAC || FgTFMCX) { // --- Sauvegarde cartes temps-frequence
POutPersist potfm(TFMoutfile);
const char* tfm_names[8]={"TFM_1H", "TFM_2H", "TFM_3H", "TFM_4H", "TFM_1V", "TFM_2V", "TFM_3V", "TFM_4V"};
for(int k=0; k<8; k++) { // loop over the 8 AutoCorr
TArray< r_4 > & tfmap = vtfmac[k];
tfmap *= (r_4)(1./((double)deltaIavg*(double)TFMfbin));
potfm << PPFNameTag(tfm_names[k]) << tfmap;
if (FgTFMAC) { // --- renormalizing and saving AutoCorr time-frequency maps
cout<<" visi2dtacx/Info: Saving 8 AutoCorr time-frequency maps to PPF file "<<TFMoutfile<<endl;
const char* tfm_names[8]={"TFM_1H", "TFM_2H", "TFM_3H", "TFM_4H", "TFM_1V", "TFM_2V", "TFM_3V", "TFM_4V"};
for(int k=0; k<8; k++) { // loop over the 8 AutoCorr
TArray< r_4 > & tfmap = vtfmac[k];
tfmap *= (r_4)(1./((double)deltaIavg*(double)TFMfbin));
potfm << PPFNameTag(tfm_names[k]) << tfmap;
}
}
}
if (FgTFMCX) { // --- renormalizing and saving H-H Cross-Corr time-frequency maps
cout<<" visi2dtacx/Info: Saving 6 H-H cross-corr time-frequency maps to PPF file "<<TFMoutfile<<endl;
POutPersist potfm(TFMoutfile);
const char* tfm_names[6]={"TFM_1H2H", "TFM_1H3H", "TFM_1H4H", "TFM_2H3H", "TFM_2H4H", "TFM_3H4H"};
for(int k=0; k<6; k++) { // loop over the 6 Xcor
TArray< complex<r_4> > & tfmap = vtfm[k];
tfmap *= complex<r_4>((r_4)(1./((double)deltaIavg*(double)TFMfbin)), 0.);
potfm << PPFNameTag(tfm_names[k]) << tfmap;
if (FgTFMCX) { // --- renormalizing and saving H-H Cross-Corr time-frequency maps
cout<<" visi2dtacx/Info: Saving 6 H-H cross-corr time-frequency maps to PPF file "<<TFMoutfile<<endl;
const char* tfm_names[6]={"TFM_1H2H", "TFM_1H3H", "TFM_1H4H", "TFM_2H3H", "TFM_2H4H", "TFM_3H4H"};
for(int k=0; k<6; k++) { // loop over the 6 Xcor
TArray< complex<r_4> > & tfmap = vtfm[k];
tfmap *= complex<r_4>((r_4)(1./((double)deltaIavg*(double)TFMfbin)), 0.);
potfm << PPFNameTag(tfm_names[k]) << tfmap;
}
}
}
} // --- FIN sauvegarde cartes temps-frequence
// resu.Update();
cout << resu; // Update est fait lors du print
}
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment