Changeset 742 for XIOS/trunk/src/node

Timestamp:

10/15/15 15:46:40 (9 years ago)

Author:

mhnguyen

Message:

Implement direct transformation with domain_ref and axis_ref

+) Add a new case in which transformations among fields can be done via domain_ref and axis_ref (not only grid_ref)
+) Fix a minor bug relating to baseFieldReference

Test
+) On Curie
+) all standard tests (client, complete) pass
+) test_remap pass and the results are correct

Location:

XIOS/trunk/src/node

Files:

• axis.cpp (modified) (3 diffs)
• axis.hpp (modified) (1 diff)
• field.cpp (modified) (4 diffs)
• grid.cpp (modified) (6 diffs)
• grid.hpp (modified) (2 diffs)

XIOS/trunk/src/node/axis.cpp

@@ -66 +66 @@
    /*!
     * Test whether the data defined on the axis can be outputted in a compressed way.
-    * 
+    *
     * \return true if and only if a mask was defined for this axis
     */
@@ -326 +326 @@
    }
 
-   void CAxis::checkAttributesOnClient(const std::vector<int>& globalDim, int orderPositionInGrid,
-                                       CServerDistributionDescription::ServerDistributionType distType)
+   void CAxis::checkAttributesOnClient()
    {
      if (this->areClientAttributesChecked_) return;
@@ -340 +339 @@
                                      CServerDistributionDescription::ServerDistributionType distType)
    {
-     if (!this->areClientAttributesChecked_) checkAttributesOnClient(globalDim,
-                                                                     orderPositionInGrid,
-                                                                     distType);
+     if (!this->areClientAttributesChecked_) checkAttributesOnClient();
      CContext* context = CContext::getCurrent();
 

XIOS/trunk/src/node/axis.hpp

@@ -102 +102 @@
          static void recvServerAttribut(CEventServer& event);
          void recvServerAttribut(CBufferIn& buffer) ;
-         void checkAttributesOnClient(const std::vector<int>& globalDim, int orderPositionInGrid,
-                                      CServerDistributionDescription::ServerDistributionType disType = CServerDistributionDescription::BAND_DISTRIBUTION);
+         void checkAttributesOnClient();
          void sendCheckedAttributes(const std::vector<int>& globalDim, int orderPositionInGrid,
                                     CServerDistributionDescription::ServerDistributionType disType = CServerDistributionDescription::BAND_DISTRIBUTION);

XIOS/trunk/src/node/field.cpp

@@ -755 +755 @@
    void CField::solveGridReference(void)
    {
+      if (grid_ref.isEmpty() && domain_ref.isEmpty() && axis_ref.isEmpty())
+      {
+        ERROR("CField::solveGridReference(void)",
+              << "At least one dimension must be defined for this field.");
+      }
+
       CDomain* domain;
       CAxis* axis;
@@ -761 +767 @@
       std::vector<std::string> domList, axisList;
 
+      // In some cases, domain_ref, axis_ref and grid_ref don't come from the current field but from his direct field_ref
+      // To make sure to use the correct *_ref, we need to compare these values to direct field_ref's
+      CField* baseFieldRef = this->getBaseFieldReference();
+      bool hasDomainFromFieldRef = false;
+      bool hasAxisFromFieldRef = false;
+      bool hasGridFromFieldRef = false;
+
       if (!domain_ref.isEmpty())
       {
+         if (0 != baseFieldRef)
+         {
+           if ((!baseFieldRef->domain_ref.isEmpty()) && (baseFieldRef->domain_ref.getValue() == domain_ref.getValue()))
+            hasDomainFromFieldRef = true;
+         }
+
          if (CDomain::has(domain_ref.getValue()))
          {
            domain = CDomain::get(domain_ref.getValue());
+           // It's a direct reference, so we use this ref as its name
+           domain->name.setValue(domain_ref.getValue());
+
            vecDom.push_back(domain);
          }
@@ -776 +798 @@
       if (!axis_ref.isEmpty())
       {
+         if (0 != baseFieldRef)
+         {
+           if ((!baseFieldRef->axis_ref.isEmpty()) && (baseFieldRef->axis_ref.getValue() == axis_ref.getValue()))
+            hasAxisFromFieldRef = true;
+         }
          if (CAxis::has(axis_ref.getValue()))
          {
            axis = CAxis::get(axis_ref.getValue());
+                      // It's a direct reference, so we use this ref as its name
+           axis->name.setValue(axis_ref.getValue());
+
            vecAxis.push_back(axis);
          }
@@ -789 +819 @@
       if (!grid_ref.isEmpty())
       {
-         if (CGrid::has(grid_ref.getValue()))
+         if (0 != baseFieldRef)
          {
-           this->grid = CGrid::get(grid_ref.getValue());
-           domList = grid->getDomainList();
-           axisList = grid->getAxisList();
-           if (domList.empty() && axisList.empty())
+           if ((!baseFieldRef->grid_ref.isEmpty()) && (baseFieldRef->grid_ref.getValue() == grid_ref.getValue()))
+            hasGridFromFieldRef = true;
+         }
+
+        if ((domain_ref.isEmpty() && axis_ref.isEmpty()) || (!hasGridFromFieldRef))
+        {
+           if (CGrid::has(grid_ref.getValue()))
            {
-             this->grid = CGrid::createGrid(vecDom, vecAxis);
+             this->grid = CGrid::get(grid_ref.getValue());
            }
-         }
-         else
-            ERROR("CField::solveGridReference(void)",
-                  << "Reference to the grid \'"
-                  << grid_ref.getValue() << "\' is wrong");
+           else
+              ERROR("CField::solveGridReference(void)",
+                    << "Reference to the grid \'"
+                    << grid_ref.getValue() << "\' is wrong");
+        }
+        else
+        {
+          this->grid = CGrid::createGrid(vecDom, vecAxis);
+          this->grid_ref.setValue(this->grid->getId());
+        }
       }
       else
       {
          this->grid = CGrid::createGrid(vecDom, vecAxis);
-      }
-
-      if (grid_ref.isEmpty() && domain_ref.isEmpty() && axis_ref.isEmpty())
-      {
-            ERROR("CField::solveGridReference(void)",
-                  << "At least one dimension must be defined for this field.");
+         this->grid_ref.setValue(this->grid->getId());
       }
    }

XIOS/trunk/src/node/grid.cpp

@@ -21 +21 @@
 namespace xios {
 
-   /// ////////////////////// Définitions ////////////////////// ///
+   /// ////////////////////// Dfinitions ////////////////////// ///
 
    CGrid::CGrid(void)
@@ -126 +126 @@
        }
      }
-     
+
      return attributesSizes;
    }
