Changeset 2385 for XIOS3/trunk/src/distribution/local_view.cpp

Timestamp:

07/26/22 09:50:48 (2 years ago)

Author:

jderouillat

Message:

Move in a new method of CLocalView the management of a without redundancy full view connector (used in hash computation of elements)

File:

• XIOS3/trunk/src/distribution/local_view.cpp (modified) (2 diffs)

XIOS3/trunk/src/distribution/local_view.cpp

@@ -3 +3 @@
 #include "array_new.hpp"
 #include "remote_connector.hpp"
+#include "grid_transform_connector.hpp"
 
 namespace xios
@@ -38 +39 @@
     connector->transferToServer(ind, client, event, message) ;
   }
+  
+  void CLocalView::createWithoutRedundancyFullViewConnector( int globalSize, MPI_Comm comm_file, shared_ptr<CGridTransformConnector>& gridTransformConnector, CArray<size_t,1>& globalIndex )
+  {
+    int comm_file_rank(0);
+    MPI_Comm_rank( comm_file, &comm_file_rank );
+    int comm_file_size(1);
+    MPI_Comm_size( comm_file, &comm_file_size );
 
+    vector<shared_ptr<CLocalView>> srcViews;
+    shared_ptr<CLocalView> srcView = shared_from_this();
+    srcViews.push_back( srcView );
+          
+    // Compute a without redundancy element FULL view to enable a consistent hash computation
+    vector<shared_ptr<CLocalView>> remoteViews;
+    shared_ptr<CLocalView> remoteView;
+    // define the remote view without redundancy (naive distribution of the remote view)
+    int localSize = globalSize/comm_file_size;
+    if ( (comm_file_rank==comm_file_size-1) && (localSize*comm_file_size != globalSize ) )
+      localSize += globalSize-localSize*comm_file_size;
+    globalIndex.resize( localSize );
+    CArray<int,1> index( localSize );
+    for (int iloc=0; iloc<localSize ; iloc++ )
+      {
+        globalIndex(iloc) = comm_file_rank*(globalSize/comm_file_size) + iloc;
+        index(iloc) = iloc;
+      }
+    shared_ptr<CLocalElement> localElement = make_shared<CLocalElement>(comm_file_rank, globalSize, globalIndex) ;
+    localElement->addView(CElementView::FULL, index) ;
+    remoteView = localElement->getView(CElementView::FULL) ;
+    remoteViews.push_back( remoteView );
+        
+    // Compute the connector between current and without redundancy FULL views
+    gridTransformConnector = make_shared<CGridTransformConnector>(srcViews, remoteViews, comm_file ) ;
+    gridTransformConnector->computeConnector(true) ; // eliminateRedondant = true
+
+  }
 
 }