Changeset 3918 in Sophya for trunk/AddOn/TAcq/racquproc.cc

Timestamp:

Dec 2, 2010, 4:41:13 PM (15 years ago)

Author:

ansari

Message:

Amelioration code thread monitoring pour tenir compte du mode fftshrink + calcul matrice saturation, cmv+Reza

File:

• trunk/AddOn/TAcq/racquproc.cc (modified) (7 diffs)

trunk/AddOn/TAcq/racquproc.cc

@@ -177 +177 @@
 
   // Histo pour valeurs des bytes dans les paquets 
-  Histo* phist[2*MAXNBFIB] = {NULL,NULL,NULL,NULL, NULL,NULL,NULL,NULL};
-  Histo2D* ph2[2*MAXNBFIB] = {NULL,NULL,NULL,NULL, NULL,NULL,NULL,NULL};
-  Histo* phexp[2*MAXNBFIB] = {NULL,NULL,NULL,NULL, NULL,NULL,NULL,NULL};
+  Histo* phist[2*MAXNBFIB];
+  Histo2D* ph2[2*MAXNBFIB];
+  Histo* phexp[2*MAXNBFIB];
+  TMatrix<r_4>* psatur[2*MAXNBFIB];  
+  for(int ifib=0; ifib<2*MAXNBFIB; ifib++) {
+    phist[ifib]=NULL; ph2[ifib]=NULL; phexp[ifib]=NULL; psatur[ifib]=NULL;
+  }
   if (par_.fgdatafft) {    // data de type FFT   
     for(sa_size_t lc=0; lc<2*memgr.NbFibres(); lc++) {
@@ -185 +189 @@
       ph2[lc]   = new Histo2D(-128.5,+128.5,257, -128.5,+128.5,257);
       phexp[lc] = new Histo(-128.5,+128.5,257);
+      psatur[lc] = new TMatrix<r_4>(3,3);
+      *psatur[lc] = 0.;
     }
   }
@@ -218 +224 @@
       }
     }
+
+    // DBG-ATTENTION : freq 0 pour expo full FFT si reordonne 
+    int idxexpo=(par_.fgfftshrink)?1023:0; // exposant FFT shrink 1024 freq, ou freq 0 pour full-FFT
     
     cpaqok_ = (uint_1)0;
@@ -230 +239 @@
         if (par_.fgdatafft) {    // Traitement data de type FFT 
           TwoByteComplex* tbcp=paq.Data1C();
-          phexp[lc]->Add((r_4)tbcp[0].realD());
+
+          //  exposant FFT 
+          phexp[lc]->Add((r_4)tbcp[idxexpo].realD());
+
           for(sa_size_t j=1; j<spectre_.NCols(); j++) {
             phist[lc]->Add((r_4)tbcp[j].realD());
             phist[lc]->Add((r_4)tbcp[j].imagD());
             ph2[lc]->Add((r_4)tbcp[j].realD(), (r_4)tbcp[j].imagD());
+            // comptage saturation 
+            int sir=0, sii=0;  // saturation partie reelle, partie imaginaire
+            if (tbcp[j].realI()==-128) sir=-1;
+            else if (tbcp[j].realI()==127) sir=1;
+            if (tbcp[j].imagI()==-128) sii=-1;
+            else if (tbcp[j].imagI()==127) sii=1;
+            (*psatur[lc])(sir+1,sii+1)+=1.;
+
             spectre_(lc,j) += Zmod2(tbcp[j]);
           }
           nzm_[lc]++;   
           tbcp=paq.Data2C();   lc++;
-          phexp[lc]->Add((r_4)tbcp[0].realD());
+          phexp[lc]->Add((r_4)tbcp[idxexpo].realD());
           for(sa_size_t j=1; j<spectre_.NCols(); j++) {
             phist[lc]->Add((r_4)tbcp[j].realD());
             phist[lc]->Add((r_4)tbcp[j].imagD());
             ph2[lc]->Add((r_4)tbcp[j].realD(), (r_4)tbcp[j].imagD());
+            // comptage saturation 
+            int sir=0, sii=0;  // saturation partie reelle, partie imaginaire
+            if (tbcp[j].realI()==-128) sir=-1;
+            else if (tbcp[j].realI()==127) sir=1;
+            if (tbcp[j].imagI()==-128) sii=-1;
+            else if (tbcp[j].imagI()==127) sii=1;
+            (*psatur[lc])(sir+1,sii+1)+=1.;
+
             spectre_(lc,j) += Zmod2(tbcp[j]);
           }
@@ -303 +331 @@
           po << PPFNameTag(buftag) << (*phexp[lc]);
           phexp[lc]->Zero();
+          sprintf(buftag,"saturV%d",(int)lc);
+          po << PPFNameTag(buftag) << (*psatur[lc]);
+          cout << " Matrice Saturation Voie[" << lc << "] : " <<  (*psatur[lc]);
+          if (nzm_[lc]>0)  { 
+            (*psatur[lc]) /= ((r_4)nzm_[lc]*0.01);  
+            cout << " Matrice Saturation Pourcentage Voie[" << lc << "] : " <<  (*psatur[lc]);
+          }
+          (*psatur[lc])=0.;
         }
         nzm_[lc]=0;
@@ -351 +387 @@
         po << PPFNameTag(buftag) << (*phexp[lc]);
         delete phexp[lc];
+        sprintf(buftag,"saturV%d",(int)lc);
+        po << PPFNameTag(buftag) << (*psatur[lc]);
+        cout << " Matrice Saturation Voie[" << lc << "] : " <<  (*psatur[lc]);
+        if (nzm_[lc]>0)  { 
+          (*psatur[lc]) /= ((r_4)nzm_[lc]*0.01);  
+          cout << " Matrice Saturation Pourcentage Voie[" << lc << "] : " <<  (*psatur[lc]);
+        }
+        (*psatur[lc])=0.;
       }
       nzm_[lc]=0;
@@ -369 +413 @@
          << " -> " << fracsame << " % )" << endl;
   }
+
+  for(sa_size_t lc=0; lc<2*memgr.NbFibres(); lc++) {
+    if (phist[lc]) delete phist[lc];
+    if (ph2[lc]) delete ph2[lc]; 
+    if (phexp[lc]) delete phexp[lc]; 
+    if (psatur[lc]) delete psatur[lc]; 
+  }
+
   return 0;
 }