Changeset 3893 in Sophya

Timestamp:

Sep 28, 2010, 5:49:39 PM (15 years ago)

Author:

ansari

Message:

Correction bug calcul du numero de paires (V1,V2) lors de la boucle sur le calcul de visibilites, Rz+cmv, 28/09/2010

Location:

trunk/AddOn/TAcq

Files:

• brviscalc.cc (modified) (3 diffs)
• brviscalc.h (modified) (1 diff)
• vismfib.cc (modified) (1 diff)

trunk/AddOn/TAcq/brviscalc.cc

@@ -326 +326 @@
   }
 
-  sa_size_t k=0;
+// kpair=numero sequentiel de la paire: 0->(0,0), 1->(0,1), 2->(0,2), 3->(0,3), 4->(1,1), 5->(1,2) ...
+  sa_size_t kpair=0;  
   for(size_t i=0; i<vpdata_.size(); i++) {
     for(size_t j=i; j<vpdata_.size(); j++) {
-      size_t kpair=i*vpdata_.size()+j;
       if (kpair<pairst_)  continue;
       if (kpair>=(pairst_+nbpairs_))  break;
-      TVector< complex<r_4> > vis = vismtx_.Row(k);   k++;
+      TVector< complex<r_4> > vis = vismtx_.Row(kpair);   
       for(sa_size_t f=1; f<vis.Size(); f++) {
         vis(f) += complex<r_4>((r_4)vpdata_[i][f].realB(), (r_4)vpdata_[i][f].imagB()) * 
           complex<r_4>((r_4)vpdata_[j][f].realB(), -(r_4)vpdata_[j][f].imagB());
       }
+      kpair++;
       nb_flop_ += (7.*(r_8)vis.Size());
     }
@@ -399 +400 @@
       pvpdata[2*fib+1] = ppaq.Data2C();
     }
-    sa_size_t k=0;
+// kpair=numero sequentiel de la paire: 0->(0,0), 1->(0,1), 2->(0,2), 3->(0,3), 4->(1,1), 5->(1,2) ...
+    sa_size_t kpair=0;  
     for(size_t i=0; i<vpdata_.size(); i++) {
       for(size_t j=i; j<vpdata_.size(); j++) {
@@ -405 +407 @@
         if (kpair<pairst_)  continue;
         if (kpair>=(pairst_+nbpairs_))  break;
-        TVector< complex<r_4> > vis = vismtx_.Row(k);   k++;
+        TVector< complex<r_4> > vis = vismtx_.Row(kpair);  
         for(sa_size_t f=fdeb; f<ffin; f++) {
           vis(f) += complex<r_4>((r_4)pvpdata[i][f].realB(), (r_4)pvpdata[i][f].imagB()) * 
             complex<r_4>((r_4)pvpdata[j][f].realB(), -(r_4)pvpdata[j][f].imagB());
         }
+        kpair++;
         nb_flop_ += (7.*(r_8)(ffin-fdeb));
       }

trunk/AddOn/TAcq/brviscalc.h

@@ -40 +40 @@
   void DefineRank(uint_4 nbc, uint_4 cid);
 
-  // Pour limiter le calcul des visibilites a certaines fibres 
+  // Pour limiter le calcul des visibilites a certaines voies, determine par les numeros de sequences des paires 
+  // Exemple pour 2 fibres et 4 voies 
+  // kpair=numero sequentiel de la paire: 0->(0,0), 1->(0,1), 2->(0,2), 3->(0,3), 4->(1,1), 5->(1,2) ...
   uint_4 SelectPairs(uint_4 pair1=0, uint_4 nbpairs=99999999);
   // pour definir le binning en frequence pour la VisibDataTable 

trunk/AddOn/TAcq/vismfib.cc

@@ -74 +74 @@
     }
     BRMultiFitsReader reader(mmgr, par.dirlist_, true, par.imin_, par.imax_, par.istep_);
-    reader.SetPrintLevel(par.prtlevel_);
+    reader.SetPrintLevel(par.prtlevel_,par.prtmodulo_);
 
     /*    BRVisibilityCalculator proc(mmgr, par.outpath_, par.nmean_, par.nthreads_);