Changeset 3921 in Sophya

Timestamp:

Dec 3, 2010, 10:34:23 PM (15 years ago)

Author:

ansari

Message:

tentative d'optimisation du calcul des visibilites, reza 03/12/2010

File:

• trunk/AddOn/TAcq/brviscalc.cc (modified) (1 diff)

trunk/AddOn/TAcq/brviscalc.cc

@@ -415 +415 @@
       if (kpair>=(pairst_+nbpairs_+1))  break;
       if (fgpimp_&&(i!=j)&&((i+j)%2==0))  continue;  // calcul des visib avec numero pair-impair + autocorrel
-      TVector< complex<r_4> > vis = vismtx_.Row(k);  k++;   
+      TVector< complex<r_4> > vis = vismtx_.Row(k);  k++; 
+      TwoByteComplex* zv1=vpdata_[i];
+      TwoByteComplex* zv2=vpdata_[j]; 
+      for(sa_size_t f=1; f<vis.Size(); f++) {
+        vis(f) += complex<r_4>((r_4)zv1->realB(), (r_4)zv1->imagB()) * 
+          complex<r_4>((r_4)zv2->realB(), -(r_4)zv2->imagB());
+        zv1++;  zv2++;
+      }
+      /*     
       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());
       }
+      */
       nb_flop_ += (8.*(r_8)(vis.Size()-1));
     }