Changeset 342 in PSPA for Interface_Web/trunk/pspaWT/sources/controler/src/softwareParmela.cc

Timestamp:

Feb 26, 2013, 2:35:45 PM (11 years ago)

Author:

touze

Message:

nvx element snapshot

File:

• Interface_Web/trunk/pspaWT/sources/controler/src/softwareParmela.cc (modified) (10 diffs)

Interface_Web/trunk/pspaWT/sources/controler/src/softwareParmela.cc

@@ -8 +8 @@
 
 softwareParmela::softwareParmela() : abstractSoftware()
-{
-  ;
-}
+{;}
 
 softwareParmela::softwareParmela(string inputFileName, globalParameters* globals, dataManager* dt) : abstractSoftware(inputFileName, globals, dt)
-{
-  ;
-}
+{;}
 
 bool softwareParmela::createInputFile( particleBeam* beamBefore, unsigned int numeroDeb, unsigned int numeroFin, string workingDir)
@@ -26 +22 @@
     return false;
   }
-
 
   ofstream outfile;
@@ -54 +49 @@
   }
 
-
+  
   outfile << "TITLE" << endl;
   outfile << " titre provisoire " << endl;
@@ -66 +61 @@
     outfile << "INPUT 0 /NP=" << beamBefore->getNbParticles() << endl;
   }
+
+  
   for ( k = premier; k <= numeroFin; k++)
     {
       elPtr = dataManager_->getElementPointerFromNumero(k);
+      if(elPtr->getNomdElement().getElementType() == snapshot) continue;
       outfile << elPtr->parmelaOutputFlow() << endl;
     }
-
+  
   outfile << "ZOUT" << endl;
   outfile << "START /wt=0.0 /dwt=" << globParamPtr_->getIntegrationStep() << "  /nsteps=" << globParamPtr_->getNbSteps() << "  nsc=" << globParamPtr_->getScPeriod() << "  /nout=10" << endl;
