Changeset 619 for XIOS/trunk/src/node/field.cpp

Timestamp:

06/22/15 13:36:12 (9 years ago)

Author:

mhnguyen

Message:

Implementing the first prototype of filter

+) Create new class filter
+) Implement class for specific algorithm
+) Implement inversing algorithm

Test
+) On Curie
+) Grid with one axis: passed

File:

• XIOS/trunk/src/node/field.cpp (modified) (9 diffs)

XIOS/trunk/src/node/field.cpp

@@ -13 +13 @@
 #include "context_server.hpp"
 #include <set>
+#include "axis_filter.hpp"
+#include "invert_algorithm.hpp"
 
 namespace xios{
@@ -27 +29 @@
       , foperation(), hasInstantData(false), hasExpression(false)
       , active(false) , hasOutputFile(false),hasFieldOut(false), slotUpdateDate(NULL)
-      , processed(false), domAxisIds_("", ""), areAllReferenceSolved(false), areAllExpressionBuilt(false)
+      , processed(false), domAxisIds_("", ""), areAllReferenceSolved(false), areAllExpressionBuilt(false), filter(0)
       , isReadDataRequestPending(false)
       { setVirtualVariableGroup(); }
@@ -40 +42 @@
       , foperation(), hasInstantData(false), hasExpression(false)
       , active(false), hasOutputFile(false), hasFieldOut(false), slotUpdateDate(NULL)
-      , processed(false), domAxisIds_("", ""), areAllReferenceSolved(false), areAllExpressionBuilt(false)
+      , processed(false), domAxisIds_("", ""), areAllReferenceSolved(false), areAllExpressionBuilt(false), filter(0)
       , isReadDataRequestPending(false)
    { setVirtualVariableGroup(); }
@@ -51 +53 @@
       if (hasExpression) delete expression;
       if (slotUpdateDate != NULL) delete slotUpdateDate;
+      if (0 != filter) delete filter;
 
    }
@@ -172 +175 @@
             int rank=(*it).first ;
             CArray<int,1>& index = *(it->second) ;
+                       std::cout << "rank " << index << std::endl;
+                       std::cout << "data " << data << std::endl;
             CArray<double,1> data_tmp(index.numElements()) ;
             for(int n=0;n<data_tmp.numElements();n++) data_tmp(n)=data(index(n)) ;
@@ -527 +532 @@
       return(this->last_operation);
    }
-
-   //----------------------------------------------------------------
-
-//   void CField::processEnabledField(void)
-//   {
-//      if (!processed)
-//      {
-//        processed = true;
-//        solveRefInheritance(true);
-//        solveBaseReference();
-//        solveOperation();
-//        solveGridReference();
-//
-//        if (hasDirectFieldReference()) baseRefObject->processEnabledField();
-//        buildExpression();
-//        active = true;
-//      }
-//    }
 
    void CField::solveAllReferenceEnabledField(bool doSending2Sever)
@@ -562 +549 @@
      }
      solveGridDomainAxisRef(doSending2Sever);
+     if (context->hasClient)
+     {
+       solveTransformedGrid();
+     }
      solveCheckMaskIndex(doSending2Sever);
    }
@@ -576 +567 @@
      {
        areAllExpressionBuilt = true;
-//       solveCheckMaskIndex(true);
-//       solveCheckMaskIndex();
-       if (hasDirectFieldReference()) baseRefObject->buildAllExpressionEnabledField();
+       if (hasDirectFieldReference() && (grid_ref.isEmpty())) baseRefObject->buildAllExpressionEnabledField();
        buildExpression();
        active = true;
@@ -770 +759 @@
    }
 
+   CGrid* CField::getGridRefOfBaseReference()
+   {
+     solveRefInheritance(true);
+     solveBaseReference();
+     baseRefObject->solveGridReference();
+
+     return baseRefObject->grid;
+   }
+
+   void CField::solveTransformedGrid()
+   {
+     if (!grid_ref.isEmpty() && (!field_ref.isEmpty()))
+     {
+       CField* fieldRef  = this;
+       CGrid* gridRefOfFieldRef = 0;
+       while (fieldRef->hasDirectFieldReference())
+       {
+         if ((!(fieldRef->grid_ref.isEmpty())) &&
+             (fieldRef->grid_ref.getValue() != grid_ref.getValue()))
+         {
+           gridRefOfFieldRef = fieldRef->getGridRefOfBaseReference();
+           fieldRef->addReference(this);
+           fieldRef->solveGridDomainAxisRef(false);
+           fieldRef->solveCheckMaskIndex(false);
+           break;
+         }
+         CField* tmp = fieldRef->getDirectFieldReference();
+         fieldRef = tmp;
+       }
+
+       if ((0 == gridRefOfFieldRef) &&
+           (!(fieldRef->grid_ref.isEmpty())) &&
+           (fieldRef->grid_ref.getValue() != grid_ref.getValue()))
+       {
+         gridRefOfFieldRef = fieldRef->getGridRefOfBaseReference();
+         fieldRef->addReference(this);
+         fieldRef->solveGridDomainAxisRef(false);
+         fieldRef->solveCheckMaskIndex(false);
+       }
+
+       CGrid* relGridRef = CGrid::get(grid_ref.getValue());
+       if ((0 != gridRefOfFieldRef) && (relGridRef != gridRefOfFieldRef) && (!(relGridRef->isTransformed())))
+       {
+         gridRefOfFieldRef->transformGrid(relGridRef);
+         filterSources_.push_back(fieldRef);
+         transformations_ = relGridRef->getTransformations();
+         switch (gridRefOfFieldRef->getGridElementType()) {
+         case CGrid::GRID_ONLY_AXIS:
+           filter = new CAxisFilter(gridRefOfFieldRef, relGridRef);
+           break;
+         default:
+           break;
+         }
+         setAlgorithms();
+       }
+     }
+   }
+
+
+  void CField::setAlgorithms()
+  {
+    std::vector<ETransformationType>::iterator itTrans  = transformations_.begin(),
+                                               iteTrans = transformations_.end();
+    std::set<ETransformationType> tmp;
+    for (; itTrans != iteTrans; ++itTrans)
+    {
+      if (tmp.end() == tmp.find(*itTrans))
+      {
+        switch (*itTrans) {
+        case eInverse:
+          algorithms_.push_back(new CInvertAlgorithm());
+          break;
+        default:
+          break;
+        }
+      }
+      tmp.insert(*itTrans);
+    }
+  }
+
+   const std::vector<CField*>& CField::getFilterSources()
+   {
+     return filterSources_;
+   }
+
+   void CField::applyFilter()
+   {
+     filter->setInputs(filterSources_);
+     filter->setOutput(this);
+     filter->apply(algorithms_);
+
+   }
+
    ///-------------------------------------------------------------------