source: XIOS/trunk/src/node/axis.cpp @ 568

Last change on this file since 568 was 568, checked in by mhnguyen, 7 years ago

Implementing discovering algorithm of server index

+) Implement the algorithm with only one level
+) Remove some redundant functions, corrects some interface

Test
+) On Curie
+) Test passed and results are correct

  • 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:executable set to *
File size: 10.1 KB
Line 
1#include "axis.hpp"
2
3#include "attribute_template.hpp"
4#include "object_template.hpp"
5#include "group_template.hpp"
6#include "message.hpp"
7#include "type.hpp"
8#include "context.hpp"
9#include "context_client.hpp"
10#include "xmlioserver_spl.hpp"
11
12namespace xios {
13
14   /// ////////////////////// Définitions ////////////////////// ///
15
16   CAxis::CAxis(void)
17      : CObjectTemplate<CAxis>()
18      , CAxisAttributes(), isChecked(false), relFiles(), baseRefObject(), areClientAttributesChecked_(false)
19   { /* Ne rien faire de plus */ }
20
21   CAxis::CAxis(const StdString & id)
22      : CObjectTemplate<CAxis>(id)
23      , CAxisAttributes(), isChecked(false), relFiles(), baseRefObject(), areClientAttributesChecked_(false)
24   { /* Ne rien faire de plus */ }
25
26   CAxis::~CAxis(void)
27   { /* Ne rien faire de plus */ }
28
29   ///---------------------------------------------------------------
30
31   const std::set<StdString> & CAxis::getRelFiles(void) const
32   {
33      return (this->relFiles);
34   }
35
36   bool CAxis::IsWritten(const StdString & filename) const
37   {
38      return (this->relFiles.find(filename) != this->relFiles.end());
39   }
40
41   void CAxis::addRelFile(const StdString & filename)
42   {
43      this->relFiles.insert(filename);
44   }
45
46   //----------------------------------------------------------------
47
48   StdString CAxis::GetName(void)   { return (StdString("axis")); }
49   StdString CAxis::GetDefName(void){ return (CAxis::GetName()); }
50   ENodeType CAxis::GetType(void)   { return (eAxis); }
51
52   //----------------------------------------------------------------
53
54   void CAxis::checkAttributes(void)
55   {
56      if (this->size.isEmpty())
57         ERROR("CAxis::checkAttributes(void)",
58               << "Attribute <size> of the axis [ id = '" << getId() << "' , context = '" << CObjectFactory::GetCurrentContextId() << "' ] must be specified");
59      StdSize size = this->size.getValue();
60
61      if (!this->ibegin.isEmpty())
62      {
63        StdSize ibegin = this->ibegin.getValue();
64        if ((ibegin < 0) || (ibegin > size-1))
65          ERROR("CAxis::checkAttributes(void)",
66                << "Attribute <ibegin> of the axis [ id = '" << getId() << "' , context = '" << CObjectFactory::GetCurrentContextId() << "' ] must be non-negative and smaller than size-1");
67      }
68      else this->ibegin.setValue(0);
69
70      if (!this->ni.isEmpty())
71      {
72        StdSize ni = this->ni.getValue();
73        if ((ni < 0) || (ni > size))
74          ERROR("CAxis::checkAttributes(void)",
75                << "Attribute <ni> of the axis [ id = '" << getId() << "' , context = '" << CObjectFactory::GetCurrentContextId() << "' ] must be non-negative and smaller than size");
76      }
77      else this->ni.setValue(size);
78
79              << "One or more attributes among <zoom_begin>, <zoom_end>, <zoom_size> of axis [ id = '" << getId() << "' , context = '" << CObjectFactory::GetCurrentContextId() << "' ] are not well specified");
80      StdSize true_size = value.numElements();
81      if (size != true_size)
82         ERROR("CAxis::checkAttributes(void)",
83               << "The array \'value\' of axis [ id = '" << getId() << "' , context = '" << CObjectFactory::GetCurrentContextId() << "' ] has a different size that the one defined by the \'size\' attribute");
84
85      this->checkData();
86      this->checkMask();
87      this->checkZoom();
88   }
89
90   void CAxis::checkData()
91   {
92      if (data_begin.isEmpty()) data_begin.setValue(0);
93      if (!data_n.isEmpty() && data_n.getValue() <= 0)
94      {
95        ERROR("CAxis::checkData(void)",
96              << "Data dimension is negative (data_n).");
97      }
98      else if (data_n.isEmpty())
99        data_n.setValue(ni.getValue());
100
101      if (data_index.isEmpty())
102      {
103        int dn = data_n.getValue();
104        data_index.resize(dn);
105        for (int i = 0; i < dn; ++i) data_index(i) = (i+1);
106      }
107   }
108
109   void CAxis::checkZoom(void)
110   {
111      StdSize zoom_begin,zoom_end, zoom_size, axisSize;
112
113      zoom_begin = (this->zoom_begin.isEmpty()) ?  0 : this->zoom_begin.getValue() ;
114      zoom_size  = (this->zoom_size.isEmpty()) ?  size.getValue() : this->zoom_size.getValue() ;
115      zoom_end   = (this->zoom_end.isEmpty()) ?  (size.getValue() - 1) : this->zoom_end.getValue() ;
116
117      if (this->zoom_begin.isEmpty()) zoom_begin=zoom_end-zoom_size+1 ;
118      if (this->zoom_end.isEmpty()) zoom_end=zoom_begin+zoom_size-1 ;
119      if (this->zoom_size.isEmpty()) zoom_size=zoom_end-zoom_begin+1 ;
120      axisSize = size.getValue();
121
122      if ( (zoom_begin < 0) || (zoom_begin > axisSize-1) || (zoom_end<0) || (zoom_end>axisSize-1) || (zoom_size<1) || (zoom_size>axisSize) || (zoom_begin>zoom_end))
123        ERROR("CAxis::checkAttributes(void)",<< "One or more attribut of <zoom_begin>, <zoom_end>, <zoom_size>, are not well specified") ;
124
125      this->zoom_begin.setValue(zoom_begin) ;
126      this->zoom_end.setValue(zoom_end) ;
127      this->zoom_size.setValue(zoom_size) ;
128
129      // compute client zoom indices
130//      zoom_begin_client = ibegin_client > zoom_begin ? begin_client : zoom_begin ;
131//      zoom_end_client   = iend_client < zoom_end ? iend_client : zoom_end ;
132//      zoom_size_client  = zoom_end_client-zoom_begin_client+1 ;
133//      if (zoom_ni_client<0) zoom_ni_client=0 ;
134   }
135
136   void CAxis::checkMask()
137   {
138      int begin_mask = 0,
139          end_mask = ni.getValue()-1;
140
141      if (!zoom_begin.isEmpty())
142      {
143         int zoom_end = zoom_begin.getValue() + zoom_size.getValue() - 1;
144
145         begin_mask = std::max(ibegin.getValue(), zoom_begin.getValue());
146         end_mask   = std::min(ibegin.getValue() + ni.getValue()-1, zoom_end);
147
148         begin_mask -= ibegin.getValue();
149         end_mask   -= ibegin.getValue();
150      }
151
152
153      if (!mask.isEmpty())
154      {
155         if (mask.extent(0) != ni)
156            ERROR("CAxis::checkMask(void)",
157                  <<"the mask has not the same size than the local axis"<<endl
158                  <<"Local size is "<<ni<<"x"<<endl
159                  <<"Mask size is "<<mask.extent(0)<<"x");
160         for (int i = 0; i < ni; ++i)
161         {
162           if (i < begin_mask && i > end_mask)  mask(i) = false;
163         }
164      }
165      else // (!mask.hasValue())
166      { // Si aucun masque n'est défini,
167        // on en crée un nouveau qui valide l'intégralité du domaine.
168         mask.resize(ni);
169         for (int i = 0; i < ni.getValue(); ++i)
170         {
171               if (i >= begin_mask && i <= end_mask)
172                 mask(i) = true;
173               else  mask(i) = false;
174         }
175      }
176   }
177
178  bool CAxis::dispatchEvent(CEventServer& event)
179   {
180      if (SuperClass::dispatchEvent(event)) return true ;
181      else
182      {
183        switch(event.type)
184        {
185           case EVENT_ID_SERVER_ATTRIBUT :
186             recvServerAttribut(event) ;
187             return true ;
188             break ;
189           default :
190             ERROR("bool CContext::dispatchEvent(CEventServer& event)",
191                    <<"Unknown Event") ;
192           return false ;
193         }
194      }
195   }
196
197   void CAxis::checkAttributesOnClient(const std::vector<int>& globalDim, int orderPositionInGrid,
198                                       CServerDistributionDescription::ServerDistributionType distType)
199   {
200     if (this->areClientAttributesChecked_) return;
201     this->checkAttributes();
202
203     CContext* context=CContext::getCurrent() ;
204     if (context->hasClient)
205     {
206       computeServerIndex(globalDim, orderPositionInGrid, distType);
207     }
208
209     this->areClientAttributesChecked_ = true;
210   }
211
212   void CAxis::computeServerIndex(const std::vector<int>& globalDim, int orderPositionInGrid,
213                                  CServerDistributionDescription::ServerDistributionType distType)
214   {
215     CServerDistributionDescription serverDescription(globalDim);
216
217     CContext* context=CContext::getCurrent() ;
218     CContextClient* client=context->client ;
219     int nbServer=client->serverSize ;
220     int serverRank=client->getServerLeader() ;
221
222     serverDescription.computeServerDistribution(nbServer, false, distType);
223     std::vector<std::vector<int> > serverIndexBegin = serverDescription.getServerIndexBegin();
224     std::vector<std::vector<int> > serverDimensionSizes = serverDescription.getServerDimensionSizes();
225     begin_srv = (serverIndexBegin[serverRank])[orderPositionInGrid];
226     ni_srv = serverDimensionSizes[serverRank][orderPositionInGrid];
227     end_srv = begin_srv+ni_srv-1;
228   }
229
230   // Send all checked attributes to server
231   void CAxis::sendCheckedAttributes(const std::vector<int>& globalDim, int orderPositionInGrid,
232                                     CServerDistributionDescription::ServerDistributionType distType)
233   {
234     if (!this->areClientAttributesChecked_) checkAttributesOnClient(globalDim,
235                                                                     orderPositionInGrid,
236                                                                     distType);
237     CContext* context=CContext::getCurrent() ;
238
239     if (this->isChecked) return;
240     if (context->hasClient)
241     {
242       sendServerAttribut() ;
243     }
244
245     this->isChecked = true;
246   }
247
248  void CAxis::sendServerAttribut(void)
249  {
250    CContext* context=CContext::getCurrent();
251    CContextClient* client=context->client;
252
253    CEventClient event(getType(),EVENT_ID_SERVER_ATTRIBUT) ;
254    if (client->isServerLeader())
255    {
256      CMessage msg ;
257      msg<<this->getId() ;
258      msg<<ni_srv<<begin_srv<<end_srv;
259      event.push(client->getServerLeader(),1,msg) ;
260      client->sendEvent(event) ;
261    }
262    else client->sendEvent(event) ;
263  }
264
265  void CAxis::recvServerAttribut(CEventServer& event)
266  {
267    CBufferIn* buffer=event.subEvents.begin()->buffer;
268    string axisId ;
269    *buffer>>axisId ;
270    get(axisId)->recvServerAttribut(*buffer) ;
271  }
272
273  void CAxis::recvServerAttribut(CBufferIn& buffer)
274  {
275    int zoom_end = zoom_begin.getValue()+zoom_size.getValue()-1;
276    int ni_srv, begin_srv, end_srv;
277
278    buffer>>ni_srv>>begin_srv>>end_srv;
279
280    zoom_begin_srv = zoom_begin.getValue() > begin_srv ? zoom_begin.getValue() : begin_srv ;
281    zoom_end_srv   = zoom_end < end_srv ? zoom_end : end_srv ;
282    zoom_size_srv  = zoom_end_srv-zoom_begin_srv+1 ;
283
284    if (zoom_size_srv<=0)
285    {
286      zoom_begin_srv=0 ; zoom_end_srv=0 ; zoom_size_srv=0 ;
287    }
288  }
289
290   DEFINE_REF_FUNC(Axis,axis)
291
292   ///---------------------------------------------------------------
293
294} // namespace xios
Note: See TracBrowser for help on using the repository browser.