| 1 | // ArchTOIPipe (C) CEA/DAPNIA/SPP IN2P3/LAL
|
|---|
| 2 | // Eric Aubourg
|
|---|
| 3 | // Christophe Magneville
|
|---|
| 4 | // Reza Ansari
|
|---|
| 5 | // $Id: toimanager.cc,v 1.24 2003-10-13 20:48:37 aubourg Exp $
|
|---|
| 6 |
|
|---|
| 7 | #include "toimanager.h"
|
|---|
| 8 | #include <limits.h>
|
|---|
| 9 | #include <pthread.h>
|
|---|
| 10 | #include <iostream.h>
|
|---|
| 11 | #include <unistd.h>
|
|---|
| 12 | #include <map>
|
|---|
| 13 |
|
|---|
| 14 | //#ifndef NOCFITSIO
|
|---|
| 15 | #include "fitstoiwtr.h"
|
|---|
| 16 | //#endif
|
|---|
| 17 |
|
|---|
| 18 | #ifndef MAXINT
|
|---|
| 19 | #define MAXINT 2147483647
|
|---|
| 20 | #endif
|
|---|
| 21 |
|
|---|
| 22 | TOIManager::TOIManager() {
|
|---|
| 23 | reqBegin = 0;
|
|---|
| 24 | reqEnd = MAXINT;
|
|---|
| 25 |
|
|---|
| 26 | // -----------ajout cgt vf 19/08/2002
|
|---|
| 27 | // par defaut TOISegmented
|
|---|
| 28 | selectTOISegmented(1024, 20);
|
|---|
| 29 | // ----------- fin ajout cgt
|
|---|
| 30 |
|
|---|
| 31 | }
|
|---|
| 32 |
|
|---|
| 33 | TOIManager* TOIManager::instance = NULL;
|
|---|
| 34 |
|
|---|
| 35 | TOIManager* TOIManager::getManager() {
|
|---|
| 36 | if (instance == NULL) instance = new TOIManager();
|
|---|
| 37 | return instance;
|
|---|
| 38 | }
|
|---|
| 39 |
|
|---|
| 40 | // ajout vf 26/07/2002
|
|---|
| 41 |
|
|---|
| 42 | // enregistrement d'un processeur dans la liste des processeurs pour une execution en groupe
|
|---|
| 43 |
|
|---|
| 44 | void TOIManager::registerProcessor(TOIProcessor* proc) {
|
|---|
| 45 |
|
|---|
| 46 | cout << "Adding processor to TOIManager for group execution" << endl;
|
|---|
| 47 | processors.push_back(proc);
|
|---|
| 48 |
|
|---|
| 49 | }
|
|---|
| 50 |
|
|---|
| 51 | void TOIManager::registerProcessor(RingProcessor* proc) {
|
|---|
| 52 |
|
|---|
| 53 | cout << "Adding processor to TOIManager for group execution" << endl;
|
|---|
| 54 | ringProcessors.push_back(proc);
|
|---|
| 55 |
|
|---|
| 56 | }
|
|---|
| 57 |
|
|---|
| 58 |
|
|---|
| 59 | // demarrage de tous les processeurs et verification auto des samplenum pour chaque processeur parametre
|
|---|
| 60 |
|
|---|
| 61 | void TOIManager::startAll() {
|
|---|
| 62 | // verification des samplenum
|
|---|
| 63 | bool samples_ok=checkSamplesLimits(1);
|
|---|
| 64 | if (samples_ok) {
|
|---|
| 65 | cout << "All limits ok" << endl << "Starting processors" << endl;
|
|---|
| 66 | } else {
|
|---|
| 67 | cout << "One or more limits ajusted for execution" << endl << "Starting processors" << endl;
|
|---|
| 68 | }
|
|---|
| 69 |
|
|---|
| 70 | // mise a jour des limites apres verification
|
|---|
| 71 | checkSamplesLimits(2);
|
|---|
| 72 | checkSamplesLimits(3);
|
|---|
| 73 | cout<<"Fin checks"<<endl;
|
|---|
| 74 |
|
|---|
| 75 |
|
|---|
| 76 | // debogage affichage des limites apres calcul
|
|---|
| 77 | {for (vector<TOIProcessor*>::iterator i = processors.begin();
|
|---|
| 78 | i != processors.end(); i++) {
|
|---|
| 79 | TOIProcessor* proc = *i;
|
|---|
| 80 | proc->printLimits();
|
|---|
| 81 | }
|
|---|
| 82 | }
|
|---|
| 83 |
|
|---|
| 84 | // demarrage
|
|---|
| 85 | {for (vector<TOIProcessor*>::iterator i = processors.begin();
|
|---|
| 86 | i != processors.end(); i++) {
|
|---|
| 87 | TOIProcessor* proc = *i;
|
|---|
| 88 | cout << "**********************" << endl;
|
|---|
| 89 | cout << "starting processor " << endl;
|
|---|
| 90 | proc->start();
|
|---|
| 91 | cout << "processor started " << endl;
|
|---|
| 92 | }}
|
|---|
| 93 |
|
|---|
| 94 | {for (vector<RingProcessor*>::iterator i = ringProcessors.begin();
|
|---|
| 95 | i != ringProcessors.end(); i++) {
|
|---|
| 96 | RingProcessor* proc = *i;
|
|---|
| 97 | cout << "**********************" << endl;
|
|---|
| 98 | cout << "starting processor " << endl;
|
|---|
| 99 | proc->start();
|
|---|
| 100 | cout << "processor started " << endl;
|
|---|
| 101 | }}
|
|---|
| 102 | cout << "**********************" << endl;
|
|---|
| 103 | }
|
|---|
| 104 |
|
|---|
| 105 | void TOIManager::checkRingLimits()
|
|---|
| 106 | {
|
|---|
| 107 | int min=0;
|
|---|
| 108 | int max = MAXINT;
|
|---|
| 109 | for (int pass=0; pass<2; pass++) {
|
|---|
| 110 | for (vector<RingProcessor*>::iterator i = ringProcessors.begin();
|
|---|
| 111 | i != ringProcessors.end(); i++) {
|
|---|
| 112 | RingProcessor* proc = *i;
|
|---|
| 113 | proc->getRingRange(min, max);
|
|---|
| 114 | }
|
|---|
| 115 | }
|
|---|
| 116 | }
|
|---|
| 117 |
|
|---|
| 118 | bool TOIManager::checkSamplesLimits(int pass)
|
|---|
| 119 | {
|
|---|
| 120 | bool processor_ok=true;
|
|---|
| 121 | bool samples_ok=true;
|
|---|
| 122 | for (vector<TOIProcessor*>::iterator i = processors.begin();
|
|---|
| 123 | i != processors.end(); i++) {
|
|---|
| 124 | TOIProcessor* proc = *i;
|
|---|
| 125 | cout << "testing processor limits " << endl;
|
|---|
| 126 | // test du processeur
|
|---|
| 127 |
|
|---|
| 128 | // test seulement pour les processor cle
|
|---|
| 129 | //if (proc->getRequested()) {
|
|---|
| 130 | processor_ok = proc->checkSampleLimits(pass);
|
|---|
| 131 | //}
|
|---|
| 132 |
|
|---|
| 133 | if (processor_ok) {
|
|---|
| 134 | cout << "processor limits ok " << endl;
|
|---|
| 135 | } else {
|
|---|
| 136 | cout << "processor limits ajusted" << endl;
|
|---|
| 137 | samples_ok = false;
|
|---|
| 138 | }
|
|---|
| 139 | }
|
|---|
| 140 | return samples_ok;
|
|---|
| 141 | }
|
|---|
| 142 |
|
|---|
| 143 | // fin ajout vf
|
|---|
| 144 |
|
|---|
| 145 | void TOIManager::setRequestedSample(long begin, long end) {
|
|---|
| 146 | cout << "TOIManager::setRequestedSample should not be called anymore" << endl;
|
|---|
| 147 | cout << " call same method one any toiprocessor" << endl;
|
|---|
| 148 | exit(-1);
|
|---|
| 149 | reqBegin = begin;
|
|---|
| 150 | reqEnd = end;
|
|---|
| 151 | }
|
|---|
| 152 |
|
|---|
| 153 | int TOIManager::getRequestedBegin() {
|
|---|
| 154 | return reqBegin;
|
|---|
| 155 | }
|
|---|
| 156 |
|
|---|
| 157 | int TOIManager::getRequestedEnd() {
|
|---|
| 158 | return reqEnd;
|
|---|
| 159 | }
|
|---|
| 160 |
|
|---|
| 161 | void TOIManager::addThread(pthread_t* t) {
|
|---|
| 162 | // cout << "adding thread " << t << endl;
|
|---|
| 163 | threads.push_back(t);
|
|---|
| 164 | }
|
|---|
| 165 |
|
|---|
| 166 | void TOIManager::joinAll() {
|
|---|
| 167 | waitForAll();
|
|---|
| 168 | }
|
|---|
| 169 |
|
|---|
| 170 | void TOIManager::waitForAll() {
|
|---|
| 171 | for (vector<pthread_t*>::iterator i = threads.begin();
|
|---|
| 172 | i != threads.end(); i++) {
|
|---|
| 173 | pthread_t* pth = *i;
|
|---|
| 174 | cout << "joining thread " << pth << endl;
|
|---|
| 175 | pthread_join(*pth, NULL);
|
|---|
| 176 | cout << "thread joined " << pth << endl;
|
|---|
| 177 | }
|
|---|
| 178 | }
|
|---|
| 179 |
|
|---|
| 180 |
|
|---|
| 181 | // -----------ajout cgt vf 19/08/2002
|
|---|
| 182 |
|
|---|
| 183 |
|
|---|
| 184 | void TOIManager::selectTOISegmented(int bufsz, int maxseg)
|
|---|
| 185 | {
|
|---|
| 186 | fgSegmented = true;
|
|---|
| 187 | segBuffsz = bufsz;
|
|---|
| 188 | segMaxseg = maxseg;
|
|---|
| 189 | }
|
|---|
| 190 |
|
|---|
| 191 | void TOIManager::selectTOISeqBuffered(int wsz)
|
|---|
| 192 | {
|
|---|
| 193 | fgSegmented = false;
|
|---|
| 194 | segBuffsz = wsz;
|
|---|
| 195 | }
|
|---|
| 196 |
|
|---|
| 197 | RingPipe& TOIManager::connect(RingProcessor& pout, string out,
|
|---|
| 198 | RingProcessor& pin, string in) {
|
|---|
| 199 | RingPipe* pipe;
|
|---|
| 200 | char buff[128];
|
|---|
| 201 | sprintf(buff, "Ring_[%s-%s]", in.c_str(), out.c_str());
|
|---|
| 202 | if ((pipe = pout.getOutRing(out)) == NULL) {
|
|---|
| 203 | pipe = new RingPipe();
|
|---|
| 204 | pout.addRingOutput(out, pipe);
|
|---|
| 205 | }
|
|---|
| 206 |
|
|---|
| 207 | pin.addRingInput(in, pipe);
|
|---|
| 208 | return(*pipe);
|
|---|
| 209 | }
|
|---|
| 210 |
|
|---|
| 211 | // methode connect de cgt simplifiee et corrigee
|
|---|
| 212 | TOI& TOIManager::connect(TOIProcessor& prout, string out,
|
|---|
| 213 | TOIProcessor& prin, string in, string nom, int wbsz, bool withFlag)
|
|---|
| 214 | {
|
|---|
| 215 | TOI* toi;
|
|---|
| 216 | if (nom.length() < 1) {
|
|---|
| 217 | char buff[128];
|
|---|
| 218 | sprintf(buff, "TOI%s_[%s-%s]", nom.c_str(), in.c_str(), out.c_str());
|
|---|
| 219 | nom = buff;
|
|---|
| 220 | }
|
|---|
| 221 | if (wbsz < 16) wbsz = segBuffsz;
|
|---|
| 222 |
|
|---|
| 223 | // ajout test pour eviter de creer 2 tois en sortie
|
|---|
| 224 | if ((toi=prout.getOutToi(out)) == NULL) {
|
|---|
| 225 | //cout << "toi cree" << endl;
|
|---|
| 226 | if (fgSegmented) toi = new TOISegmented(nom, wbsz, segMaxseg);
|
|---|
| 227 | else toi = new TOISeqBuffered(nom, wbsz);
|
|---|
| 228 | // on ajoute le toi de sortie
|
|---|
| 229 | prout.addOutput(out, toi);
|
|---|
| 230 | } else {
|
|---|
| 231 | //cout << "toi deja cree stop" << endl;
|
|---|
| 232 | }
|
|---|
| 233 |
|
|---|
| 234 | if (withFlag) { // Si c'est un FITSTOIWriter
|
|---|
| 235 | //#ifndef NOCFITSIO
|
|---|
| 236 | FITSTOIWriter* ftw = dynamic_cast< FITSTOIWriter* >(&prin);
|
|---|
| 237 | if (ftw) ftw->addInput(in, toi, withFlag);
|
|---|
| 238 | else prin.addInput(in, toi);
|
|---|
| 239 | //#else
|
|---|
| 240 | // prin.addInput(in, toi);
|
|---|
| 241 | //#endif
|
|---|
| 242 | }
|
|---|
| 243 | else prin.addInput(in, toi);
|
|---|
| 244 | return(*toi);
|
|---|
| 245 | }
|
|---|
| 246 |
|
|---|
| 247 |
|
|---|
| 248 | TOI& TOIManager::connect(TOIProcessor& prout, const char* out,
|
|---|
| 249 | TOIProcessor& prin, const char* in, string nom, int wbsz, bool withFlag)
|
|---|
| 250 | {
|
|---|
| 251 | string outs = out;
|
|---|
| 252 | string ins = in;
|
|---|
| 253 | return connect(prout, outs, prin, ins, nom, wbsz, withFlag);
|
|---|
| 254 | }
|
|---|
| 255 |
|
|---|
| 256 | // ----------- fin ajout cgt
|
|---|
| 257 |
|
|---|
| 258 |
|
|---|
| 259 |
|
|---|
| 260 | // -----------------------------------------------------------------
|
|---|
| 261 | // Classe pour affichage de l'avancement des TOIProcessors
|
|---|
| 262 | // Reza 08/2001
|
|---|
| 263 | // -----------------------------------------------------------------
|
|---|
| 264 |
|
|---|
| 265 | RzProcSampleCounter::RzProcSampleCounter()
|
|---|
| 266 | {
|
|---|
| 267 | _msg = "SampleCounter/Info";
|
|---|
| 268 | _rate = 50;
|
|---|
| 269 | }
|
|---|
| 270 |
|
|---|
| 271 | RzProcSampleCounter::~RzProcSampleCounter()
|
|---|
| 272 | {
|
|---|
| 273 | }
|
|---|
| 274 |
|
|---|
| 275 | long RzProcSampleCounter::PrintStats()
|
|---|
| 276 | {
|
|---|
| 277 | int istart = 0;
|
|---|
| 278 | int iend = 0;
|
|---|
| 279 | long dns_print = 1000;
|
|---|
| 280 | int dns_print_fac = _rate;
|
|---|
| 281 | int nbmax_dns_print = 2;
|
|---|
| 282 |
|
|---|
| 283 | // TOIManager* mgr = TOIManager::getManager();
|
|---|
| 284 |
|
|---|
| 285 | // istart = mgr->getRequestedBegin();
|
|---|
| 286 | // iend = mgr->getRequestedEnd();
|
|---|
| 287 | istart = SampleBegin();
|
|---|
| 288 | iend = SampleEnd();
|
|---|
| 289 |
|
|---|
| 290 | dns_print = (iend-istart)/dns_print_fac;
|
|---|
| 291 | if (dns_print < 1000) dns_print = ((iend-istart) < 1000) ? (iend-istart) : 1000;
|
|---|
| 292 | if (dns_print < 1) dns_print = 1;
|
|---|
| 293 | nbmax_dns_print = (iend-istart)/dns_print;
|
|---|
| 294 |
|
|---|
| 295 | cout << "RzProcSampleCounter::PrintStats() InfoMessage=" << _msg
|
|---|
| 296 | << "\n ... " << _msg << " istart="
|
|---|
| 297 | << istart << " iend= " << iend << " dns_print= " << dns_print
|
|---|
| 298 | << " nbmax_dns_print= " << nbmax_dns_print << endl;
|
|---|
| 299 | // ------------------- Impression continu de stat ------------------------
|
|---|
| 300 | long nb_dns_print = 0;
|
|---|
| 301 | int nb_sleep = 0;
|
|---|
| 302 | long last_sample_count = 0;
|
|---|
| 303 | long processed_samples = 0;
|
|---|
| 304 | long total_sample_count = dns_print*nbmax_dns_print;
|
|---|
| 305 | bool alldone = false;
|
|---|
| 306 | double fracperc = 0.;
|
|---|
| 307 | int fperc = 0;
|
|---|
| 308 | while (!alldone) {
|
|---|
| 309 | processed_samples = ProcessedSampleCount();
|
|---|
| 310 | if ( (processed_samples-last_sample_count > dns_print) ||
|
|---|
| 311 | (processed_samples > total_sample_count-10) ) {
|
|---|
| 312 | last_sample_count = processed_samples;
|
|---|
| 313 | if (nb_dns_print == 0) cout << "\n";
|
|---|
| 314 | nb_dns_print++;
|
|---|
| 315 | fracperc = (double)processed_samples*100./(double)total_sample_count;
|
|---|
| 316 | fperc = (int)(fracperc*100);
|
|---|
| 317 | cout << ">>> " << _msg << ": ProcessedSampleCount()= " << last_sample_count
|
|---|
| 318 | << " Frac done = " << (double)fperc/100. << " %" << endl;
|
|---|
| 319 | if (last_sample_count > total_sample_count-10) alldone = true;
|
|---|
| 320 | nb_sleep = 0;
|
|---|
| 321 | }
|
|---|
| 322 | else if ((nb_sleep+1)%5 == 0) {
|
|---|
| 323 | fracperc = (double)processed_samples*100./(double)total_sample_count;
|
|---|
| 324 | fperc = (int)(fracperc*100);
|
|---|
| 325 | cout << "> " << _msg << ": ProcSamples()= " << processed_samples
|
|---|
| 326 | << " Done = " << " %" << (double)fperc/100.
|
|---|
| 327 | << " NbSleep(1) = " << nb_sleep << endl;
|
|---|
| 328 | }
|
|---|
| 329 | sleep(1); nb_sleep++;
|
|---|
| 330 | }
|
|---|
| 331 |
|
|---|
| 332 | // -----------------------------------------------------------------------
|
|---|
| 333 |
|
|---|
| 334 | return last_sample_count;
|
|---|
| 335 |
|
|---|
| 336 | }
|
|---|
| 337 |
|
|---|
| 338 |
|
|---|
| 339 |
|
|---|
| 340 |
|
|---|
| 341 |
|
|---|
| 342 |
|
|---|
| 343 |
|
|---|
| 344 |
|
|---|
| 345 |
|
|---|
| 346 |
|
|---|
| 347 |
|
|---|
| 348 |
|
|---|
| 349 |
|
|---|
| 350 |
|
|---|