[236] | 1 | #include "beam2Moments.h" |
---|
| 2 | #include <sstream> |
---|
| 3 | |
---|
| 4 | |
---|
| 5 | beam2Moments::beam2Moments() { |
---|
[348] | 6 | initializeMoments(); |
---|
[236] | 7 | } |
---|
| 8 | |
---|
| 9 | beam2Moments::~beam2Moments() {;} |
---|
| 10 | |
---|
| 11 | beam2Moments::beam2Moments( beam2Moments& tm) { |
---|
[348] | 12 | rij_secondOrderMoments_ = tm.rij_secondOrderMoments_; |
---|
[236] | 13 | } |
---|
| 14 | |
---|
| 15 | beam2Moments::beam2Moments( const beam2Moments& tm) { |
---|
[348] | 16 | rij_secondOrderMoments_ = tm.rij_secondOrderMoments_; |
---|
[236] | 17 | } |
---|
| 18 | |
---|
| 19 | beam2Moments::beam2Moments(double x, double xp, double y, double yp, double dl, double del) { |
---|
[348] | 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 | |
---|
[236] | 28 | } |
---|
| 29 | |
---|
| 30 | void beam2Moments::initializeMoments() { |
---|
[348] | 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 | } |
---|
[236] | 38 | } |
---|
| 39 | |
---|
| 40 | |
---|
[348] | 41 | // const vector< vector<double> >& transportMoments::getMoments() const |
---|
| 42 | // { |
---|
| 43 | // return rij_transportMoments_; |
---|
| 44 | // } |
---|
[236] | 45 | |
---|
[348] | 46 | vector< vector<double> >& beam2Moments::getMatrix() |
---|
| 47 | { |
---|
| 48 | return rij_secondOrderMoments_; |
---|
| 49 | } |
---|
[236] | 50 | |
---|
| 51 | |
---|
| 52 | void beam2Moments::impression() const { |
---|
[348] | 53 | unsigned j,m; |
---|
| 54 | cout << " impression des moments " << endl; |
---|
| 55 | for ( j = 0; j < 6; j++) |
---|
[236] | 56 | { |
---|
[348] | 57 | for (m=0; m <= j; m++) |
---|
| 58 | { |
---|
| 59 | cout << ( rij_secondOrderMoments_.at(j) ).at(m) << " "; |
---|
| 60 | } |
---|
| 61 | cout << endl; |
---|
[236] | 62 | } |
---|
| 63 | } |
---|
| 64 | |
---|
| 65 | void beam2Moments::raz() { |
---|
[348] | 66 | // initialisation des moments |
---|
| 67 | unsigned j,m; |
---|
| 68 | for ( j = 0; j < rij_secondOrderMoments_.size(); j++) |
---|
[236] | 69 | { |
---|
[348] | 70 | for (m=0; m <= j; m++) |
---|
| 71 | { |
---|
| 72 | ( rij_secondOrderMoments_.at(j) ).at(m) = 0.0; // element r_jm |
---|
| 73 | } |
---|
[236] | 74 | } |
---|
| 75 | } |
---|
| 76 | |
---|
| 77 | |
---|
| 78 | |
---|
| 79 | // suppose qu'on est bien positionné dans le stream (a revoir) |
---|
[386] | 80 | bool beam2Moments::readFromTransportOutput(stringstream& inp) { |
---|
[355] | 81 | string bidString, bidString1, bidString2; |
---|
| 82 | double bidon, bidon1, bidon2; |
---|
[348] | 83 | |
---|
[355] | 84 | // initialisation des moments |
---|
| 85 | raz(); |
---|
| 86 | bool notFound = true; |
---|
| 87 | string buf; |
---|
| 88 | // sauter les eventuels 'VARY' |
---|
[386] | 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; |
---|
[355] | 99 | while ( notFound ) { |
---|
| 100 | getline(inp, buf); |
---|
[386] | 101 | // cout << " beam2Moments::readFromTransportOutput ligne lue : " << buf << endl; |
---|
| 102 | stringstream str(buf); |
---|
[355] | 103 | if ( str >> bidon1 >> bidString1 >> bidon2 >> ( rij_secondOrderMoments_.at(0) ).at(0) >> bidString2 ) { |
---|
[386] | 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 | } |
---|
[348] | 170 | } |
---|
[386] | 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; |
---|
[348] | 179 | |
---|
[386] | 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; |
---|
[348] | 181 | |
---|
[386] | 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; |
---|
[236] | 185 | } |
---|
| 186 | |
---|
| 187 | |
---|
| 188 | string beam2Moments::writeToTransportInput(double P0) const { |
---|
[348] | 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; |
---|
[236] | 201 | return sortie.str(); |
---|
| 202 | } |
---|
| 203 | |
---|
| 204 | beam2Moments& beam2Moments::operator= (const beam2Moments& tpm) { |
---|
[348] | 205 | rij_secondOrderMoments_ = tpm.rij_secondOrderMoments_; |
---|
| 206 | return *this; |
---|
[236] | 207 | } |
---|