Changeset 620 for XIOS/trunk/src/node/grid.cpp

Timestamp:

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

Author:

mhnguyen

Message:

Implementing generic transformation algorithm (local commit)

+) Implement 3 important classes:

-gridTransformation to read transformation info from grid and interface with the rest of XIOS
-transformationMapping to be in charge of sending/receiving transformation info among clients
-transformationAlgorithm to represent various algorithms

+) Make some change on field to use the new classes

Test
+) Only test_new_features with inversed axis

File:

• XIOS/trunk/src/node/grid.cpp (modified) (10 diffs)

XIOS/trunk/src/node/grid.cpp

@@ -14 +14 @@
 #include "array_new.hpp"
 #include "client_server_mapping_distributed.hpp"
-#include "invert_algorithm.hpp"
-#include "visitable.hpp"
 
 namespace xios {
@@ -26 +24 @@
       , vDomainGroup_(), vAxisGroup_(), axisList_(), isAxisListSet(false), isDomListSet(false), clientDistribution_(0), isIndexSent(false)
       , serverDistribution_(0), serverDistributionDescription_(0), clientServerMap_(0), writtenDataSize_(0), globalDim_()
-      , connectedDataSize_(), connectedServerRank_(), isDataDistributed_(true), algorithms_(), transformations_(), isTransformed_(false)
+      , connectedDataSize_(), connectedServerRank_(), isDataDistributed_(true), transformations_(0), isTransformed_(false)
+      , axisPositionInGrid_()
    {
      setVirtualDomainGroup();
@@ -37 +36 @@
       , vDomainGroup_(), vAxisGroup_(), axisList_(), isAxisListSet(false), isDomListSet(false), clientDistribution_(0), isIndexSent(false)
       , serverDistribution_(0), serverDistributionDescription_(0), clientServerMap_(0), writtenDataSize_(0), globalDim_()
-      , connectedDataSize_(), connectedServerRank_(), isDataDistributed_(true), algorithms_(), transformations_(), isTransformed_(false)
+      , connectedDataSize_(), connectedServerRank_(), isDataDistributed_(true), transformations_(0), isTransformed_(false)
+      , axisPositionInGrid_()
    {
      setVirtualDomainGroup();
@@ -53 +53 @@
     if (0 != serverDistribution_) delete serverDistribution_;
     if (0 != serverDistributionDescription_) delete serverDistributionDescription_;
-    for (std::vector<CGenericAlgorithm*>::iterator it = algorithms_.begin(); it != algorithms_.end(); ++it) delete *it;
+    if (0 != transformations_) delete transformations_;
    }
 
@@ -324 +324 @@
       {
         int idx = 0;
-        std::vector<int> axisPositionMap;
+        axisPositionInGrid_.resize(0);
         for (int i = 0; i < axis_domain_order.numElements(); ++i)
         {
           if (false == axis_domain_order(i))
           {
-            axisPositionMap.push_back(idx);
+            axisPositionInGrid_.push_back(idx);
             ++idx;
           }
@@ -339 +339 @@
         {
           if (sendAtt)
-            axisListP[i]->sendCheckedAttributes(globalDim_,axisPositionMap[i]);
+            axisListP[i]->sendCheckedAttributes(globalDim_,axisPositionInGrid_[i]);
           else
-            axisListP[i]->checkAttributesOnClient(globalDim_,axisPositionMap[i]);
+            axisListP[i]->checkAttributesOnClient(globalDim_,axisPositionInGrid_[i]);
           ++idx;
         }
@@ -347 +347 @@
       }
    }
+
+   std::vector<int> CGrid::getAxisPositionInGrid() const
+   {
+     return axisPositionInGrid_;
+   }
+
 
    //---------------------------------------------------------------
@@ -759 +765 @@
   {
     return serverDistribution_;
+  }
+
+  const CDistributionClient* CGrid::getDistributionClient() const
+  {
+    return clientDistribution_;
   }
 
@@ -1103 +1114 @@
   }
 
-  CGrid::EElementType CGrid::getGridElementType()
-  {
-    EElementType gridType = GRID_ONLY_AXIS;
-    int numElements = axis_domain_order.numElements();
-    if (1 == numElements)
-    {
-      if (true == axis_domain_order(0)) gridType = GRID_ONLY_DOMAIN;
-    }
-    else gridType = GRID_AXIS_DOMAIN;
-  }
-
-  void CGrid::solveTransformations()
-  {
-    std::vector<CAxis*> axisPtr = getAxis();
-    for (std::vector<CAxis*>::iterator it = axisPtr.begin(); it != axisPtr.end(); ++it)
-    {
-      std::vector<ETransformationType> axisTransformation = (*it)->getTransformations();
-      std::vector<ETransformationType>::iterator itTrans  = axisTransformation.begin(),
-                                                 iteTrans = axisTransformation.end();
-      for (;itTrans != iteTrans; ++ itTrans)
-        transformations_.push_back(*itTrans);
-    }
-  }
-
-  void CGrid::setTransformationAlgorithms()
-  {
-    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);
-    }
-  }
-
-  std::vector<ETransformationType> CGrid::getTransformations()
+  bool CGrid::isTransformed()
+  {
+    return isTransformed_;
+  }
+
+  void CGrid::setTransformed()
+  {
+    isTransformed_ = true;
+  }
+
+  CGridTransformation* CGrid::getTransformations()
   {
     return transformations_;
-  }
-
-  bool CGrid::isTransformed()
-  {
-    return isTransformed_;
-  }
-
-  void CGrid::setTransformed()
-  {
-    isTransformed_ = true;
-  }
-
-  const std::vector<CGenericAlgorithm*>& CGrid::getTransformationAlgo()
-  {
-    return algorithms_;
   }
 
@@ -1189 +1150 @@
     }
 
-    EElementType gridType = getGridElementType();
-    CGenericTransformation* gTransform = 0;
-    switch (gridType) {
-    case GRID_ONLY_AXIS:
-      gTransform = new CAxisTransformation(this->getAxis(), transformedGrid->getAxis());
-      break;
-    case GRID_ONLY_DOMAIN:
-      break;
-    case GRID_AXIS_DOMAIN:
-      break;
-    default:
-      break;
-    }
-
-    transformedGrid->solveTransformations();
-    transformedGrid->setTransformationAlgorithms();
-    gTransform->apply(transformedGrid->algorithms_);
-    if (0 != gTransform) delete gTransform;
+    transformations_ = new CGridTransformation(transformedGrid, this);
+    transformations_->computeTransformationMapping();
+    std::cout << "send index " << *(transformations_->getLocalIndexToSendFromGridSource()[0]) << std::endl;
+    std::cout << "receive index " << *(transformations_->getLocalIndexToReceiveOnGridDest()[0][0]) << std::endl;
   }