Changeset 3893 in Sophya for trunk/AddOn/TAcq/brviscalc.cc

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

File:

• trunk/AddOn/TAcq/brviscalc.cc (modified) (3 diffs)

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));
       }