source: PSPA/Interface_Web/trunk/pspaWT/sources/controler/src/beam2Moments.cc @ 386

Last change on this file since 386 was 386, checked in by lemeur, 11 years ago

amélioration fiabilité graphique

File size: 8.0 KB
Line 
1#include "beam2Moments.h"
2#include <sstream>
3
4
5beam2Moments::beam2Moments() {
6    initializeMoments();
7}
8
9beam2Moments::~beam2Moments() {;}
10
11beam2Moments::beam2Moments( beam2Moments& tm) {
12    rij_secondOrderMoments_ = tm.rij_secondOrderMoments_;
13}
14
15beam2Moments::beam2Moments( const beam2Moments& tm) {
16    rij_secondOrderMoments_ = tm.rij_secondOrderMoments_;
17}
18
19beam2Moments::beam2Moments(double x, double xp, double y, double yp, double dl, double del) {
20    initializeMoments();
21    ( rij_secondOrderMoments_.at(0) ).at(0) = x;
22    ( rij_secondOrderMoments_.at(1) ).at(1) = xp;
23    ( rij_secondOrderMoments_.at(2) ).at(2) = y;
24    ( rij_secondOrderMoments_.at(3) ).at(3) = yp;
25    ( rij_secondOrderMoments_.at(4) ).at(4) = dl;
26    ( rij_secondOrderMoments_.at(5) ).at(5) = del;
27   
28}
29
30void beam2Moments::initializeMoments() {
31    rij_secondOrderMoments_.clear();
32    rij_secondOrderMoments_.resize(6);
33    unsigned dim=0;
34    unsigned k;
35    for ( k=0; k < 6; k++){
36        rij_secondOrderMoments_.at(k).resize(++dim, 0.0);
37    }
38}
39
40
41// const vector< vector<double> >& transportMoments::getMoments() const
42// {
43//   return rij_transportMoments_;
44// }
45
46vector< vector<double> >& beam2Moments::getMatrix()
47{
48    return rij_secondOrderMoments_;
49}
50
51
52void beam2Moments::impression() const {
53    unsigned j,m;
54    cout << " impression des moments " << endl;
55    for ( j = 0; j < 6; j++)
56    {
57        for (m=0; m <= j; m++)
58        {
59            cout  << ( rij_secondOrderMoments_.at(j) ).at(m) << " ";
60        }
61        cout << endl;
62    }
63}
64
65void beam2Moments::raz() {
66    // initialisation des moments
67    unsigned j,m;
68    for ( j = 0; j < rij_secondOrderMoments_.size(); j++)
69    {
70        for (m=0; m <= j; m++)
71        {
72            ( rij_secondOrderMoments_.at(j) ).at(m) = 0.0; // element r_jm
73        }
74    }
75}
76
77
78
79// suppose qu'on est bien positionné dans le stream (a revoir)
80bool beam2Moments::readFromTransportOutput(stringstream& inp) {
81  string bidString, bidString1, bidString2;
82  double bidon, bidon1, bidon2;
83   
84  // initialisation des moments
85  raz();
86  bool notFound = true;
87  string buf;
88  // sauter les eventuels 'VARY'
89
90  string varied = "VARIED:";
91  // while ( varied == "VARIED:" ) {
92  //   getline(inp, buf);
93  //   cout << " beam2Moments::readFromTransportOutput ligne lue : " << buf << endl;
94  //   stringstream str(buf);
95  //   str >> varied;
96  // }
97
98  string testVary;
99  while ( notFound ) {
100    getline(inp, buf);
101    //    cout << " beam2Moments::readFromTransportOutput ligne lue : " << buf << endl;
102    stringstream str(buf); 
103    if ( str >> bidon1 >>  bidString1 >>  bidon2 >>  ( rij_secondOrderMoments_.at(0) ).at(0) >> bidString2 )  {
104      // le format semble correct
105      if ( bidString1 != "M" ) {
106        // c'est bien ce qu'on attend : la longuer en metres
107        cout << " beam2Moments::readFromTransportOutput ON N'A PAS 'M' ligne lue : " << buf << endl;
108        notFound = true;
109        break;
110      } else {
111        // on a trouve la bonne ligne
112        cout << " beam2Moments::readFromTransportOutput lecture OK :  ligne lue : " << buf << endl;
113        notFound = false;
114      }
115    } else { 
116      // la ligne lue n'est pas au format
117      cout << "  beam2Moments::readFromTransportOutput format non conforme ligne lue : " << buf << endl;
118      stringstream test(buf);
119      if ( test >> testVary ) {
120        if ( testVary == varied ) {
121          // ca peut arriver... on continue
122          cout << "  beam2Moments::readFromTransportOutput VARIED RENCONTRE ligne lue : " << buf << endl;
123          notFound = true;
124        } else {
125          // on essaie de voir ce qui se passe (ce sera probablement a completer)
126          if ( testVary == "(" ) {
127          // on est trop loin, on est sur l'element suivant, on arrete
128          cout << "  beam2Moments::readFromTransportOutput ON EST TOMBE SUR ELEMNT SUIVANT ligne lue : " << buf << endl;
129            notFound = true;
130            break;
131          }
132          if ( testVary == "0" ) {
133            // on doit etre a la fin du fichier, ce qu'on verifie
134            if ( test >> bidString ) {
135              if ( bidString == "*LENGTH*" ) {
136                // c'est bien ca, on arrete
137                cout << "  beam2Moments::readFromTransportOutput ON EST A LA FIN DU FICHIER ligne lue : " << buf << endl;
138                notFound = true;
139                break;
140              } else {
141                // on n'y comprend plus grans chose, on arrete de toute facon
142                cout << "  beam2Moments::readFromTransportOutput SITUATION NON PREVUE 1 ligne lue : " << buf << endl;
143                notFound = true;
144                break;
145              }                // fin test lecture d'une chaine de caracteres
146            } else {
147              // on n'a pas pu lire la chaine, on n'y comprend rien, on arrete
148                cout << "  beam2Moments::readFromTransportOutput SITUATION NON PREVUE 2 ligne lue : " << buf << endl;
149                notFound = true;
150                break;
151            }
152          } else {
153            // le premier element lu n'est pas 0, on n'est pas a la fin
154            // mais il n'y a pas de faisceau a cet endroit, sortie normale
155                cout << "  beam2Moments::readFromTransportOutput SORTIE NORMALE ligne lue : " << buf << endl;
156            notFound = true;
157            break;
158          }
159          // si on arrive la, on est en pleine panade
160                cout << "  beam2Moments::readFromTransportOutput SITUATION NON PREVUE 4 ligne lue : " << buf << endl;
161            notFound = true;
162            break;
163        } 
164      } else {
165        // on n'a pas pu lire le premier element de la ligne (bizarre)
166                cout << "  beam2Moments::readFromTransportOutput SITUATION NON PREVUE 5 ligne lue : " << buf << endl;
167        notFound = true;
168        break;
169      }
170    }
171  } // fin de boucle sur les lignes
172
173  if ( notFound ) return false;
174
175
176  if ( ! ( inp >> bidon >> ( rij_secondOrderMoments_.at(1) ).at(1) >>  bidString >> ( rij_secondOrderMoments_.at(1) ).at(0) ) ) return false;
177  if ( ! (inp >> bidon >> ( rij_secondOrderMoments_.at(2) ).at(2) >>  bidString >> ( rij_secondOrderMoments_.at(2) ).at(0)  >> ( rij_secondOrderMoments_.at(2) ).at(1) ) ) return false;
178  if ( ! (inp >> bidon >> ( rij_secondOrderMoments_.at(3) ).at(3) >>  bidString >> ( rij_secondOrderMoments_.at(3) ).at(0)  >> ( rij_secondOrderMoments_.at(3) ).at(1) >> ( rij_secondOrderMoments_.at(3) ).at(2) ) ) return false; 
179   
180  if ( ! ( inp >> bidon >> ( rij_secondOrderMoments_.at(4) ).at(4) >>  bidString >> ( rij_secondOrderMoments_.at(4) ).at(0)  >> ( rij_secondOrderMoments_.at(4) ).at(1) >> ( rij_secondOrderMoments_.at(4) ).at(2) >> ( rij_secondOrderMoments_.at(4) ).at(3) ) ) return false;
181   
182  if ( ! (inp >> bidon >> ( rij_secondOrderMoments_.at(5) ).at(5) >>  bidString >> ( rij_secondOrderMoments_.at(5) ).at(0)  >> ( rij_secondOrderMoments_.at(5) ).at(1) >> ( rij_secondOrderMoments_.at(5) ).at(2) >> ( rij_secondOrderMoments_.at(5) ).at(3) >> ( rij_secondOrderMoments_.at(5) ).at(4) ) ) return false;
183
184  return true;   
185}
186
187
188string beam2Moments::writeToTransportInput(double P0) const {
189    ostringstream sortie;
190    sortie << " BEAM, X=" << (rij_secondOrderMoments_.at(0)).at(0) << ", XP="<< (rij_secondOrderMoments_.at(1)).at(1) << ", Y=" << (rij_secondOrderMoments_.at(2)).at(2)<< ", YP=" << (rij_secondOrderMoments_.at(3)).at(3) << ", & " << endl;
191   
192    sortie << " DL=" << (rij_secondOrderMoments_.at(4)).at(4) << ", DEL="<< (rij_secondOrderMoments_.at(5)).at(5) << ", P0=" << P0 << " ; " << endl;
193   
194    sortie << " CORR, C21=" << (rij_secondOrderMoments_.at(1)).at(0) << ", C31="<< (rij_secondOrderMoments_.at(2)).at(0) << ", C32=" << (rij_secondOrderMoments_.at(2)).at(1) << ", & " << endl;
195   
196    sortie << " C41=" << (rij_secondOrderMoments_.at(3)).at(0) << ", C42="<< (rij_secondOrderMoments_.at(3)).at(1) << ", C43=" << (rij_secondOrderMoments_.at(3)).at(2) << ", C51=" << (rij_secondOrderMoments_.at(4)).at(0) << ", & " << endl;
197   
198    sortie << " C52=" << (rij_secondOrderMoments_.at(4)).at(1) << ", C53="<< (rij_secondOrderMoments_.at(4)).at(2) << ", C54=" << (rij_secondOrderMoments_.at(4)).at(3) << ", C61=" << (rij_secondOrderMoments_.at(5)).at(0) << ", & " << endl;
199   
200    sortie << " C62=" << (rij_secondOrderMoments_.at(5)).at(1) << ", C63="<< (rij_secondOrderMoments_.at(5)).at(2) << ", C64=" << (rij_secondOrderMoments_.at(5)).at(3) << ", C65=" << (rij_secondOrderMoments_.at(5)).at(4) << " ; " << endl;
201    return sortie.str();
202}
203
204beam2Moments& beam2Moments::operator= (const beam2Moments& tpm) {
205    rij_secondOrderMoments_ = tpm.rij_secondOrderMoments_;
206    return *this;
207}
Note: See TracBrowser for help on using the repository browser.