Changeset 3703 in Sophya for trunk/AddOn/TAcq/svv2mtx.cc

Timestamp:

Dec 6, 2009, 12:19:14 AM (16 years ago)

Author:

ansari

Message:

Correction bug ds svv2mtx.cc, Reza 5/12/2009

File:

• trunk/AddOn/TAcq/svv2mtx.cc (modified) (4 diffs)

trunk/AddOn/TAcq/svv2mtx.cc

@@ -153 +153 @@
 {
   Timer tm("ProcSVFiles");
-
-  vector< TMatrix< complex<r_4> > > vmtf(rowlist.size());
-
+  
   DataTable visdt;
   visdt.AddDoubleColumn("mfc");
@@ -164 +162 @@
   visdt.AddFloatColumn("vim");
 
+  vector< TMatrix< complex<r_4> > > vmtf;
+
   if (jf1<1)  jf1=1;
   if ((jf2<1)||(jf2<jf1))  jf2=jf1; 
@@ -169 +169 @@
   int djf=(jf2-jf1)/nfreq;
   if (djf<1) djf=0;
-  
+
   char fname[1024];
 
+  cout << " ---- ProcSVFiles()-Begin - Reading chanum vector" << endl;
   TVector< uint_4 > chanum;
+  {
   sprintf(fname, "%s/chanum.ppf",inoutpath.c_str());
-  {
   PInPersist pic(fname);
   pic >> chanum;
   }
 
-  sa_size_t nrows = (imax-imin+1)/istep;\
+  cout << " ---- ProcSVFiles()-Read chanum done " << endl;
+  sa_size_t nrows = (imax-imin+1)/istep;
   sa_size_t kr=0;
 
@@ -193 +195 @@
       cout << " ProcSVFilesVJun09/Info: Output Time-Frequency matrices NRows=NFiles" 
            << nrows << " NCols=NFreq=" << ncols << endl;
-      for(size_t j=0; j<rowlist.size(); j++)   vmtf[j].SetSize(nrows, ncols);
+      for(size_t j=0; j<rowlist.size(); j++)   
+        vmtf.push_back(TMatrix< complex<r_4> >(nrows, ncols) );
     }
     for(size_t j=0; j<rowlist.size(); j++)