Changeset 3965 in Sophya for trunk

Timestamp:

Mar 4, 2011, 5:53:32 PM (15 years ago)

Author:

ansari

Message:

Bug reel/imag ds accumulation en frequence visfits2dt.cc, Cedric+Reza, 04/03/2011

File:

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

trunk/AddOn/TAcq/visfits2dt.cc

@@ -200 +200 @@
       ncols = vismtx.NCols();
       cout << " ProcVisMtxFiles/Info: Output Time-Frequency matrices NRows(time) " 
-           << nrows/mtf_bintime+1 << " NCols(freq) =" << nfreq << endl;
+           << nrows/mtf_bintime+1 << " NCols=2*NFreq) =" << nfreq << endl;
       for(size_t j=0; j<rowlist.size(); j++)   
-        vmtf.push_back(TMatrix< r_4 >(nrows/mtf_bintime+1, nfreq));
+        vmtf.push_back(TMatrix< r_4 >(nrows/mtf_bintime+1, nfreq*2));
     }
 
@@ -208 +208 @@
       for(size_t j=0; j<rowlist.size(); j++) {
         sa_size_t jfreb=0;
-        for(sa_size_t jf=jf1; jf<jf2; jf+=djf) {
+        for(sa_size_t jf=jf1; jf<jf2; jf++) {
           sa_size_t jfreb=(jf-jf1)/djf;
           if (jfreb>=nfreq)  break;
-          vmtf[j](ktime/mtf_bintime,jfreb)+=vismtx(rowlist[j],jf);
+          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 
         }
       }