Changeset 2663

Timestamp:

10/25/24 08:59:55 (3 months ago)

Author:

jderouillat

Message:

Fix the ordering of the interpolation's weights to ensure reproductibilty of interpolation between compute or read modes

File:

• XIOS2/trunk/src/transformation/domain_algorithm_interpolate.cpp (modified) (2 diffs)

XIOS2/trunk/src/transformation/domain_algorithm_interpolate.cpp

@@ -579 +579 @@
     }      
   }      
+
+  std::map<int, std::vector<double> > transWeightOrdered(transWeight.begin(), transWeight.end());
+  transWeight.clear();
+  std::map<int, std::vector<double> >::const_iterator itWeightO = transWeightOrdered.begin();
+  for(; itWeightO!=transWeightOrdered.end() ; ++itWeightO) {
+      transWeight[itWeightO->first] = itWeightO->second;
+  }
+
+  std::map<int, std::vector<int> > transMapOrdered(transMap.begin(), transMap.end());
+  transMap.clear();
+  std::map<int, std::vector<int> >::const_iterator itMapO = transMapOrdered.begin();
+  for(; itMapO!=transMapOrdered.end() ; ++itMapO) {
+      transMap[itMapO->first] = itMapO->second;
+  }
 }
 CATCH
@@ -790 +804 @@
   }
   
+  std::map<int, std::vector<double> > transWeightOrdered(transWeight.begin(), transWeight.end());
+  transWeight.clear();
+  std::map<int, std::vector<double> >::const_iterator itWeightO = transWeightOrdered.begin();
+  for(; itWeightO!=transWeightOrdered.end() ; ++itWeightO) {
+      transWeight[itWeightO->first] = itWeightO->second;
+  }
+
+  std::map<int, std::vector<int> > transMapOrdered(transMap.begin(), transMap.end());
+  transMap.clear();
+  std::map<int, std::vector<int> >::const_iterator itMapO = transMapOrdered.begin();
+  for(; itMapO!=transMapOrdered.end() ; ++itMapO) {
+      transMap[itMapO->first] = itMapO->second;
+  }
+  
   MPI_Waitall(sendRequest.size(), &sendRequest[0], MPI_STATUS_IGNORE);