Changeset 2606 for XIOS3/trunk

Timestamp:

01/29/24 10:38:53 (12 months ago)

Author:

jderouillat

Message:

Compute the GridTransformConnector? (and the associated GridRemoteConnector?) only in a distributed context for the identification of templated elements

Location:

XIOS3/trunk/src/node

Files:

• axis.cpp (modified) (1 diff)
• domain.cpp (modified) (1 diff)

XIOS3/trunk/src/node/axis.cpp

@@ -308 +308 @@
    size_t CAxis::computeAttributesHash( MPI_Comm comm )
    {
-     // Compute the hash of distributed attributs (value ...)
-     int globalSize = this->n_glo.getValue();
-     CArray<size_t,1> globalIndex; // No redundancy globalIndex will be computed with the connector
-     shared_ptr<CGridTransformConnector> gridTransformConnector;
-     // Compute a without redundancy element FULL view to enable a consistent hash computation (and a distributed globalIndex)
-     this->getLocalView(CElementView::FULL)->createWithoutRedundancyFullViewConnector( globalSize, comm, gridTransformConnector, globalIndex );
-     int localSize = globalIndex.numElements();
-
-     CArray<double,1> distributedValue ;
-     gridTransformConnector->transfer(this->value, distributedValue );
-        
-     size_t localHash = 0;
-     for (int iloc=0; iloc<localSize ; iloc++ ) localHash+=globalIndex(iloc)*distributedValue(iloc);
+     int sz(1);
+     MPI_Comm_size( comm, &sz );
      size_t distributedHash = 0;
-     MPI_Allreduce( &localHash, &distributedHash, 1, MPI_UNSIGNED_LONG, MPI_SUM, comm  );
+     if (sz!=1) // compute the connector only if the element is distributed
+     {
+       // Compute the hash of distributed attributs (value ...)
+       int globalSize = this->n_glo.getValue();
+       CArray<size_t,1> globalIndex; // No redundancy globalIndex will be computed with the connector
+       shared_ptr<CGridTransformConnector> gridTransformConnector;
+       // Compute a without redundancy element FULL view to enable a consistent hash computation (and a distributed globalIndex)
+       this->getLocalView(CElementView::FULL)->createWithoutRedundancyFullViewConnector( globalSize, comm, gridTransformConnector, globalIndex );
+       int localSize = globalIndex.numElements();
+  
+       CArray<double,1> distributedValue ;
+       gridTransformConnector->transfer(this->value, distributedValue );
+          
+       size_t localHash = 0;
+       for (int iloc=0; iloc<localSize ; iloc++ ) localHash+=globalIndex(iloc)*distributedValue(iloc);
+       distributedHash = 0;
+       MPI_Allreduce( &localHash, &distributedHash, 1, MPI_UNSIGNED_LONG, MPI_SUM, comm  );
+     }
+     else // if the element is not distributed, the local hash is valid
+     {
+       int globalSize = this->n_glo.getValue();
+       int localSize = globalSize;
+       size_t localHash = 0;
+       for (int iloc=0; iloc<localSize ; iloc++ ) localHash+=iloc*this->value(iloc);
+       distributedHash = localHash;
+     }
 
      // Compute the hash of global attributs (unit, prec ...)

XIOS3/trunk/src/node/domain.cpp

@@ -1807 +1807 @@
    size_t CDomain::computeAttributesHash( MPI_Comm comm )
    {
-     // Compute the hash of distributed attributs (value ...)
-     int globalSize = this->ni_glo.getValue()*this->nj_glo.getValue();
-     CArray<size_t,1> globalIndex; // No redundancy globalIndex will be computed with the connector
-     shared_ptr<CGridTransformConnector> gridTransformConnector;
-     // Compute a without redundancy element FULL view to enable a consistent hash computation
-     this->getLocalView(CElementView::FULL)->createWithoutRedundancyFullViewConnector( globalSize, comm, gridTransformConnector, globalIndex );
-     int localSize = globalIndex.numElements();
-           
-     CArray<double,1> lon_distributedValue, lat_distributedValue ;
-     gridTransformConnector->transfer(this->lonvalue, lon_distributedValue );
-     gridTransformConnector->transfer(this->latvalue, lat_distributedValue );
-
-     // Compute the distributed hash (v0) of the element
-     // it will be associated to the default element name (= map key), and to the name really written
-     size_t localHash = 0;
-     for (int iloc=0; iloc<localSize ; iloc++ ) localHash+=globalIndex(iloc)*lon_distributedValue(iloc)*lat_distributedValue(iloc);
+     int sz(1);
+     MPI_Comm_size( comm, &sz );
      size_t distributedHash = 0;
-     MPI_Allreduce( &localHash, &distributedHash, 1, MPI_UNSIGNED_LONG, MPI_SUM, comm  );
-     
+     if (sz!=1) // compute the connector only if the element is distributed
+     {
+       // Compute the hash of distributed attributs (value ...)
+       int globalSize = this->ni_glo.getValue()*this->nj_glo.getValue();
+       CArray<size_t,1> globalIndex; // No redundancy globalIndex will be computed with the connector
+       shared_ptr<CGridTransformConnector> gridTransformConnector;
+       // Compute a without redundancy element FULL view to enable a consistent hash computation
+       this->getLocalView(CElementView::FULL)->createWithoutRedundancyFullViewConnector( globalSize, comm, gridTransformConnector, globalIndex );
+       int localSize = globalIndex.numElements();
+             
+       CArray<double,1> lon_distributedValue, lat_distributedValue ;
+       gridTransformConnector->transfer(this->lonvalue, lon_distributedValue );
+       gridTransformConnector->transfer(this->latvalue, lat_distributedValue );
+  
+       // Compute the distributed hash (v0) of the element
+       // it will be associated to the default element name (= map key), and to the name really written
+       size_t localHash = 0;
+       for (int iloc=0; iloc<localSize ; iloc++ ) localHash+=globalIndex(iloc)*lon_distributedValue(iloc)*lat_distributedValue(iloc);
+       distributedHash = 0;
+       MPI_Allreduce( &localHash, &distributedHash, 1, MPI_UNSIGNED_LONG, MPI_SUM, comm  );
+     }
+     else // if the element is not distributed, the local hash is valid
+     {
+       int globalSize = this->ni_glo.getValue()*this->nj_glo.getValue();
+       int localSize = globalSize;
+       size_t localHash = 0;
+       for (int iloc=0; iloc<localSize ; iloc++ ) localHash+=iloc*this->lonvalue(iloc)*this->latvalue(iloc);
+       distributedHash = localHash;
+     }
+
      // Compute the hash of global attributs (unit, prec ...)
      vector<StdString> excludedAttr;