@@ -342 +342 @@
             axisListP[i]->sendCheckedAttributes(globalDim_,axisPositionInGrid_[i]);
           else
-            axisListP[i]->checkAttributesOnClient(globalDim_,axisPositionInGrid_[i]);
+            axisListP[i]->checkAttributesOnClient();
           ++idx;
         }
@@ -500 +500 @@
       }
 
-      grid->computeGridGlobalDimension(domains, axis, grid->axis_domain_order);
+
+      grid->solveDomainAxisRefInheritance(true);
 
       return grid;
@@ -1307 +1308 @@
     if (isDomListSet) return;
     std::vector<CDomain*> domList = this->getVirtualDomainGroup()->getAllChildren();
-    if (!domains.empty() && domList.empty()) domList = domains;
+    if (!domains.empty() && domList.empty())
+    {
+      for (int i = 0; i < domains.size(); ++i)
+        this->getVirtualDomainGroup()->addChild(domains[i]);
+      domList = this->getVirtualDomainGroup()->getAllChildren();
+    }
+
     if (!domList.empty())
     {
@@ -1329 +1336 @@
     if (isAxisListSet) return;
     std::vector<CAxis*> aList = this->getVirtualAxisGroup()->getAllChildren();
-    if (!axis.empty() && aList.empty()) aList = axis;
+    if (!axis.empty() && aList.empty())
+    {
+      for (int i = 0; i < axis.size(); ++i)
+        this->getVirtualAxisGroup()->addChild(axis[i]);
+      aList = this->getVirtualAxisGroup()->getAllChildren();
+    }
+
     if (!aList.empty())
     {

XIOS/trunk/src/node/grid.hpp

@@ -202 +202 @@
          void modifyMask(const CArray<int,1>& indexToModify);
 
+         void computeGridGlobalDimension(const std::vector<CDomain*>& domains,
+                                         const std::vector<CAxis*>& axis,
+                                         const CArray<bool,1>& axisDomainOrder);
+
       private:
        template<int N>
@@ -218 +222 @@
         void setAxisList(const std::vector<CAxis*> axis = std::vector<CAxis*>());
         void setDomainList(const std::vector<CDomain*> domains = std::vector<CDomain*>());
-
-        void computeGridGlobalDimension(const std::vector<CDomain*>& domains,
-                                        const std::vector<CAxis*>& axis,
-                                        const CArray<bool,1>& axisDomainOrder);
 
         CDomain* addDomain(const std::string& id);