source: XIOS/trunk/src/server.cpp @ 1587

Last change on this file since 1587 was 1587, checked in by ymipsl, 3 years ago

XIOS - OASIS interaction : due to many problem occurring in the oasis/XIOS initialization phase due to a bad order of intialization call from both, you have now the possibily to explicitly inform xios that the servers must call oasis_enddef().
New rules : On model side, before calling oasis_enddef, you must add a call to "xios_oasis_enddef()"
Old rules : oasis_enddef must be call before any call to "xios_context_initialize" otherwise it may lead to a deadlock.
You can use the old rules if the variable <call_oasis_enddef> is set to false (default value is true), and by this way no need to modify the source code of the models

YM

  • Property copyright set to
    Software name : XIOS (Xml I/O Server)
    http://forge.ipsl.jussieu.fr/ioserver
    Creation date : January 2009
    Licence : CeCCIL version2
    see license file in root directory : Licence_CeCILL_V2-en.txt
    or http://www.cecill.info/licences/Licence_CeCILL_V2-en.html
    Holder : CEA/LSCE (Laboratoire des Sciences du CLimat et de l'Environnement)
    CNRS/IPSL (Institut Pierre Simon Laplace)
    Project Manager : Yann Meurdesoif
    yann.meurdesoif@cea.fr
  • Property svn:eol-style set to native
File size: 32.4 KB
Line 
1#include "globalScopeData.hpp"
2#include "xios_spl.hpp"
3#include "cxios.hpp"
4#include "server.hpp"
5#include "client.hpp"
6#include "type.hpp"
7#include "context.hpp"
8#include "object_template.hpp"
9#include "oasis_cinterface.hpp"
10#include <boost/functional/hash.hpp>
11#include <boost/algorithm/string.hpp>
12#include "mpi.hpp"
13#include "tracer.hpp"
14#include "timer.hpp"
15#include "event_scheduler.hpp"
16#include "string_tools.hpp"
17
18namespace xios
19{
20    MPI_Comm CServer::intraComm ;
21    std::list<MPI_Comm> CServer::interCommLeft ;
22    std::list<MPI_Comm> CServer::interCommRight ;
23    std::list<MPI_Comm> CServer::contextInterComms;
24    std::list<MPI_Comm> CServer::contextIntraComms;
25    int CServer::serverLevel = 0 ;
26    int CServer::nbContexts = 0;
27    bool CServer::isRoot = false ;
28    int CServer::rank_ = INVALID_RANK;
29    StdOFStream CServer::m_infoStream;
30    StdOFStream CServer::m_errorStream;
31    map<string,CContext*> CServer::contextList ;
32    vector<int> CServer::sndServerGlobalRanks;
33    bool CServer::finished=false ;
34    bool CServer::is_MPI_Initialized ;
35    CEventScheduler* CServer::eventScheduler = 0;
36
37//---------------------------------------------------------------
38/*!
39 * \fn void CServer::initialize(void)
40 * Creates intraComm for each possible type of servers (classical, primary or secondary).
41 * Creates interComm and stores them into the following lists:
42 *   classical server -- interCommLeft
43 *   primary server -- interCommLeft and interCommRight
44 *   secondary server -- interCommLeft for each pool.
45 *   IMPORTANT: CXios::usingServer2 should NOT be used beyond this function. Use CServer::serverLevel instead.
46 */
47    void CServer::initialize(void)
48    {
49      int initialized ;
50      MPI_Initialized(&initialized) ;
51      if (initialized) is_MPI_Initialized=true ;
52      else is_MPI_Initialized=false ;
53      int rank ;
54
55      // Not using OASIS
56      if (!CXios::usingOasis)
57      {
58
59        if (!is_MPI_Initialized)
60        {
61          MPI_Init(NULL, NULL);
62        }
63        CTimer::get("XIOS").resume() ;
64
65        boost::hash<string> hashString ;
66        unsigned long hashServer = hashString(CXios::xiosCodeId);
67
68        unsigned long* hashAll ;
69        unsigned long* srvLevelAll ;
70
71        int size ;
72        int myColor ;
73        int i,c ;
74        MPI_Comm newComm;
75
76        MPI_Comm_size(CXios::globalComm, &size) ;
77        MPI_Comm_rank(CXios::globalComm, &rank_);
78
79        hashAll=new unsigned long[size] ;
80        MPI_Allgather(&hashServer, 1, MPI_LONG, hashAll, 1, MPI_LONG, CXios::globalComm) ;
81
82        map<unsigned long, int> colors ;
83        map<unsigned long, int> leaders ;
84        map<unsigned long, int>::iterator it ;
85
86        // (1) Establish client leaders, distribute processes between two server levels
87        std::vector<int> srvRanks;
88        for(i=0,c=0;i<size;i++)
89        {
90          if (colors.find(hashAll[i])==colors.end())
91          {
92            colors[hashAll[i]]=c ;
93            leaders[hashAll[i]]=i ;
94            c++ ;
95          }
96          if (CXios::usingServer2)
97            if (hashAll[i] == hashServer)
98              srvRanks.push_back(i);
99        }
100
101        if (CXios::usingServer2)
102        {
103          int reqNbProc = srvRanks.size()*CXios::ratioServer2/100.;
104          if (reqNbProc<1 || reqNbProc==srvRanks.size())
105          {
106            error(0)<<"WARNING: void CServer::initialize(void)"<<endl
107                << "It is impossible to dedicate the requested number of processes = "<<reqNbProc
108                <<" to secondary server. XIOS will run in the classical server mode."<<endl;
109          }
110          else
111          {
112            if (CXios::nbPoolsServer2 == 0) CXios::nbPoolsServer2 = reqNbProc;
113            int firstSndSrvRank = srvRanks.size()*(100.-CXios::ratioServer2)/100. ;
114            int poolLeader = firstSndSrvRank;
115//*********** (1) Comment out the line below to set one process per pool
116            sndServerGlobalRanks.push_back(srvRanks[poolLeader]);
117            int nbPools = CXios::nbPoolsServer2;
118            if ( nbPools > reqNbProc || nbPools < 1)
119            {
120              error(0)<<"WARNING: void CServer::initialize(void)"<<endl
121                  << "It is impossible to allocate the requested number of pools = "<<nbPools
122                  <<" on the secondary server. It will be set so that there is one process per pool."<<endl;
123              nbPools = reqNbProc;
124            }
125            int remainder = ((int) (srvRanks.size()*CXios::ratioServer2/100.)) % nbPools;
126            int procsPerPool = ((int) (srvRanks.size()*CXios::ratioServer2/100.)) / nbPools;
127            for (i=0; i<srvRanks.size(); i++)
128            {
129              if (i >= firstSndSrvRank)
130              {
131                if (rank_ == srvRanks[i])
132                {
133                  serverLevel=2;
134                }
135                poolLeader += procsPerPool;
136                if (remainder != 0)
137                {
138                  ++poolLeader;
139                  --remainder;
140                }
141//*********** (2) Comment out the two lines below to set one process per pool
142                if (poolLeader < srvRanks.size())
143                  sndServerGlobalRanks.push_back(srvRanks[poolLeader]);
144//*********** (3) Uncomment the line below to set one process per pool
145//                sndServerGlobalRanks.push_back(srvRanks[i]);
146              }
147              else
148              {
149                if (rank_ == srvRanks[i]) serverLevel=1;
150              }
151            }
152            if (serverLevel==2)
153            {
154              info(50)<<"The number of secondary server pools is "<< sndServerGlobalRanks.size() <<endl ;
155              for (i=0; i<sndServerGlobalRanks.size(); i++)
156              {
157                if (rank_>= sndServerGlobalRanks[i])
158                {
159                  if ( i == sndServerGlobalRanks.size()-1)
160                  {
161                    myColor = colors.size() + sndServerGlobalRanks[i];
162                  }
163                  else if (rank_< sndServerGlobalRanks[i+1])
164                  {
165                    myColor = colors.size() + sndServerGlobalRanks[i];
166                    break;
167                  }
168                }
169              }
170            }
171          }
172        }
173
174        // (2) Create intraComm
175        if (serverLevel != 2) myColor=colors[hashServer];
176        MPI_Comm_split(CXios::globalComm, myColor, rank_, &intraComm) ;
177
178        // (3) Create interComm
179        if (serverLevel == 0)
180        {
181          int clientLeader;
182          for(it=leaders.begin();it!=leaders.end();it++)
183          {
184            if (it->first!=hashServer)
185            {
186              clientLeader=it->second ;
187              int intraCommSize, intraCommRank ;
188              MPI_Comm_size(intraComm,&intraCommSize) ;
189              MPI_Comm_rank(intraComm,&intraCommRank) ;
190              info(50)<<"intercommCreate::server (classical mode) "<<rank_<<" intraCommSize : "<<intraCommSize
191                       <<" intraCommRank :"<<intraCommRank<<"  clientLeader "<< clientLeader<<endl ;
192
193              MPI_Intercomm_create(intraComm, 0, CXios::globalComm, clientLeader, 0, &newComm) ;
194              interCommLeft.push_back(newComm) ;
195            }
196          }
197        }
198        else if (serverLevel == 1)
199        {
200          int clientLeader, srvSndLeader;
201          int srvPrmLeader ;
202
203          for (it=leaders.begin();it!=leaders.end();it++)
204          {
205            if (it->first != hashServer)
206            {
207              clientLeader=it->second ;
208              int intraCommSize, intraCommRank ;
209              MPI_Comm_size(intraComm, &intraCommSize) ;
210              MPI_Comm_rank(intraComm, &intraCommRank) ;
211              info(50)<<"intercommCreate::server (server level 1) "<<rank_<<" intraCommSize : "<<intraCommSize
212                       <<" intraCommRank :"<<intraCommRank<<"  clientLeader "<< clientLeader<<endl ;
213              MPI_Intercomm_create(intraComm, 0, CXios::globalComm, clientLeader, 0, &newComm) ;
214              interCommLeft.push_back(newComm) ;
215            }
216          }
217
218          for (int i = 0; i < sndServerGlobalRanks.size(); ++i)
219          {
220            int intraCommSize, intraCommRank ;
221            MPI_Comm_size(intraComm, &intraCommSize) ;
222            MPI_Comm_rank(intraComm, &intraCommRank) ;
223            info(50)<<"intercommCreate::client (server level 1) "<<rank_<<" intraCommSize : "<<intraCommSize
224                <<" intraCommRank :"<<intraCommRank<<"  clientLeader "<< sndServerGlobalRanks[i]<<endl ;
225            MPI_Intercomm_create(intraComm, 0, CXios::globalComm, sndServerGlobalRanks[i], 1, &newComm) ;
226            interCommRight.push_back(newComm) ;
227          }
228        }
229        else
230        {
231          int clientLeader;
232          clientLeader = leaders[hashString(CXios::xiosCodeId)];
233          int intraCommSize, intraCommRank ;
234          MPI_Comm_size(intraComm, &intraCommSize) ;
235          MPI_Comm_rank(intraComm, &intraCommRank) ;
236          info(50)<<"intercommCreate::server (server level 2) "<<rank_<<" intraCommSize : "<<intraCommSize
237                   <<" intraCommRank :"<<intraCommRank<<"  clientLeader "<< clientLeader<<endl ;
238
239          MPI_Intercomm_create(intraComm, 0, CXios::globalComm, clientLeader, 1, &newComm) ;
240          interCommLeft.push_back(newComm) ;
241        }
242
243        delete [] hashAll ;
244
245      }
246      // using OASIS
247      else
248      {
249        int size;
250        int myColor;
251        int* srvGlobalRanks;
252        if (!is_MPI_Initialized) oasis_init(CXios::xiosCodeId);
253
254        CTimer::get("XIOS").resume() ;
255        MPI_Comm localComm;
256        oasis_get_localcomm(localComm);
257        MPI_Comm_rank(localComm,&rank_) ;
258
259//      (1) Create server intraComm
260        if (!CXios::usingServer2)
261        {
262          MPI_Comm_dup(localComm, &intraComm);
263        }
264        else
265        {
266          int globalRank;
267          MPI_Comm_size(localComm,&size) ;
268          MPI_Comm_rank(CXios::globalComm,&globalRank) ;
269          srvGlobalRanks = new int[size] ;
270          MPI_Allgather(&globalRank, 1, MPI_INT, srvGlobalRanks, 1, MPI_INT, localComm) ;
271
272          int reqNbProc = size*CXios::ratioServer2/100.;
273          if (reqNbProc < 1 || reqNbProc == size)
274          {
275            error(0)<<"WARNING: void CServer::initialize(void)"<<endl
276                << "It is impossible to dedicate the requested number of processes = "<<reqNbProc
277                <<" to secondary server. XIOS will run in the classical server mode."<<endl;
278            MPI_Comm_dup(localComm, &intraComm);
279          }
280          else
281          {
282            int firstSndSrvRank = size*(100.-CXios::ratioServer2)/100. ;
283            int poolLeader = firstSndSrvRank;
284//*********** (1) Comment out the line below to set one process per pool
285//            sndServerGlobalRanks.push_back(srvGlobalRanks[poolLeader]);
286            int nbPools = CXios::nbPoolsServer2;
287            if ( nbPools > reqNbProc || nbPools < 1)
288            {
289              error(0)<<"WARNING: void CServer::initialize(void)"<<endl
290                  << "It is impossible to allocate the requested number of pools = "<<nbPools
291                  <<" on the secondary server. It will be set so that there is one process per pool."<<endl;
292              nbPools = reqNbProc;
293            }
294            int remainder = ((int) (size*CXios::ratioServer2/100.)) % nbPools;
295            int procsPerPool = ((int) (size*CXios::ratioServer2/100.)) / nbPools;
296            for (int i=0; i<size; i++)
297            {
298              if (i >= firstSndSrvRank)
299              {
300                if (globalRank == srvGlobalRanks[i])
301                {
302                  serverLevel=2;
303                }
304                poolLeader += procsPerPool;
305                if (remainder != 0)
306                {
307                  ++poolLeader;
308                  --remainder;
309                }
310//*********** (2) Comment out the two lines below to set one process per pool
311//                if (poolLeader < size)
312//                  sndServerGlobalRanks.push_back(srvGlobalRanks[poolLeader]);
313//*********** (3) Uncomment the line below to set one process per pool
314                sndServerGlobalRanks.push_back(srvGlobalRanks[i]);
315              }
316              else
317              {
318                if (globalRank == srvGlobalRanks[i]) serverLevel=1;
319              }
320            }
321            if (serverLevel==2)
322            {
323              info(50)<<"The number of secondary server pools is "<< sndServerGlobalRanks.size() <<endl ;
324              for (int i=0; i<sndServerGlobalRanks.size(); i++)
325              {
326                if (globalRank>= sndServerGlobalRanks[i])
327                {
328                  if (i == sndServerGlobalRanks.size()-1)
329                  {
330                    myColor = sndServerGlobalRanks[i];
331                  }
332                  else if (globalRank< sndServerGlobalRanks[i+1])
333                  {
334                    myColor = sndServerGlobalRanks[i];
335                    break;
336                  }
337                }
338              }
339            }
340            if (serverLevel != 2) myColor=0;
341            MPI_Comm_split(localComm, myColor, rank_, &intraComm) ;
342          }
343        }
344
345        string codesId=CXios::getin<string>("oasis_codes_id") ;
346        vector<string> oasisCodeId=splitRegex(codesId,"\\s*,\\s*") ;
347 
348        vector<string>::iterator it ;
349
350        MPI_Comm newComm ;
351        int globalRank ;
352        MPI_Comm_rank(CXios::globalComm,&globalRank);
353
354//      (2) Create interComms with models
355        for(it=oasisCodeId.begin();it!=oasisCodeId.end();it++)
356        {
357          oasis_get_intercomm(newComm,*it) ;
358          if ( serverLevel == 0 || serverLevel == 1)
359          {
360            interCommLeft.push_back(newComm) ;
361            if (rank_==0) MPI_Send(&globalRank,1,MPI_INT,0,0,newComm) ;
362          }
363        }
364
365//      (3) Create interComms between primary and secondary servers
366        int intraCommSize, intraCommRank ;
367        MPI_Comm_size(intraComm,&intraCommSize) ;
368        MPI_Comm_rank(intraComm, &intraCommRank) ;
369
370        if (serverLevel == 1)
371        {
372          for (int i = 0; i < sndServerGlobalRanks.size(); ++i)
373          {
374            int srvSndLeader = sndServerGlobalRanks[i];
375            info(50)<<"intercommCreate::client (server level 1) "<<globalRank<<" intraCommSize : "<<intraCommSize
376                <<" intraCommRank :"<<intraCommRank<<"  clientLeader "<< srvSndLeader<<endl ;
377            MPI_Intercomm_create(intraComm, 0, CXios::globalComm, srvSndLeader, 0, &newComm) ;
378            interCommRight.push_back(newComm) ;
379          }
380        }
381        else if (serverLevel == 2)
382        {
383          info(50)<<"intercommCreate::server (server level 2)"<<globalRank<<" intraCommSize : "<<intraCommSize
384                   <<" intraCommRank :"<<intraCommRank<<"  clientLeader "<< srvGlobalRanks[0] <<endl ;
385          MPI_Intercomm_create(intraComm, 0, CXios::globalComm, srvGlobalRanks[0], 0, &newComm) ;
386          interCommLeft.push_back(newComm) ;
387        }
388        if (CXios::usingServer2) delete [] srvGlobalRanks ;
389
390        bool oasisEnddef=CXios::getin<bool>("call_oasis_enddef",true) ;
391        if (!oasisEnddef) oasis_enddef() ;
392      }
393
394
395      MPI_Comm_rank(intraComm, &rank) ;
396      if (rank==0) isRoot=true;
397      else isRoot=false;
398     
399      eventScheduler = new CEventScheduler(intraComm) ;
400    }
401
402    void CServer::finalize(void)
403    {
404      CTimer::get("XIOS").suspend() ;
405     
406      delete eventScheduler ;
407
408      for (std::list<MPI_Comm>::iterator it = contextInterComms.begin(); it != contextInterComms.end(); it++)
409        MPI_Comm_free(&(*it));
410
411      for (std::list<MPI_Comm>::iterator it = contextIntraComms.begin(); it != contextIntraComms.end(); it++)
412        MPI_Comm_free(&(*it));
413
414//      for (std::list<MPI_Comm>::iterator it = interComm.begin(); it != interComm.end(); it++)
415//        MPI_Comm_free(&(*it));
416
417//        for (std::list<MPI_Comm>::iterator it = interCommLeft.begin(); it != interCommLeft.end(); it++)
418//          MPI_Comm_free(&(*it));
419
420        for (std::list<MPI_Comm>::iterator it = interCommRight.begin(); it != interCommRight.end(); it++)
421          MPI_Comm_free(&(*it));
422
423      MPI_Comm_free(&intraComm);
424
425      if (!is_MPI_Initialized)
426      {
427        if (CXios::usingOasis) oasis_finalize();
428        else MPI_Finalize() ;
429      }
430      report(0)<<"Performance report : Time spent for XIOS : "<<CTimer::get("XIOS server").getCumulatedTime()<<endl  ;
431      report(0)<<"Performance report : Time spent in processing events : "<<CTimer::get("Process events").getCumulatedTime()<<endl  ;
432      report(0)<<"Performance report : Ratio : "<<CTimer::get("Process events").getCumulatedTime()/CTimer::get("XIOS server").getCumulatedTime()*100.<<"%"<<endl  ;
433      report(100)<<CTimer::getAllCumulatedTime()<<endl ;
434    }
435
436     void CServer::eventLoop(void)
437     {
438       bool stop=false ;
439
440       CTimer::get("XIOS server").resume() ;
441       while(!stop)
442       {
443         if (isRoot)
444         {
445           listenContext();
446           listenRootContext();
447           listenOasisEnddef() ;
448           listenRootOasisEnddef() ;
449           if (!finished) listenFinalize() ;
450         }
451         else
452         {
453           listenRootContext();
454           listenRootOasisEnddef() ;
455           if (!finished) listenRootFinalize() ;
456         }
457
458         contextEventLoop() ;
459         if (finished && contextList.empty()) stop=true ;
460         eventScheduler->checkEvent() ;
461       }
462       CTimer::get("XIOS server").suspend() ;
463     }
464
465     void CServer::listenFinalize(void)
466     {
467        list<MPI_Comm>::iterator it, itr;
468        int msg ;
469        int flag ;
470
471        for(it=interCommLeft.begin();it!=interCommLeft.end();it++)
472        {
473           MPI_Status status ;
474           traceOff() ;
475           MPI_Iprobe(0,0,*it,&flag,&status) ;
476           traceOn() ;
477           if (flag==true)
478           {
479              MPI_Recv(&msg,1,MPI_INT,0,0,*it,&status) ;
480              info(20)<<" CServer : Receive client finalize"<<endl ;
481              // Sending server finalize message to secondary servers (if any)
482              for(itr=interCommRight.begin();itr!=interCommRight.end();itr++)
483              {
484                MPI_Send(&msg,1,MPI_INT,0,0,*itr) ;
485              }
486              MPI_Comm_free(&(*it));
487              interCommLeft.erase(it) ;
488              break ;
489            }
490         }
491
492         if (interCommLeft.empty())
493         {
494           int i,size ;
495           MPI_Comm_size(intraComm,&size) ;
496           MPI_Request* requests= new MPI_Request[size-1] ;
497           MPI_Status* status= new MPI_Status[size-1] ;
498
499           for(int i=1;i<size;i++) MPI_Isend(&msg,1,MPI_INT,i,4,intraComm,&requests[i-1]) ;
500           MPI_Waitall(size-1,requests,status) ;
501
502           finished=true ;
503           delete [] requests ;
504           delete [] status ;
505         }
506     }
507
508
509     void CServer::listenRootFinalize()
510     {
511        int flag ;
512        MPI_Status status ;
513        int msg ;
514
515        traceOff() ;
516        MPI_Iprobe(0,4,intraComm, &flag, &status) ;
517        traceOn() ;
518        if (flag==true)
519        {
520           MPI_Recv(&msg,1,MPI_INT,0,4,intraComm,&status) ;
521           finished=true ;
522        }
523      }
524
525
526   /*!
527    * Root process is listening for an order sent by client to call "oasis_enddef".
528    * The root client of a compound send the order (tag 5). It is probed and received.
529    * When the order has been received from each coumpound, the server root process ping the order to the root processes of the secondary levels of servers (if any).
530    * After, it also inform (asynchronous call) other processes of the communicator that the oasis_enddef call must be done
531    */
532   
533     void CServer::listenOasisEnddef(void)
534     {
535        int flag ;
536        MPI_Status status ;
537        list<MPI_Comm>::iterator it;
538        int msg ;
539        static int nbCompound=0 ;
540        int size ;
541        static bool sent=false ;
542        static MPI_Request* allRequests ;
543        static MPI_Status* allStatus ;
544
545
546        if (sent)
547        {
548          MPI_Comm_size(intraComm,&size) ;
549          MPI_Testall(size,allRequests, &flag, allStatus) ;
550          if (flag==true)
551          {
552            delete [] allRequests ;
553            delete [] allStatus ;
554            sent=false ;
555          }
556        }
557       
558
559        for(it=interCommLeft.begin();it!=interCommLeft.end();it++)
560        {
561           MPI_Status status ;
562           traceOff() ;
563           MPI_Iprobe(0,5,*it,&flag,&status) ;  // tags oasis_endded = 5
564           traceOn() ;
565           if (flag==true)
566           {
567              MPI_Recv(&msg,1,MPI_INT,0,5,*it,&status) ; // tags oasis_endded = 5
568              nbCompound++ ;
569              if (nbCompound==interCommLeft.size())
570              {
571                for (std::list<MPI_Comm>::iterator it = interCommRight.begin(); it != interCommRight.end(); it++)
572                {
573                   MPI_Send(&msg,1,MPI_INT,0,5,*it) ; // tags oasis_endded = 5
574                }
575                MPI_Comm_size(intraComm,&size) ;
576                allRequests= new MPI_Request[size] ;
577                allStatus= new MPI_Status[size] ;
578                for(int i=0;i<size;i++) MPI_Isend(&msg,1,MPI_INT,i,5,intraComm,&allRequests[i]) ; // tags oasis_endded = 5
579                sent=true ;
580              }
581           }
582        }
583     }
584     
585   /*!
586    * Processes probes message from root process if oasis_enddef call must be done.
587    * When the order is received it is scheduled to be treated in a synchronized way by all server processes of the communicator
588    */
589     void CServer::listenRootOasisEnddef(void)
590     {
591       int flag ;
592       MPI_Status status ;
593       const int root=0 ;
594       int msg ;
595       static bool eventSent=false ;
596
597       if (eventSent)
598       {
599         boost::hash<string> hashString;
600         size_t hashId = hashString("oasis_enddef");
601         if (eventScheduler->queryEvent(0,hashId))
602         {
603           oasis_enddef() ;
604           eventSent=false ;
605         }
606       }
607         
608       traceOff() ;
609       MPI_Iprobe(root,5,intraComm, &flag, &status) ;
610       traceOn() ;
611       if (flag==true)
612       {
613         MPI_Recv(&msg,1,MPI_INT,root,5,intraComm,&status) ; // tags oasis_endded = 5
614         boost::hash<string> hashString;
615         size_t hashId = hashString("oasis_enddef");
616         eventScheduler->registerEvent(0,hashId);
617         eventSent=true ;
618       }
619     }
620
621
622
623     
624
625     void CServer::listenContext(void)
626     {
627
628       MPI_Status status ;
629       int flag ;
630       static char* buffer ;
631       static MPI_Request request ;
632       static bool recept=false ;
633       int rank ;
634       int count ;
635
636       if (recept==false)
637       {
638         traceOff() ;
639         MPI_Iprobe(MPI_ANY_SOURCE,1,CXios::globalComm, &flag, &status) ;
640         traceOn() ;
641         if (flag==true)
642         {
643           rank=status.MPI_SOURCE ;
644           MPI_Get_count(&status,MPI_CHAR,&count) ;
645           buffer=new char[count] ;
646           MPI_Irecv((void*)buffer,count,MPI_CHAR,rank,1,CXios::globalComm,&request) ;
647           recept=true ;
648         }
649       }
650       else
651       {
652         traceOff() ;
653         MPI_Test(&request,&flag,&status) ;
654         traceOn() ;
655         if (flag==true)
656         {
657           rank=status.MPI_SOURCE ;
658           MPI_Get_count(&status,MPI_CHAR,&count) ;
659           recvContextMessage((void*)buffer,count) ;
660           delete [] buffer ;
661           recept=false ;
662         }
663       }
664     }
665
666     void CServer::recvContextMessage(void* buff,int count)
667     {
668       static map<string,contextMessage> recvContextId;
669       map<string,contextMessage>::iterator it ;
670       CBufferIn buffer(buff,count) ;
671       string id ;
672       int clientLeader ;
673       int nbMessage ;
674
675       buffer>>id>>nbMessage>>clientLeader ;
676
677       it=recvContextId.find(id) ;
678       if (it==recvContextId.end())
679       {
680         contextMessage msg={0,0} ;
681         pair<map<string,contextMessage>::iterator,bool> ret ;
682         ret=recvContextId.insert(pair<string,contextMessage>(id,msg)) ;
683         it=ret.first ;
684       }
685       it->second.nbRecv+=1 ;
686       it->second.leaderRank+=clientLeader ;
687
688       if (it->second.nbRecv==nbMessage)
689       {
690         int size ;
691         MPI_Comm_size(intraComm,&size) ;
692//         MPI_Request* requests= new MPI_Request[size-1] ;
693//         MPI_Status* status= new MPI_Status[size-1] ;
694         MPI_Request* requests= new MPI_Request[size] ;
695         MPI_Status* status= new MPI_Status[size] ;
696
697         CMessage msg ;
698         msg<<id<<it->second.leaderRank;
699         int messageSize=msg.size() ;
700         void * sendBuff = new char[messageSize] ;
701         CBufferOut sendBuffer(sendBuff,messageSize) ;
702         sendBuffer<<msg ;
703
704         // Include root itself in order not to have a divergence
705         for(int i=0; i<size; i++)
706         {
707           MPI_Isend(sendBuff,sendBuffer.count(),MPI_CHAR,i,2,intraComm,&requests[i]) ;
708         }
709
710         recvContextId.erase(it) ;
711         delete [] requests ;
712         delete [] status ;
713
714       }
715     }
716
717     void CServer::listenRootContext(void)
718     {
719       MPI_Status status ;
720       int flag ;
721       static std::vector<void*> buffers;
722       static std::vector<MPI_Request> requests ;
723       static std::vector<int> counts ;
724       static std::vector<bool> isEventRegistered ;
725       static std::vector<bool> isEventQueued ;
726       MPI_Request request;
727
728       int rank ;
729       const int root=0 ;
730       boost::hash<string> hashString;
731       size_t hashId = hashString("RegisterContext");
732
733       // (1) Receive context id from the root, save it into a buffer
734       traceOff() ;
735       MPI_Iprobe(root,2,intraComm, &flag, &status) ;
736       traceOn() ;
737       if (flag==true)
738       {
739         counts.push_back(0);
740         MPI_Get_count(&status,MPI_CHAR,&(counts.back())) ;
741         buffers.push_back(new char[counts.back()]) ;
742         requests.push_back(request);
743         MPI_Irecv((void*)(buffers.back()),counts.back(),MPI_CHAR,root,2,intraComm,&(requests.back())) ;
744         isEventRegistered.push_back(false);
745         isEventQueued.push_back(false);
746         nbContexts++;
747       }
748
749       for (int ctxNb = 0; ctxNb < nbContexts; ctxNb++ )
750       {
751         // (2) If context id is received, register an event
752         MPI_Test(&requests[ctxNb],&flag,&status) ;
753         if (flag==true && !isEventRegistered[ctxNb])
754         {
755           eventScheduler->registerEvent(ctxNb,hashId);
756           isEventRegistered[ctxNb] = true;
757         }
758         // (3) If event has been scheduled, call register context
759         if (eventScheduler->queryEvent(ctxNb,hashId) && !isEventQueued[ctxNb])
760         {
761           registerContext(buffers[ctxNb],counts[ctxNb]) ;
762           isEventQueued[ctxNb] = true;
763           delete [] buffers[ctxNb] ;
764         }
765       }
766
767     }
768
769     void CServer::registerContext(void* buff, int count, int leaderRank)
770     {
771       string contextId;
772       CBufferIn buffer(buff, count);
773//       buffer >> contextId;
774       buffer >> contextId>>leaderRank;
775       CContext* context;
776
777       info(20) << "CServer : Register new Context : " << contextId << endl;
778
779       if (contextList.find(contextId) != contextList.end())
780         ERROR("void CServer::registerContext(void* buff, int count, int leaderRank)",
781               << "Context '" << contextId << "' has already been registred");
782
783       context=CContext::create(contextId);
784       contextList[contextId]=context;
785
786       // Primary or classical server: create communication channel with a client
787       // (1) create interComm (with a client)
788       // (2) initialize client and server (contextClient and contextServer)
789       MPI_Comm inter;
790       if (serverLevel < 2)
791       {
792         MPI_Comm contextInterComm;
793         MPI_Intercomm_create(intraComm, 0, CXios::globalComm, leaderRank, 10+leaderRank, &contextInterComm);
794         MPI_Intercomm_merge(contextInterComm,1,&inter);
795         MPI_Barrier(inter);
796         MPI_Comm_free(&inter);
797         context->initServer(intraComm,contextInterComm);
798         contextInterComms.push_back(contextInterComm);
799
800       }
801       // Secondary server: create communication channel with a primary server
802       // (1) duplicate interComm with a primary server
803       // (2) initialize client and server (contextClient and contextServer)
804       // Remark: in the case of the secondary server there is no need to create an interComm calling MPI_Intercomm_create,
805       //         because interComm of CContext is defined on the same processes as the interComm of CServer.
806       //         So just duplicate it.
807       else if (serverLevel == 2)
808       {
809         MPI_Comm_dup(interCommLeft.front(), &inter);
810         contextInterComms.push_back(inter);
811         context->initServer(intraComm, contextInterComms.back());
812       }
813
814       // Primary server:
815       // (1) send create context message to secondary servers
816       // (2) initialize communication channels with secondary servers (create contextClient and contextServer)
817       if (serverLevel == 1)
818       {
819         int i = 0, size;
820         MPI_Comm_size(intraComm, &size) ;
821         for (std::list<MPI_Comm>::iterator it = interCommRight.begin(); it != interCommRight.end(); it++, ++i)
822         {
823           StdString str = contextId +"_server_" + boost::lexical_cast<string>(i);
824           CMessage msg;
825           int messageSize;
826           msg<<str<<size<<rank_ ;
827           messageSize = msg.size() ;
828           buff = new char[messageSize] ;
829           CBufferOut buffer(buff,messageSize) ;
830           buffer<<msg ;
831           MPI_Send(buff, buffer.count(), MPI_CHAR, sndServerGlobalRanks[i], 1, CXios::globalComm) ;
832           MPI_Comm_dup(*it, &inter);
833           contextInterComms.push_back(inter);
834           MPI_Comm_dup(intraComm, &inter);
835           contextIntraComms.push_back(inter);
836           context->initClient(contextIntraComms.back(), contextInterComms.back()) ;
837           delete [] buff ;
838         }
839       }
840     }
841
842     void CServer::contextEventLoop(bool enableEventsProcessing /*= true*/)
843     {
844       bool isFinalized ;
845       map<string,CContext*>::iterator it ;
846
847       for(it=contextList.begin();it!=contextList.end();it++)
848       {
849         isFinalized=it->second->isFinalized();
850         if (isFinalized)
851         {
852           contextList.erase(it) ;
853           break ;
854         }
855         else
856           it->second->checkBuffersAndListen(enableEventsProcessing);
857       }
858     }
859
860     //! Get rank of the current process in the intraComm
861     int CServer::getRank()
862     {
863       int rank;
864       MPI_Comm_rank(intraComm,&rank);
865       return rank;
866     }
867
868     vector<int>& CServer::getSecondaryServerGlobalRanks()
869     {
870       return sndServerGlobalRanks;
871     }
872
873    /*!
874    * Open a file specified by a suffix and an extension and use it for the given file buffer.
875    * The file name will be suffix+rank+extension.
876    *
877    * \param fileName[in] protype file name
878    * \param ext [in] extension of the file
879    * \param fb [in/out] the file buffer
880    */
881    void CServer::openStream(const StdString& fileName, const StdString& ext, std::filebuf* fb)
882    {
883      StdStringStream fileNameClient;
884      int numDigit = 0;
885      int size = 0;
886      int id;
887      MPI_Comm_size(CXios::globalComm, &size);
888      while (size)
889      {
890        size /= 10;
891        ++numDigit;
892      }
893      id = rank_; //getRank();
894
895      fileNameClient << fileName << "_" << std::setfill('0') << std::setw(numDigit) << id << ext;
896      fb->open(fileNameClient.str().c_str(), std::ios::out);
897      if (!fb->is_open())
898        ERROR("void CServer::openStream(const StdString& fileName, const StdString& ext, std::filebuf* fb)",
899              << std::endl << "Can not open <" << fileNameClient.str() << "> file to write the server log(s).");
900    }
901
902    /*!
903    * \brief Open a file stream to write the info logs
904    * Open a file stream with a specific file name suffix+rank
905    * to write the info logs.
906    * \param fileName [in] protype file name
907    */
908    void CServer::openInfoStream(const StdString& fileName)
909    {
910      std::filebuf* fb = m_infoStream.rdbuf();
911      openStream(fileName, ".out", fb);
912
913      info.write2File(fb);
914      report.write2File(fb);
915    }
916
917    //! Write the info logs to standard output
918    void CServer::openInfoStream()
919    {
920      info.write2StdOut();
921      report.write2StdOut();
922    }
923
924    //! Close the info logs file if it opens
925    void CServer::closeInfoStream()
926    {
927      if (m_infoStream.is_open()) m_infoStream.close();
928    }
929
930    /*!
931    * \brief Open a file stream to write the error log
932    * Open a file stream with a specific file name suffix+rank
933    * to write the error log.
934    * \param fileName [in] protype file name
935    */
936    void CServer::openErrorStream(const StdString& fileName)
937    {
938      std::filebuf* fb = m_errorStream.rdbuf();
939      openStream(fileName, ".err", fb);
940
941      error.write2File(fb);
942    }
943
944    //! Write the error log to standard error output
945    void CServer::openErrorStream()
946    {
947      error.write2StdErr();
948    }
949
950    //! Close the error log file if it opens
951    void CServer::closeErrorStream()
952    {
953      if (m_errorStream.is_open()) m_errorStream.close();
954    }
955}
Note: See TracBrowser for help on using the repository browser.