source: XIOS/trunk/src/context_server.cpp @ 998

Last change on this file since 998 was 998, checked in by ymipsl, 5 years ago
  • Clarify the situation about filter chaining order between spatial and arithmetics operation :
    • Arithmetic expression is computed first on the source grid
    • Spatial transform is computed after, transforming the result of arithmetics expression on source grid to the target grid.
    • Time filter will be involved on the target grid after spatial transform
  • @this is working now.
  • Add a field attribute "expr" that is equivalent to the XML content describing the arithmetic expression. This attribute can be inherited.

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: 7.9 KB
Line 
1#include "context_server.hpp"
2#include "buffer_in.hpp"
3#include "type.hpp"
4#include "context.hpp"
5#include "object_template.hpp"
6#include "group_template.hpp"
7#include "attribute_template.hpp"
8#include "domain.hpp"
9#include "field.hpp"
10#include "file.hpp"
11#include "grid.hpp"
12#include "mpi.hpp"
13#include "tracer.hpp"
14#include "timer.hpp"
15#include "cxios.hpp"
16#include "event_scheduler.hpp"
17#include "server.hpp"
18#include <boost/functional/hash.hpp>
19
20
21
22namespace xios
23{
24
25  CContextServer::CContextServer(CContext* parent,MPI_Comm intraComm_,MPI_Comm interComm_)
26  {
27    context=parent;
28    intraComm=intraComm_;
29    MPI_Comm_size(intraComm,&intraCommSize);
30    MPI_Comm_rank(intraComm,&intraCommRank);
31    interComm=interComm_;
32    int flag;
33    MPI_Comm_test_inter(interComm,&flag);
34    if (flag) MPI_Comm_remote_size(interComm,&commSize);
35    else  MPI_Comm_size(interComm,&commSize);
36    currentTimeLine=0;
37    scheduled=false;
38    finished=false;
39
40    boost::hash<string> hashString;
41    hashId=hashString(context->getId());
42
43  }
44  void CContextServer::setPendingEvent(void)
45  {
46    pendingEvent=true;
47  }
48
49  bool CContextServer::hasPendingEvent(void)
50  {
51    return pendingEvent;
52  }
53
54  bool CContextServer::hasFinished(void)
55  {
56    return finished;
57  }
58
59  bool CContextServer::eventLoop(void)
60  {
61    listen();
62    checkPendingRequest();
63    processEvents();
64    return finished;
65  }
66
67  void CContextServer::listen(void)
68  {
69    int rank;
70    int flag;
71    int count;
72    char * addr;
73    MPI_Status status;
74    map<int,CServerBuffer*>::iterator it;
75
76    for(rank=0;rank<commSize;rank++)
77    {
78      if (pendingRequest.find(rank)==pendingRequest.end())
79      {
80        traceOff();
81        MPI_Iprobe(rank,20,interComm,&flag,&status);
82        traceOn();
83        if (flag==true)
84        {
85          it=buffers.find(rank);
86          if (it==buffers.end()) // Receive the buffer size and allocate the buffer
87          {
88            StdSize buffSize = 0;
89            MPI_Recv(&buffSize, 1, MPI_LONG, rank, 20, interComm, &status);
90            mapBufferSize_.insert(std::make_pair(rank, buffSize));
91            it=(buffers.insert(pair<int,CServerBuffer*>(rank,new CServerBuffer(buffSize)))).first;
92          }
93          else
94          {
95            MPI_Get_count(&status,MPI_CHAR,&count);
96            if (it->second->isBufferFree(count))
97            {
98              addr=(char*)it->second->getBuffer(count);
99              MPI_Irecv(addr,count,MPI_CHAR,rank,20,interComm,&pendingRequest[rank]);
100              bufferRequest[rank]=addr;
101            }
102          }
103        }
104      }
105    }
106  }
107
108  void CContextServer::checkPendingRequest(void)
109  {
110    map<int,MPI_Request>::iterator it;
111    list<int> recvRequest;
112    list<int>::iterator itRecv;
113    int rank;
114    int flag;
115    int count;
116    MPI_Status status;
117
118    for(it=pendingRequest.begin();it!=pendingRequest.end();it++)
119    {
120      rank=it->first;
121      traceOff();
122      MPI_Test(& it->second, &flag, &status);
123      traceOn();
124      if (flag==true)
125      {
126        recvRequest.push_back(rank);
127        MPI_Get_count(&status,MPI_CHAR,&count);
128        processRequest(rank,bufferRequest[rank],count);
129      }
130    }
131
132    for(itRecv=recvRequest.begin();itRecv!=recvRequest.end();itRecv++)
133    {
134      pendingRequest.erase(*itRecv);
135      bufferRequest.erase(*itRecv);
136    }
137  }
138
139  void CContextServer::processRequest(int rank, char* buff,int count)
140  {
141
142    CBufferIn buffer(buff,count);
143    char* startBuffer,endBuffer;
144    int size, offset;
145    size_t timeLine;
146    map<size_t,CEventServer*>::iterator it;
147
148    while(count>0)
149    {
150      char* startBuffer=(char*)buffer.ptr();
151      CBufferIn newBuffer(startBuffer,buffer.remain());
152      newBuffer>>size>>timeLine;
153
154      it=events.find(timeLine);
155      if (it==events.end()) it=events.insert(pair<int,CEventServer*>(timeLine,new CEventServer)).first;
156      it->second->push(rank,buffers[rank],startBuffer,size);
157
158      buffer.advance(size);
159      count=buffer.remain();
160    }
161
162  }
163
164  void CContextServer::processEvents(void)
165  {
166    map<size_t,CEventServer*>::iterator it;
167    CEventServer* event;
168
169    it=events.find(currentTimeLine);
170    if (it!=events.end())
171    {
172      event=it->second;
173
174      if (event->isFull())
175      {
176        if (!scheduled && CServer::eventScheduler) // Skip event scheduling for attached mode and reception on client side
177        {
178          CServer::eventScheduler->registerEvent(currentTimeLine,hashId);
179          scheduled=true;
180        }
181        else if (!CServer::eventScheduler || CServer::eventScheduler->queryEvent(currentTimeLine,hashId) )
182        {
183         // When using attached mode, synchronise the processes to avoid that differents event be scheduled by differents processes
184         // The best way to properly solve this problem will be to use the event scheduler also in attached mode
185         // for now just set up a MPI barrier
186         if (!CServer::eventScheduler && CXios::isServer) MPI_Barrier(intraComm) ;
187
188         CTimer::get("Process events").resume();
189         dispatchEvent(*event);
190         CTimer::get("Process events").suspend();
191         pendingEvent=false;
192         delete event;
193         events.erase(it);
194         currentTimeLine++;
195         scheduled = false;
196        }
197      }
198    }
199  }
200
201  CContextServer::~CContextServer()
202  {
203    map<int,CServerBuffer*>::iterator it;
204    for(it=buffers.begin();it!=buffers.end();++it) delete it->second;
205  }
206
207
208  void CContextServer::dispatchEvent(CEventServer& event)
209  {
210    string contextName;
211    string buff;
212    int MsgSize;
213    int rank;
214    list<CEventServer::SSubEvent>::iterator it;
215    CContext::setCurrent(context->getId());
216
217    if (event.classId==CContext::GetType() && event.type==CContext::EVENT_ID_CONTEXT_FINALIZE)
218    {
219      finished=true;
220      info(20)<<"Server Side context <"<<context->getId()<<"> finalized"<<endl;
221      std::map<int, StdSize>::const_iterator itbMap = mapBufferSize_.begin(),
222                                             iteMap = mapBufferSize_.end(), itMap;
223      StdSize totalBuf = 0;
224      for (itMap = itbMap; itMap != iteMap; ++itMap)
225      {
226        report(10)<< " Memory report : Context <"<<context->getId()<<"> : server side : memory used for buffer of each connection to client" << endl
227                  << "  +) With client of rank " << itMap->first << " : " << itMap->second << " bytes " << endl;
228        totalBuf += itMap->second;
229      }
230      context->finalize();
231      report(0)<< " Memory report : Context <"<<context->getId()<<"> : server side : total memory used for buffer "<<totalBuf<<" bytes"<<endl;
232    }
233    else if (event.classId==CContext::GetType()) CContext::dispatchEvent(event);
234    else if (event.classId==CContextGroup::GetType()) CContextGroup::dispatchEvent(event);
235    else if (event.classId==CCalendarWrapper::GetType()) CCalendarWrapper::dispatchEvent(event);
236    else if (event.classId==CDomain::GetType()) CDomain::dispatchEvent(event);
237    else if (event.classId==CDomainGroup::GetType()) CDomainGroup::dispatchEvent(event);
238    else if (event.classId==CAxis::GetType()) CAxis::dispatchEvent(event);
239    else if (event.classId==CAxisGroup::GetType()) CAxisGroup::dispatchEvent(event);
240    else if (event.classId==CScalar::GetType()) CScalar::dispatchEvent(event);
241    else if (event.classId==CScalarGroup::GetType()) CScalarGroup::dispatchEvent(event);
242    else if (event.classId==CGrid::GetType()) CGrid::dispatchEvent(event);
243    else if (event.classId==CGridGroup::GetType()) CGridGroup::dispatchEvent(event);
244    else if (event.classId==CField::GetType()) CField::dispatchEvent(event);
245    else if (event.classId==CFieldGroup::GetType()) CFieldGroup::dispatchEvent(event);
246    else if (event.classId==CFile::GetType()) CFile::dispatchEvent(event);
247    else if (event.classId==CFileGroup::GetType()) CFileGroup::dispatchEvent(event);
248    else if (event.classId==CVariable::GetType()) CVariable::dispatchEvent(event);
249    else
250    {
251      ERROR("void CContextServer::dispatchEvent(CEventServer& event)",<<" Bad event class Id"<<endl);
252    }
253  }
254}
Note: See TracBrowser for help on using the repository browser.