Changeset 3966 in Sophya for trunk

Timestamp:

Mar 4, 2011, 6:16:35 PM (15 years ago)

Author:

ansari

Message:

Ajout calcul module^{2 moyenne ds visfits2dt.cc, Cedric+Reza, 04/03/2011}

File:

• trunk/AddOn/TAcq/visfits2dt.cc (modified) (5 diffs)

trunk/AddOn/TAcq/visfits2dt.cc

@@ -161 +161 @@
 
   vector< TMatrix< r_4 > > vmtf;
+  vector< TMatrix< r_4 > > vmtfa;
 
   if (jf1<1)  jf1=1;
@@ -194 +195 @@
     FitsInOutFile fin(fname, FitsInOutFile::Fits_RO);
     TMatrix< r_4 > vismtx;
+
     fin >> vismtx;
     //    cout << vismtx.Info();
@@ -201 +203 @@
       cout << " ProcVisMtxFiles/Info: Output Time-Frequency matrices NRows(time) " 
            << nrows/mtf_bintime+1 << " NCols=2*NFreq) =" << nfreq << endl;
-      for(size_t j=0; j<rowlist.size(); j++)   
+      for(size_t j=0; j<rowlist.size(); j++) {
         vmtf.push_back(TMatrix< r_4 >(nrows/mtf_bintime+1, nfreq*2));
+        vmtfa.push_back(TMatrix< r_4 >(nrows/mtf_bintime+1, nfreq));
+      }
     }
 
@@ -213 +217 @@
           vmtf[j](ktime/mtf_bintime,jfreb*2)+=vismtx(rowlist[j],jf*2); // partie reelle 
           vmtf[j](ktime/mtf_bintime,jfreb*2+1)+=vismtx(rowlist[j],jf*2+1); // partie imaginaire 
+          vmtfa[j](ktime/mtf_bintime,jfreb)+=vismtx(rowlist[j],jf*2)*vismtx(rowlist[j],jf*2); // module carre - z.real^2 
+          vmtfa[j](ktime/mtf_bintime,jfreb)+=vismtx(rowlist[j],jf*2+1)*vismtx(rowlist[j],jf*2+1); // module carre z.imag^2 
         }
       }
@@ -233 +239 @@
       fo << vmtf[j];
     }
+    sprintf(fname, "!%s/acortfreqmtx.fits",inoutpath.c_str());
+    cout << "ProcVisMtxFiles: Opening file " << fname << " for writing Visi(Freq,Time) matrices" << endl;  
+    FitsInOutFile foa(fname, FitsInOutFile::Fits_Create);
+    for(size_t j=0; j<rowlist.size(); j++)  {
+      cout << " Writing Time-Freq modulesquare matrix[" << j << "] for Mod2IJ= " << chanum(rowlist[j]) << endl;
+      foa << vmtfa[j];
+    }
+
   }
   if (fgfilldt)  {