Changeset 3956 in Sophya for trunk/AddOn/TAcq/vismfib.cc

Timestamp:

Mar 2, 2011, 3:31:39 PM (15 years ago)

Author:

ansari

Message:

Amelioration du processeur de calcul de visibilite (BRVisibilityCalculator) et du programme vismfib.cc pour permettre la prise en charge des donnees raw-2c pour le calcul des visibilites et ajout de la possibilite d ecrire les fichiers de sortie (matrices de visibilites) au format FITS, Reza 02/03/2011

File:

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

trunk/AddOn/TAcq/vismfib.cc

@@ -58 +58 @@
     ResourceUsage resu;
 
-    cout << " visimfib: Creating MemZoneMgr/threads - PaqSz= " << par.paqsize_ << endl;
+    BRPaquet paq(par.paqsize_);
+    uint_4 procsz=sizeof(float)*(paq.DataSize()+4);
+    if ((par.fgdatafft_)||(par.action_=="cube3d")) procsz = 0;
+    cout << " visimfib: Creating MemZoneMgr/processing threads - PaqSz= " << par.paqsize_ 
+         << " ProcPaqSz=" << procsz << endl;
 
-    RAcqMemZoneMgr mmgr(par.nzones_, par.dirlist_.size(), par.npaqinzone_, par.paqsize_);
-    if (par.action_ == "cube3d")  mmgr.SetFinalizedMask((uint_4)MemZS_Saved);
-    else  {
-      if (par.nbcalgrp_==2) 
-        mmgr.SetFinalizedMask((uint_4)MemZS_ProcA|(uint_4)MemZS_ProcB);
-      else if (par.nbcalgrp_==3) 
-        mmgr.SetFinalizedMask((uint_4)MemZS_ProcA|(uint_4)MemZS_ProcB|(uint_4)MemZS_ProcC);
-      else if (par.nbcalgrp_==4) 
-        mmgr.SetFinalizedMask((uint_4)MemZS_ProcA|(uint_4)MemZS_ProcB|(uint_4)MemZS_ProcC|(uint_4)MemZS_ProcD);
-      else if (par.nbcalgrp_==5) 
-        mmgr.SetFinalizedMask((uint_4)MemZS_ProcA|(uint_4)MemZS_ProcB|(uint_4)MemZS_ProcC|(uint_4)MemZS_ProcD|(uint_4)MemZS_ProcE);
-      else if (par.nbcalgrp_==6) 
-        mmgr.SetFinalizedMask((uint_4)MemZS_ProcA|(uint_4)MemZS_ProcB|(uint_4)MemZS_ProcC|
-                              (uint_4)MemZS_ProcD|(uint_4)MemZS_ProcE|(uint_4)MemZS_ProcF);
-      else mmgr.SetFinalizedMask((uint_4)MemZS_ProcA);  
-    }
+    RAcqMemZoneMgr mmgr(par.nzones_, par.dirlist_.size(), par.npaqinzone_, par.paqsize_, procsz);
+    
     BRMultiFitsReader reader(mmgr, par.dirlist_, par.rdsamefc_, par.imin_, par.imax_, par.istep_);
     reader.SetPrintLevel(par.prtlevel_,par.prtmodulo_);
@@ -85 +75 @@
 
     if (par.action_ == "viscktt")  procg.ActivateTimeTagCheck(par.TotalNPaquets());
+    if (par.fgfitsout_) procg.SetFitsOutput();
+
+    BRFFTCalculator procfft(mmgr, par.fgsinglechannel_);
+
+    // On determine MemZaction pour chaque processeur et le finalizemask en fonction du traitement demande 
+    if (par.action_ == "cube3d")  mmgr.SetFinalizedMask((uint_4)MemZS_Saved);
+    else {
+      MemZStatus mfmask=MemZS_ProcA;
+      if (par.fgdatafft_) mfmask=procg.SetMemZAction(MemZA_ProcA);
+      else mfmask=procg.SetMemZAction(MemZA_ProcB);
+      mmgr.SetFinalizedMask((uint_4)mfmask);
+    }
 
     FitsCubeWriter wrt(mmgr, par.outpath_, par.nbloc_);
@@ -92 +94 @@
     reader.start();
     if (par.action_ == "cube3d")   wrt.start();
-    else  procg.start();
+    else  { 
+      if (!par.fgdatafft_)  procfft.start(); 
+      procg.start();
+    }
     usleep(200000);
     reader.join();
     if (par.action_ == "cube3d")  wrt.join();
-    else procg.join();
+    else {
+      if (!par.fgdatafft_)  procfft.join(); 
+      procg.join();
+    }
     mmgr.Print(cout);
     cout << resu ;