Changeset 2291 for XIOS/dev/dev_ym/XIOS_COUPLING/src/transformation/scalar_algorithm/scalar_algorithm_reduce_axis.cpp

Timestamp:

02/01/22 15:28:48 (2 years ago)

Author:

ymipsl

Message:

Improve reduction transformation

• make the difference between reduction over geometry or reduction between process.

• geometrical reduction :

domain -> axis
axis -> scalar
domain -> scalar

• reduction across processes for redondant geometrical cell :

axis -> axis
scalar -> scalar

Reduction can be local (only for the geometrical cell owned by current process) or global, using the "local" attribute (bool) over the reduction.

YM

File:

• XIOS/dev/dev_ym/XIOS_COUPLING/src/transformation/scalar_algorithm/scalar_algorithm_reduce_axis.cpp (modified) (1 diff)

XIOS/dev/dev_ym/XIOS_COUPLING/src/transformation/scalar_algorithm/scalar_algorithm_reduce_axis.cpp

@@ -79 +79 @@
 
   } 
+
+  bool local=false ;
+  if (!algo->local.isEmpty()) local=algo->local ;
   
-  int globalIndexSize = axisSource-> n_glo;
-  for (int idx = 0; idx < globalIndexSize; ++idx)  transformationMapping_[0].push_back(idx);
+  auto& transMap = this->transformationMapping_;  
+  
+  if (local)
+  {
+    scalarDestination->n=1 ; // no mask
+    scalarDestination->mask.reset() ;
+    
+    CArray<size_t,1> srcGlobalIndex ;
+    axisSource->getLocalView(CElementView::WORKFLOW)->getGlobalIndexView(srcGlobalIndex) ;
+    size_t nbIdx = srcGlobalIndex.numElements();
+    if (nbIdx==0) scalarDestination->n=0 ;
+    scalarDestination->checkAttributes() ;
 
-  scalarDestination->checkAttributes() ;
+    for (size_t idx = 0; idx < nbIdx; ++idx)
+    {
+      size_t globalIdx = srcGlobalIndex(idx);
+      transformationMapping_[0].push_back(globalIdx);
+    }
+
+  }
+  else
+  {
+   scalarDestination->checkAttributes() ;
+
+   int globalIndexSize = axisSource->getLocalView(CElementView::WORKFLOW)->getGlobalSize();
+   CArray<size_t,1> dstGlobalIndex ;
+   scalarDestination->getLocalView(CElementView::WORKFLOW)->getGlobalIndexView(dstGlobalIndex) ;
+   if (dstGlobalIndex.numElements()!=0)
+     for (int idx = 0; idx < globalIndexSize; ++idx)  transformationMapping_[0].push_back(idx);
+  }
+
   this->computeAlgorithm(axisSource->getLocalView(CElementView::WORKFLOW), scalarDestination->getLocalView(CElementView::WORKFLOW)) ;