source: XIOS/trunk/src/node/field.cpp @ 553

Last change on this file since 553 was 551, checked in by mhnguyen, 9 years ago

Redesigning grid structure

+) Add an intermediate class to calculate distribution on client and servers
+) Change all index of attributes to zero (0), instead of one(1)

Test
+) On Curie
+) Test new features passes but some data are still shifted
  • 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
File size: 32.1 KB
Line 
1#include "field.hpp"
2
3#include "attribute_template.hpp"
4#include "object_template.hpp"
5#include "group_template.hpp"
6
7#include "node_type.hpp"
8#include "calendar_util.hpp"
9#include "message.hpp"
10#include "xmlioserver_spl.hpp"
11#include "type.hpp"
12#include "context_client.hpp"
13#include <set>
14
15namespace xios{
16
17   /// ////////////////////// Définitions ////////////////////// ///
18
19   CField::CField(void)
20      : CObjectTemplate<CField>(), CFieldAttributes()
21      , refObject(), baseRefObject()
22      , grid(), file()
23      , freq_operation(), freq_write()
24      , nstep(0)
25      , last_Write(), last_operation()
26      , foperation(), hasInstantData(false), hasExpression(false)
27      , active(false) , hasOutputFile(false),hasFieldOut(false), slotUpdateDate(NULL)
28      , processed(false), domAxisIds_("",""), areAllReferenceSolved(false), areAllExpressionBuilt(false)
29      { setVirtualVariableGroup() ; }
30
31   CField::CField(const StdString & id)
32      : CObjectTemplate<CField>(id), CFieldAttributes()
33      , refObject(), baseRefObject()
34      , grid(), file()
35      , freq_operation(), freq_write()
36      , nstep(0)
37      , last_Write(), last_operation()
38      , foperation(), hasInstantData(false), hasExpression(false)
39      , active(false), hasOutputFile(false), hasFieldOut(false), slotUpdateDate(NULL)
40      , processed(false), domAxisIds_("",""), areAllReferenceSolved(false), areAllExpressionBuilt(false)
41   { setVirtualVariableGroup() ; }
42
43   CField::~CField(void)
44   {
45//      this->grid.reset() ;
46//      this->file.reset() ;
47      this->foperation.reset() ;
48      if (hasExpression) delete expression ;
49      if (slotUpdateDate!=NULL) delete slotUpdateDate ;
50
51   }
52
53
54  //----------------------------------------------------------------
55
56   void CField::setVirtualVariableGroup(CVariableGroup* newVVariableGroup)
57   {
58      this->vVariableGroup = newVVariableGroup;
59   }
60
61   void CField::setVirtualVariableGroup(void)
62   {
63      this->setVirtualVariableGroup(CVariableGroup::create());
64   }
65
66   CVariableGroup* CField::getVirtualVariableGroup(void) const
67   {
68      return (this->vVariableGroup);
69   }
70
71
72   std::vector<CVariable*> CField::getAllVariables(void) const
73   {
74      return (this->vVariableGroup->getAllChildren());
75   }
76
77   void CField::solveDescInheritance(bool apply, const CAttributeMap * const parent)
78   {
79      SuperClassAttribute::setAttributes(parent,apply);
80      this->getVirtualVariableGroup()->solveDescInheritance(apply, NULL);
81   }
82   //----------------------------------------------------------------
83   //----------------------------------------------------------------
84
85   bool CField::updateDataServer
86      (const CDate & currDate,
87       const std::deque< CArray<double, 1>* > storedClient)
88   {
89      const CDate opeDate      = *last_operation + freq_operation;
90      const CDate writeDate    = *last_Write     + freq_write;
91
92      if (opeDate <= currDate)
93      {
94         if (this->data.numElements() != this->grid->storeIndex[0]->numElements())
95         {
96            this->data.resize(this->grid->storeIndex[0] ->numElements());
97         }
98         CArray<double,1> input(data.numElements()) ;
99         this->grid->inputFieldServer(storedClient, input);
100         (*this->foperation)(input);
101         *last_operation = currDate;
102      }
103      if (writeDate < (currDate + freq_operation))
104      {
105         this->foperation->final();
106         this->incrementNStep();
107         *last_Write = writeDate;
108         return (true);
109      }
110      return (false);
111   }
112
113   bool CField::dispatchEvent(CEventServer& event)
114  {
115
116    if (SuperClass::dispatchEvent(event)) return true ;
117    else
118    {
119      switch(event.type)
120      {
121        case EVENT_ID_UPDATE_DATA :
122          recvUpdateData(event) ;
123          return true ;
124          break ;
125
126            case EVENT_ID_ADD_VARIABLE :
127             recvAddVariable(event) ;
128             return true ;
129             break ;
130
131           case EVENT_ID_ADD_VARIABLE_GROUP :
132             recvAddVariableGroup(event) ;
133             return true ;
134             break ;
135
136        default :
137          ERROR("bool CField::dispatchEvent(CEventServer& event)",<<"Unknown Event") ;
138          return false ;
139      }
140    }
141  }
142
143  void CField::sendUpdateData(void)
144  {
145    CContext* context = CContext::getCurrent() ;
146    CContextClient* client=context->client ;
147
148    CEventClient event(getType(),EVENT_ID_UPDATE_DATA) ;
149
150    map<int,CArray<int, 1>* >::iterator it ;
151    list<shared_ptr<CMessage> > list_msg ;
152    list< CArray<double,1>* > list_data ;
153
154     for(it=grid->storeIndex_toSrv.begin();it!=grid->storeIndex_toSrv.end();it++)
155    {
156      int rank=(*it).first ;
157      CArray<int,1>& index = *(it->second) ;
158      CArray<double,1> data_tmp(index.numElements()) ;
159      for(int n=0;n<data_tmp.numElements();n++) data_tmp(n)=data(index(n)) ;
160      list_msg.push_back(shared_ptr<CMessage>(new CMessage)) ;
161      list_data.push_back(new CArray<double,1>(data_tmp)) ;
162      *list_msg.back()<<getId()<<*list_data.back() ;
163      event.push(rank,grid->nbSenders[rank],*list_msg.back()) ;
164    }
165    client->sendEvent(event) ;
166
167    for(list< CArray<double,1>* >::iterator it=list_data.begin();it!=list_data.end();it++) delete *it ;
168  }
169
170  void CField::recvUpdateData(CEventServer& event)
171  {
172    vector<int> ranks ;
173    vector<CBufferIn*> buffers ;
174
175    list<CEventServer::SSubEvent>::iterator it ;
176    string fieldId ;
177
178    for (it=event.subEvents.begin();it!=event.subEvents.end();++it)
179    {
180      int rank=it->rank;
181      CBufferIn* buffer=it->buffer;
182      *buffer>>fieldId ;
183      ranks.push_back(rank) ;
184      buffers.push_back(buffer) ;
185    }
186    get(fieldId)->recvUpdateData(ranks,buffers) ;
187  }
188
189  void  CField::recvUpdateData(vector<int>& ranks, vector<CBufferIn*>& buffers)
190  {
191
192    if (data_srv.empty())
193    {
194//      for(map<int, CArray<int, 1>* >::iterator it=grid->out_i_fromClient.begin();it!=grid->out_i_fromClient.end();it++)
195      for(map<int, CArray<size_t, 1>* >::iterator it=grid->outIndexFromClient.begin();it!=grid->outIndexFromClient.end();++it)
196      {
197        int rank=it->first ;
198        CArray<double,1> data_tmp(it->second->numElements()) ;
199        data_srv.insert( pair<int, CArray<double,1>* >(rank, new CArray<double,1>(data_tmp) ) ) ;
200        foperation_srv.insert(pair<int,boost::shared_ptr<func::CFunctor> >(rank,boost::shared_ptr<func::CFunctor>(new func::CInstant(*data_srv[rank])))) ;
201      }
202    }
203
204    CContext* context = CContext::getCurrent() ;
205    const CDate & currDate = context->getCalendar()->getCurrentDate();
206    const CDate opeDate      = *last_operation_srv + freq_operation_srv;
207    const CDate writeDate    = *last_Write_srv     + freq_write_srv;
208
209
210
211    if (opeDate <= currDate)
212    {
213      for(int n=0;n<ranks.size();n++)
214      {
215        CArray<double,1> data_tmp ;
216        *buffers[n]>>data_tmp ;
217        (*foperation_srv[ranks[n]])(data_tmp) ;
218      }
219      *last_operation_srv = currDate;
220    }
221
222    if (writeDate < (currDate + freq_operation_srv))
223    {
224      for(int n=0;n<ranks.size();n++)
225      {
226        this->foperation_srv[ranks[n]]->final();
227      }
228
229      *last_Write_srv = writeDate;
230      writeField() ;
231      *lastlast_Write_srv=*last_Write_srv;
232    }
233  }
234
235  void CField::writeField(void)
236  {
237    if (! getRelFile()->allDomainEmpty )
238      if (! grid->domain->isEmpty() || getRelFile()->type == CFile::type_attr::one_file)
239      {
240        getRelFile()->checkFile();
241        this->incrementNStep();
242        getRelFile()->getDataOutput()->writeFieldData(CField::get(this));
243      }
244  }
245   //----------------------------------------------------------------
246
247   void CField::setRelFile(CFile* _file)
248   {
249      this->file = _file;
250      hasOutputFile=true ;
251   }
252
253   //----------------------------------------------------------------
254
255   StdString CField::GetName(void)   { return (StdString("field")); }
256   StdString CField::GetDefName(void){ return (CField::GetName()); }
257   ENodeType CField::GetType(void)   { return (eField); }
258
259   //----------------------------------------------------------------
260
261   CGrid* CField::getRelGrid(void) const
262   {
263      return (this->grid);
264   }
265
266   //----------------------------------------------------------------
267
268   CFile* CField::getRelFile(void) const
269   {
270      return (this->file);
271   }
272
273   StdSize CField::getNStep(void) const
274   {
275      return (this->nstep);
276   }
277
278   void CField::incrementNStep(void)
279   {
280      this->nstep++;
281   }
282
283   void CField::resetNStep(void)
284   {
285      this->nstep=0;
286   }
287
288   //----------------------------------------------------------------
289
290   const CDuration & CField::getFreqOperation(void) const
291   {
292      return (this->freq_operation);
293   }
294
295   //----------------------------------------------------------------
296   const CDuration & CField::getFreqWrite(void) const
297   {
298      return (this->freq_write);
299   }
300
301   //----------------------------------------------------------------
302
303   boost::shared_ptr<func::CFunctor> CField::getFieldOperation(void) const
304   {
305      return (this->foperation);
306   }
307
308
309   bool CField::isActive(void) const
310   {
311      return (!this->refObject.empty());
312   }
313   //----------------------------------------------------------------
314
315   CArray<double, 1> CField::getData(void) const
316   {
317      return(this->data);
318   }
319
320   //----------------------------------------------------------------
321
322   boost::shared_ptr<CDate> CField::getLastWriteDate(void) const
323   {
324      return(this->last_Write);
325   }
326
327   //----------------------------------------------------------------
328
329   boost::shared_ptr<CDate> CField::getLastOperationDate(void) const
330   {
331      return(this->last_operation);
332   }
333
334   //----------------------------------------------------------------
335
336//   void CField::processEnabledField(void)
337//   {
338//      if (!processed)
339//      {
340//        processed=true ;
341//        solveRefInheritance(true) ;
342//        solveBaseReference() ;
343//        solveOperation() ;
344//        solveGridReference() ;
345//
346//        if (hasDirectFieldReference()) baseRefObject->processEnabledField() ;
347//        buildExpression();
348//        active=true;
349//      }
350//    }
351
352   void CField::solveAllReferenceEnabledField(bool doSending2Sever)
353   {
354     CContext* context = CContext::getCurrent();
355     if (!areAllReferenceSolved)
356     {
357        areAllReferenceSolved = true;
358        if (context->hasClient)
359        {
360          solveRefInheritance(true);
361          solveBaseReference();
362        }
363
364        solveOperation();
365        solveGridReference();
366     }
367     solveGridDomainAxisRef(doSending2Sever);
368     solveCheckMaskIndex(doSending2Sever);
369   }
370
371   std::map<int, StdSize> CField::getGridDataSize()
372   {
373     return grid->getConnectedServerDataSize();
374   }
375
376   void CField::buildAllExpressionEnabledField()
377   {
378     if (!areAllReferenceSolved) solveAllReferenceEnabledField(true);
379     if (!areAllExpressionBuilt)
380     {
381       areAllExpressionBuilt = true;
382//       solveCheckMaskIndex(true);
383//       solveCheckMaskIndex();
384       if (hasDirectFieldReference()) baseRefObject->buildAllExpressionEnabledField();
385       buildExpression();
386       active=true;
387     }
388   }
389
390   //----------------------------------------------------------------
391
392   void  CField::solveOperation(void)
393   {
394      using namespace func;
395
396      if (!hasOutputFile && !hasFieldOut) return ;
397
398      StdString id ;
399      if (hasId()) id=getId();
400      else if (!name.isEmpty()) id=name ;
401      else if (hasDirectFieldReference()) id=baseRefObject->getId() ;
402
403      CContext* context = CContext::getCurrent();
404
405      if (freq_op.isEmpty()) freq_op.setValue(TimeStep);
406
407      if (operation.isEmpty() )
408      {
409         ERROR("CField::solveOperation(void)",
410               << "[ id = " << id << "]"
411               << "Impossible to define an operation for this field !");
412      }
413
414      if (freq_offset.isEmpty())
415        freq_offset.setValue(NoneDu);
416
417//      if (CXIOSManager::GetStatus() == CXIOSManager::LOC_SERVER)
418      if (context->hasServer)
419      {
420         if (hasOutputFile)
421         {
422           this->freq_operation_srv = this->file->output_freq.getValue();
423           this->freq_write_srv = this->file->output_freq.getValue();
424         }
425         this->lastlast_Write_srv = boost::shared_ptr<CDate>
426                        (new CDate(context->getCalendar()->getInitDate()));
427         this->last_Write_srv     = boost::shared_ptr<CDate>
428                        (new CDate(context->getCalendar()->getInitDate()));
429         this->last_operation_srv = boost::shared_ptr<CDate>
430                        (new CDate(context->getCalendar()->getInitDate()));
431//         this->foperation_srv     =
432//             boost::shared_ptr<func::CFunctor>(new CInstant(this->data_srv));
433
434         if (hasOutputFile)
435         {
436           const CDuration toffset = this->freq_operation_srv - freq_offset.getValue() - context->getCalendar()->getTimeStep();
437           *this->last_operation_srv   = *this->last_operation_srv - toffset;
438         }
439      }
440
441//      if (context->hasClient)
442//      {
443         this->freq_operation = freq_op.getValue();
444         if (hasOutputFile) this->freq_write = this->file->output_freq.getValue();
445         if (hasFieldOut)
446         {
447           this->freq_write = this->fieldOut->freq_op.getValue();
448         }
449         this->last_Write     = boost::shared_ptr<CDate>
450                        (new CDate(context->getCalendar()->getInitDate()));
451         this->last_operation = boost::shared_ptr<CDate>
452                        (new CDate(context->getCalendar()->getInitDate()));
453
454         const CDuration toffset = this->freq_operation - freq_offset.getValue() - context->getCalendar()->getTimeStep();
455         *this->last_operation   = *this->last_operation - toffset;
456
457        if (operation.get()=="once") isOnceOperation=true ;
458        else isOnceOperation=false;
459        isFirstOperation=true;
460
461
462#define DECLARE_FUNCTOR(MType, mtype)              \
463   if  (operation.getValue().compare(#mtype) == 0) \
464   {                                               \
465      if (!detect_missing_value.isEmpty() && !default_value.isEmpty() && detect_missing_value==true) \
466      { \
467        boost::shared_ptr<func::CFunctor> foperation_(new C##MType(this->data,default_value)); \
468        this->foperation = foperation_; \
469      } \
470      else \
471      { \
472        boost::shared_ptr<func::CFunctor> foperation_(new C##MType(this->data)); \
473        this->foperation = foperation_;  \
474      } \
475      return;                                      \
476   }
477
478#include "functor_type.conf"
479
480         ERROR("CField::solveOperation(void)",
481               << "[ operation = " << operation.getValue() << "]"
482               << "The operation is not defined !");
483//      }
484
485
486   }
487
488   //----------------------------------------------------------------
489/*
490   void CField::fromBinary(StdIStream & is)
491   {
492      SuperClass::fromBinary(is);
493#define CLEAR_ATT(name_)\
494      SuperClassAttribute::operator[](#name_)->reset()
495
496         CLEAR_ATT(domain_ref);
497         CLEAR_ATT(axis_ref);
498#undef CLEAR_ATT
499
500   }
501*/
502   //----------------------------------------------------------------
503
504   void CField::solveGridReference(void)
505   {
506      CDomain* domain;
507      CAxis* axis;
508      std::vector<CDomain*> vecDom;
509      std::vector<CAxis*> vecAxis;
510      std::vector<std::string> domList, axisList;
511
512      if (!domain_ref.isEmpty())
513      {
514         if (CDomain::has(domain_ref.getValue()))
515         {
516           domain = CDomain::get(domain_ref.getValue()) ;
517           vecDom.push_back(domain);
518         }
519         else
520            ERROR("CField::solveGridReference(void)",
521                  << "Reference to the domain \'"
522                  << domain_ref.getValue() << "\' is wrong") ;
523      }
524
525      if (!axis_ref.isEmpty())
526      {
527         if (CAxis::has(axis_ref.getValue()))
528         {
529           axis = CAxis::get(axis_ref.getValue()) ;
530           vecAxis.push_back(axis);
531         }
532         else
533            ERROR("CField::solveGridReference(void)",
534                  << "Reference to the axis \'"
535                  << axis_ref.getValue() <<"\' is wrong") ;
536      }
537
538      if (!grid_ref.isEmpty())
539      {
540         if (CGrid::has(grid_ref.getValue()))
541         {
542           this->grid = CGrid::get(grid_ref.getValue()) ;
543           domList = grid->getDomainList();
544           axisList = grid->getAxisList();
545         }
546         else
547            ERROR("CField::solveGridReference(void)",
548                  << "Reference to the grid \'"
549                  << grid_ref.getValue() << "\' is wrong");
550      }
551
552      if (grid_ref.isEmpty() &&  domain_ref.isEmpty() && axis_ref.isEmpty())
553      {
554            ERROR("CField::solveGridReference(void)",
555                  << "At least one dimension must be defined for this field.");
556      }
557
558//     if (!grid_ref.isEmpty())
559//     {
560//       domain = grid->domain;
561//       axis = grid->axis;
562//     }
563
564//     CType<string> goodDomain ;
565//     CType<string> goodAxis ;
566//     if (!grid_ref.isEmpty())
567//     {
568//       if (!grid->domain_ref.isEmpty()) goodDomain=grid->domain_ref ;
569//       if (!grid->axis_ref.isEmpty()) goodAxis=grid->axis_ref ;
570//     }
571//     if (!domain_ref.isEmpty()) goodDomain=domain_ref ;
572//     if (!axis_ref.isEmpty()) goodAxis=axis_ref ;
573
574//     CArray<std::string,1> domListTmp = grid->domainList.getValue();
575//     CArray<std::string,1> axisListTmp = grid->axisList.getValue();
576
577     if (domList.empty() && axisList.empty())
578     {
579       this->grid = CGrid::createGrid(vecDom, vecAxis);
580     }
581
582//     std::string goodDomain = domListTmp[0];
583//     std::string goodAxis = axisListTmp[0];
584
585//     if (goodDomain.isEmpty())
586//     if (goodDomain.empty())
587//     {
588//       ERROR("CField::solveGridReference(void)", << "The horizontal domain for this field is not defined");
589//     }
590//     else
591//     {
592//       if (CDomain::has(goodDomain)) domain = CDomain::get(goodDomain) ;
593//       else ERROR("CField::solveGridReference(void)",<< "Reference to the domain \'"
594//                  <<goodDomain << "\' is wrong") ;
595////                  <<goodDomain.get() << "\' is wrong") ;
596//     }
597//
598////     if (!goodAxis.isEmpty())
599//     if (!goodAxis.empty())
600//     {
601//       if (CAxis::has(goodAxis))  axis = CAxis::get(goodAxis) ;
602//       else  ERROR("CField::solveGridReference(void)", << "Reference to the axis \'"
603//                  << goodAxis <<"\' is wrong") ;
604//                  << goodAxis.get() <<"\' is wrong") ;
605//     }
606
607//     bool nothingToDo=false ;
608//
609//     if (!grid_ref.isEmpty())
610//     {
611//       if (!grid->domain_ref.isEmpty() && goodDomain.get() == grid->domain_ref.get())
612//         if (goodAxis.isEmpty()) nothingToDo=true ;
613//         else if (!grid->axis_ref.isEmpty())
614//                 if (grid->axis_ref.get()==goodAxis.get()) nothingToDo=true ;
615//     }
616//
617//     nothingToDo = true;
618//     if (!nothingToDo)
619//     {
620//       if (!goodAxis.isEmpty())
621//       {
622//         this->grid = CGrid::createGrid(domain, axis) ;
623//         this->grid_ref.setValue(this->grid->getId());
624//       }
625//       else
626//       {
627//         this->grid = CGrid::createGrid(domain) ;
628//         this->grid_ref.setValue(this->grid->getId());
629//       }
630//     }
631
632//     grid->solveReference() ;
633//     grid->solveDomainAxisRef();
634//     grid->checkMaskIndex();
635   }
636
637   void CField::solveGridDomainAxisRef(bool checkAtt)
638   {
639     grid->solveDomainAxisRef(checkAtt);
640   }
641
642   void CField::solveCheckMaskIndex(bool doSendingIndex)
643   {
644     grid->checkMaskIndex(doSendingIndex);
645   }
646
647   ///-------------------------------------------------------------------
648
649   template <>
650      void CGroupTemplate<CField, CFieldGroup, CFieldAttributes>::solveRefInheritance(void)
651   {
652      if (this->group_ref.isEmpty()) return;
653      StdString gref = this->group_ref.getValue();
654
655      if (!CFieldGroup::has(gref))
656         ERROR("CGroupTemplate<CField, CFieldGroup, CFieldAttributes>::solveRefInheritance(void)",
657               << "[ gref = " << gref << "]"
658               << " invalid group name !");
659
660      CFieldGroup* group = CFieldGroup::get(gref);
661      CFieldGroup* owner = CFieldGroup::get(boost::polymorphic_downcast<CFieldGroup*>(this));
662
663      std::vector<CField*> allChildren  = group->getAllChildren();
664      std::vector<CField*>::iterator
665         it = allChildren.begin(), end = allChildren.end();
666
667      for (; it != end; it++)
668      {
669         CField* child = *it;
670         if (child->hasId()) owner->createChild()->field_ref.setValue(child->getId()) ;
671
672      }
673   }
674
675   void CField::scaleFactorAddOffset(double scaleFactor, double addOffset)
676   {
677     map<int, CArray<double,1>* >::iterator it;
678     for(it=data_srv.begin();it!=data_srv.end();it++) *it->second = (*it->second -addOffset) * 1./scaleFactor  ;
679   }
680
681   void CField::outputField(CArray<double,3>& fieldOut)
682   {
683      map<int, CArray<double,1>* >::iterator it;
684      for(it=data_srv.begin();it!=data_srv.end();it++)
685      {
686        grid->outputField(it->first,*it->second, fieldOut.dataFirst()) ;
687      }
688
689//         grid->outputField(it->first,*it->second, fieldOut.) ;
690
691   }
692
693   void CField::outputField(CArray<double,2>& fieldOut)
694   {
695      map<int, CArray<double,1>* >::iterator it;
696
697      for(it=data_srv.begin();it!=data_srv.end();it++)
698      {
699         grid->outputField(it->first,*it->second, fieldOut) ;
700      }
701   }
702
703   ///-------------------------------------------------------------------
704
705   void CField::parse(xml::CXMLNode & node)
706   {
707      SuperClass::parse(node);
708      if (! node.getContent(this->content))
709      {
710        if (node.goToChildElement())
711        {
712          do
713          {
714            if (node.getElementName()=="variable" || node.getElementName()=="variable_group") this->getVirtualVariableGroup()->parseChild(node);
715          } while (node.goToNextElement()) ;
716          node.goToParentElement();
717        }
718      }
719    }
720
721  CArray<double,1>* CField::getInstantData(void)
722  {
723    if (!hasInstantData)
724    {
725      instantData.resize(grid->storeIndex_client.numElements()) ;
726      hasInstantData=true ;
727    }
728    return &instantData ;
729  }
730
731  void CField::addDependency(CField* field, int slotId)
732  {
733    fieldDependency.push_back(pair<CField*,int>(field,slotId)) ;
734  }
735
736  void CField::buildExpression(void)
737  {
738    if (content.size() > 0)
739    {
740      CSimpleNodeExpr* simpleExpr=parseExpr(content+'\0') ;
741      expression=CFieldNode::newNode(simpleExpr) ;
742      delete simpleExpr ;
743      set<string> instantFieldIds ;
744      map<string,CField*> associatedInstantFieldIds ;
745      expression->getInstantFieldIds(instantFieldIds) ;
746      for (set<string>::iterator it=instantFieldIds.begin() ; it!=instantFieldIds.end();++it)
747      {
748        if (*it!="this")
749        {
750          if (CField::has(*it))
751          {
752            CField* field=CField::get(*it) ;
753//            field->processEnabledField() ;
754            field->buildAllExpressionEnabledField();
755            associatedInstantFieldIds[*it]=field ;
756          }
757          else  ERROR("void CField::buildExpression(void)",<<" Field "<<*it<<" does not exist") ;
758        }
759      }
760
761      set<string> averageFieldIds ;
762      map<string,CField*> associatedAverageFieldIds ;
763
764      expression->getAverageFieldIds(averageFieldIds) ;
765      for (set<string>::iterator it=averageFieldIds.begin() ; it!=averageFieldIds.end();++it)
766      {
767        if (CField::has(*it))
768        {
769           CFieldGroup* root=CFieldGroup::get("field_definition") ;
770           CField* averageField=root->createChild() ;
771           CField* instantField=root->createChild() ;
772           averageField->field_ref=*it ;
773           averageField->hasFieldOut=true ;
774           averageField->fieldOut=instantField ;
775           instantField->freq_op=freq_op ;
776//           averageField-> processEnabledField() ;
777           averageField->buildAllExpressionEnabledField();
778           instantField->SuperClassAttribute::setAttributes(averageField, true);
779           instantField->field_ref.reset() ;
780           instantField->operation.reset() ;
781
782//           instantField-> processEnabledField() ;
783           instantField->buildAllExpressionEnabledField();
784           associatedAverageFieldIds[*it]=instantField  ;
785        }
786        else ERROR("void CField::buildExpression(void)",<<" Field "<<*it<<" does not exist") ;
787      }
788
789      expression->reduce(this,associatedInstantFieldIds,associatedAverageFieldIds) ;
790
791      slots.resize(instantFieldIds.size()+averageFieldIds.size()) ;
792      resetSlots() ;
793      int slotId=0 ;
794      set<CField*> fields ;
795      expression->getFields(fields) ;
796      for (set<CField*>::iterator it=fields.begin() ; it!=fields.end();++it,++slotId) (*it)->addDependency(this,slotId) ;
797      hasExpression=true;
798    }
799  }
800
801  void CField::resetSlots(void)
802  {
803    for(vector<bool>::iterator it=slots.begin();it!=slots.end();++it) *it=false ;
804  }
805
806  bool CField::slotsFull(void)
807  {
808    bool ret=true ;
809    for(vector<bool>::iterator it=slots.begin();it!=slots.end();++it) ret &= *it;
810    return ret ;
811  }
812
813
814  void CField::setSlot(int slotId)
815  {
816    CContext* context = CContext::getCurrent() ;
817    const CDate & currDate = context->getCalendar()->getCurrentDate();
818    if (slotUpdateDate==NULL || currDate!=*slotUpdateDate)
819    {
820      resetSlots() ;
821      if (slotUpdateDate==NULL) slotUpdateDate=new CDate(currDate) ;
822      else *slotUpdateDate=currDate ;
823    }
824    slots[slotId]=true ;
825    if (slotsFull())
826    {
827      CArray<double,1> expr(expression->compute()) ;
828
829      if (hasInstantData)
830      {
831        instantData=expr ;
832        for(list< pair<CField *,int> >::iterator it=fieldDependency.begin(); it!=fieldDependency.end(); ++it)
833          if (it->first!=this) it->first->setSlot(it->second) ;
834      }
835
836      if (hasOutputFile) updateDataFromExpression(expr) ;
837
838    }
839  }
840
841   /*!
842     This function retrieves Id of corresponding domain_ref and axis_ref (if any)
843   of a field. In some cases, only domain exists but axis doesn't
844   \return pair of Domain and Axis id
845   */
846   const std::pair<StdString,StdString>& CField::getDomainAxisIds()
847   {
848     CGrid* cgPtr = getRelGrid();
849     if (NULL != cgPtr)
850     {
851       if (NULL != cgPtr->getRelDomain()) domAxisIds_.first = cgPtr->getRelDomain()->getId();
852       if (NULL != cgPtr->getRelAxis()) domAxisIds_.second = cgPtr->getRelAxis()->getId();
853     }
854
855     return (domAxisIds_);
856   }
857
858   CVariable* CField::addVariable(const string& id)
859   {
860     return vVariableGroup->createChild(id) ;
861   }
862
863   CVariableGroup* CField::addVariableGroup(const string& id)
864   {
865     return vVariableGroup->createChildGroup(id) ;
866   }
867
868
869   void CField::sendAddAllVariables()
870   {
871     if (!getAllVariables().empty())
872     {
873       // Firstly, it's necessary to add virtual variable group
874       sendAddVariableGroup(getVirtualVariableGroup()->getId());
875
876       // Okie, now we can add to this variable group
877       std::vector<CVariable*> allVar = getAllVariables();
878       std::vector<CVariable*>::const_iterator it = allVar.begin();
879       std::vector<CVariable*>::const_iterator itE = allVar.end();
880
881       for (; it != itE; ++it)
882       {
883         this->sendAddVariable((*it)->getId());
884         (*it)->sendAllAttributesToServer();
885         (*it)->sendValue();
886       }
887     }
888   }
889
890   void CField::sendAddVariable(const string& id)
891   {
892    CContext* context=CContext::getCurrent() ;
893
894    if (! context->hasServer )
895    {
896       CContextClient* client=context->client ;
897
898       CEventClient event(this->getType(),EVENT_ID_ADD_VARIABLE) ;
899       if (client->isServerLeader())
900       {
901         CMessage msg ;
902         msg<<this->getId() ;
903         msg<<id ;
904         event.push(client->getServerLeader(),1,msg) ;
905         client->sendEvent(event) ;
906       }
907       else client->sendEvent(event) ;
908    }
909
910   }
911
912
913   void CField::sendAddVariableGroup(const string& id)
914   {
915    CContext* context=CContext::getCurrent() ;
916    if (! context->hasServer )
917    {
918       CContextClient* client=context->client ;
919
920       CEventClient event(this->getType(),EVENT_ID_ADD_VARIABLE_GROUP) ;
921       if (client->isServerLeader())
922       {
923         CMessage msg ;
924         msg<<this->getId() ;
925         msg<<id ;
926         event.push(client->getServerLeader(),1,msg) ;
927         client->sendEvent(event) ;
928       }
929       else client->sendEvent(event) ;
930    }
931
932   }
933
934   void CField::recvAddVariable(CEventServer& event)
935   {
936
937      CBufferIn* buffer=event.subEvents.begin()->buffer;
938      string id;
939      *buffer>>id ;
940      get(id)->recvAddVariable(*buffer) ;
941   }
942
943
944   void CField::recvAddVariable(CBufferIn& buffer)
945   {
946      string id ;
947      buffer>>id ;
948      addVariable(id) ;
949   }
950
951   void CField::recvAddVariableGroup(CEventServer& event)
952   {
953
954      CBufferIn* buffer=event.subEvents.begin()->buffer;
955      string id;
956      *buffer>>id ;
957      get(id)->recvAddVariableGroup(*buffer) ;
958   }
959
960
961   void CField::recvAddVariableGroup(CBufferIn& buffer)
962   {
963      string id ;
964      buffer>>id ;
965      addVariableGroup(id) ;
966   }
967
968   DEFINE_REF_FUNC(Field,field)
969
970//  void CField::addReference(CField* field)
971//  {
972//    refObject.push_back(field) ;
973//  }
974//
975//   //----------------------------------------------------------------
976//
977//   bool CField::hasDirectFieldReference(void) const
978//   {
979//     return (!this->field_ref.isEmpty());
980//   }
981//
982//   //----------------------------------------------------------------
983//
984//   const StdString & CField::getBaseFieldId(void) const
985//   {
986//      return (this->getBaseFieldReference()->getId());
987//   }
988//
989//
990//   //----------------------------------------------------------------
991//
992//   /*!
993//   \brief Get pointer to direct field to which the current field refers.
994//   */
995//   CField* CField::getDirectFieldReference(void) const
996//   {
997//      if (this->field_ref.isEmpty())
998//         return (this->getBaseFieldReference());
999//
1000//      if (! CField::has(this->field_ref.getValue()))
1001//         ERROR("CField::getDirectFieldReference(void)",
1002//               << "[ ref_name = " << this->field_ref.getValue() << "]"
1003//               << " invalid field name !");
1004//
1005//      return (CField::get(this->field_ref.getValue()));
1006//   }
1007//
1008//   //----------------------------------------------------------------
1009//
1010//   CField* CField::getBaseFieldReference(void) const
1011//   {
1012//      return (baseRefObject);
1013//   }
1014//
1015//   //----------------------------------------------------------------
1016//
1017//   const std::vector<CField*>& CField::getAllReference(void) const
1018//   {
1019//      return (refObject);
1020//   }
1021//
1022//
1023//   /*!
1024//   \brief Searching for all reference of a field
1025//   If a field refers to (an)other field(s), we will search for all its referenced parents.
1026//   Moreover, if any father, direct or indirect (e.g: two levels up), has non-empty attributes,
1027//   all its attributes will be added to the current field
1028//   \param [in] apply Flag to specify whether current field uses attributes of its father
1029//               in case the attribute is empty (true) or its attributes are replaced by ones of its father (false)
1030//   */
1031//   void CField::solveRefInheritance(bool apply)
1032//   {
1033//      std::set<CField *> sset;
1034//      CField* refer_sptr;
1035//      CField * refer_ptr = this;
1036//
1037//      while (refer_ptr->hasDirectFieldReference())
1038//      {
1039//         refer_sptr = refer_ptr->getDirectFieldReference();
1040//         refer_ptr  = refer_sptr;
1041//
1042//         if(sset.end() != sset.find(refer_ptr))
1043//         {
1044//            DEBUG (<< "Circular dependency stopped for field object on "
1045//                   << "\"" + refer_ptr->getId() + "\" !");
1046//            break;
1047//         }
1048//
1049//         SuperClassAttribute::setAttributes(refer_ptr, apply);
1050//         sset.insert(refer_ptr);
1051//      }
1052//   }
1053//
1054//   /*!
1055//   \brief Only on SERVER side. Remove all field_ref from current field
1056//   On creating a new field on server side, redundant "field_ref" is still kept in the attribute list
1057//   of the current field. This function removes this from current field
1058//   */
1059//   void CField::removeRefInheritance()
1060//   {
1061//     if (this->field_ref.isEmpty()) return;
1062//     this->clearAttribute("field_ref");
1063//   }
1064//
1065//   void CField::solveBaseReference(void)
1066//   {
1067//      std::set<CField *> sset;
1068//      CField* refer_sptr;
1069//      CField * refer_ptr = this;
1070//
1071//      if (this->hasDirectFieldReference())  baseRefObject = getDirectFieldReference();
1072//      else  baseRefObject = CField::get(this);
1073//
1074//      while (refer_ptr->hasDirectFieldReference())
1075//      {
1076//         refer_sptr = refer_ptr->getDirectFieldReference();
1077//         refer_ptr  = refer_sptr;
1078//
1079//         if(sset.end() != sset.find(refer_ptr))
1080//         {
1081//            DEBUG (<< "Circular dependency stopped for field object on "
1082//                   << "\"" + refer_ptr->getId() + "\" !");
1083//            break;
1084//         }
1085//
1086//         sset.insert(refer_ptr);
1087//      }
1088//
1089//      if (hasDirectFieldReference()) baseRefObject->addReference(this) ;
1090//   }
1091//
1092
1093} // namespace xios
Note: See TracBrowser for help on using the repository browser.