@@ -81 +79 @@
 
 
-bool  softwareParmela::execute(unsigned int numeroDeb,unsigned int numeroFin,string workingDir)
-{
+bool softwareParmela::execute(unsigned int numeroDeb,unsigned int numeroFin,string workingDir)
+{
+  bool ExecuteStatus = true;
+  
   ostringstream sortie;
-  bool ExecuteStatus = true;
-  //  resul.clear();
   sortie << " EXECUTION DE PARMELA DE l'ELEMENT " << numeroDeb << " A L'ELEMENT " << numeroFin << endl;
 
-  char buf[132];
   string parmelaJob = workingDir + "parmela";
   parmelaJob += string("   ");
   parmelaJob += workingDir;
-  //  cout << " job parmela= " << parmelaJob << endl;
-
+  
   string resultOfRun;
-  bool success = launchJob(parmelaJob, resultOfRun);
+  bool success = launchJob(parmelaJob,resultOfRun);
   sortie << resultOfRun << endl;
-  if ( !success) {
+  if ( !success ) {
     sortie << " launching of parmela failed " << endl;
     ExecuteStatus = false;
-  }
-  else {
+  } else {
     sortie << " successful launching of parmela " << endl;
     cout << " execution parmela MARCHE " << endl;
@@ -112 +107 @@
   }
 
-  //  resul =  sortie.str();
-    dataManager_->consoleMessage(sortie.str());
-  
+  dataManager_->consoleMessage(sortie.str());
   return ExecuteStatus;
 }
 
-
-bool  softwareParmela::buildBeamAfterElements(unsigned int numeroDeb,unsigned int numeroFin, vector<particleBeam>& beams, string workingDir) {
+bool softwareParmela::buildBeamAfterElements(unsigned int numeroDeb,unsigned int numeroFin, vector<particleBeam>& beams, string workingDir) 
+{
   bool result = true;
-  unsigned k;
-
-  for ( k= numeroDeb; k <= numeroFin; k++)
-    {
+  cout << "debug:: debut " << numeroDeb << ", fin " << numeroFin << endl;
+
+  // index du premier element de parmela
+  int id= numeroDeb-offsetNumElem_;
+
+  for(unsigned k = numeroDeb; k <= numeroFin; k++) 
+    {
+      abstractElement* elem = dataManager_->getElementPointerFromNumero(k);
+
+      // si l'element est un snapshot, recuperer la sortie precedente
+      if(elem->getNomdElement().getElementType() == snapshot) {
+        int avantDernier = beams.size() - 1;
+        cout<<"["<<k<<"] : ecrit sur fichier le contenu de beam["<<avantDernier<<"]"<<endl;
+        continue;
+      }
+
+      // sinon c'est un element de parmela d'index id
       beams.push_back(particleBeam());
       vector<double> centroid;
       bareParticle refPart;
       vector<bareParticle> particles;
-      if (!beamFromParmela(workingDir,k - offsetNumElem_ , globParamPtr_->getFrequency(), centroid, refPart,particles ))
+      double freq= globParamPtr_->getFrequency();
+
+      if(!beamFromParmela(workingDir,id,freq,centroid,refPart,particles)) 
         {
-          abstractElement* elem = dataManager_->getElementPointerFromNumero(k);
-          if ( elem->is_accepted_by_software(nomDeLogiciel::parmela) == warning) {
+          if(elem->is_accepted_by_software(nomDeLogiciel::parmela) == warning) {
             int avantDernier = beams.size() -2;
             beams.back() = beams.at(avantDernier);
           } else {
-            // sortie << " reading parmdesz  failed " << endl;
-            dataManager_->consoleMessage(" softwareParmela::buildBeamAfterElements : reading parmdesz  failed " );
+            dataManager_->consoleMessage("softwareParmela::buildBeamAfterElements : reading parmdesz failed");
             result = false;
             break;
           }
-        }
-      else {
-        beams.back().setWithParticles(centroid, refPart,particles); 
-      }
+        } else {
+        beams.back().setWithParticles(centroid,refPart,particles); 
+      }
+
+      // l'element de parmela suivant
+      id++;
     }
+
   return result;
 }
 
-
-
-
-bool softwareParmela::beamFromParmela(string workingDir,unsigned numeroParmel, double referencefrequency, vector<double>& centroid, bareParticle& refPart,vector<bareParticle>& particles ) {
-  unsigned  k;
-  FILE* filefais;
+bool softwareParmela::beamFromParmela(string workingDir,unsigned numeroParmel, double referencefrequency, vector<double>& centroid, bareParticle& refPart,vector<bareParticle>& particles ) 
+{
   string nomfilefais = workingDir + "parmdesz";
   cout << " nom fichier desz : " << nomfilefais << endl;
-  filefais = fopen(nomfilefais.c_str(), "r");
+  FILE *filefais = fopen(nomfilefais.c_str(), "r");
 
   if ( filefais == (FILE*)0 ) {
     dataManager_->consoleMessage(" beamFromParmela() erreur a l'ouverture du fichier 'parmdesz'");
-    cerr << " beamFromParmela() erreur a l'ouverture du fichier" << nomfilefais  << endl;; 
-    return false;
-  }
-  else cout << " beamFromParmela() : ouverture du fichier " << nomfilefais << endl; 
-
+    cerr << " beamFromParmela() erreur a l'ouverture du fichier" << nomfilefais  << endl; 
+    return false;
+  } else 
+    cout << " beamFromParmela() : ouverture du fichier " << nomfilefais << endl; 
+  
   parmelaParticle partic;
   std::vector<parmelaParticle> faisceau;
-
-  //  unsigned indexElement = numeroElement-1;
- 
-
-
-
   int testNombrePartRef =0;
   double phaseRef = 0.0;
 
-  while( partic.readFromParmelaFile(filefais) > 0 ) {
-    if ( partic.ne == (int)numeroParmel )
-      {
+  while( partic.readFromParmelaFile(filefais) > 0 ) 
+    {
+
+      if ( partic.ne == (int)numeroParmel ) {
         faisceau.push_back(partic);
-
         if ( partic.np == 1 ) {
           // en principe on est sur la particule de reference
@@ -200 +199 @@
         }
       }
-  }
+    } //while
 
   if ( faisceau.size() == 0) 
@@ -213 +212 @@
   // avec freq en Mhz et dphi en degres et résultat en cm: 
   double FACTEUR =  83.3333;  // ameliorer la precision
-
-
 
   // pour l'instant on choisit un centroid nul;
@@ -232 +229 @@
   // contrairement a ce qu'indique la notice PARMELA, dans parmdesz, les xp et yp 
   // sont donnes en radians
-  for ( k=0; k < faisceau.size(); k++) {
-    x=faisceau.at(k).xx;
+  for (unsigned k = 0; k < faisceau.size(); k++) {
+    x= faisceau.at(k).xx;
     xp=faisceau.at(k).xxp;
-    y=faisceau.at(k).yy;
+    y= faisceau.at(k).yy;
     yp=faisceau.at(k).yyp;
-
     // dephasage par rapport a la reference  
     dephas = faisceau.at(k).phi - phaseRef; // degrés
@@ -250 +246 @@
     particles.at(k) = bareParticle(pos,betagamma);
   }
+  
   return true;
 }