Changeset 2386 for XIOS3/trunk/src/io/nc4_data_output.cpp

Timestamp:

07/26/22 15:26:52 (2 years ago)

Author:

jderouillat

Message:

Set the code structure to compute the hash value of an element based on its attributs, use for now before writing an element in a file

File:

• XIOS3/trunk/src/io/nc4_data_output.cpp (modified) (2 diffs)

XIOS3/trunk/src/io/nc4_data_output.cpp

@@ -78 +78 @@
          if (!domain->lonvalue.isEmpty() )
          {
-           int globalSize = domain->ni_glo.getValue()*domain->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
-           domain->getLocalView(CElementView::FULL)->createWithoutRedundancyFullViewConnector( globalSize, comm_file, gridTransformConnector, globalIndex );
-           int localSize = globalIndex.numElements();
-           
-           CArray<double,1> lon_distributedValue, lat_distributedValue ;
-           gridTransformConnector->transfer(domain->lonvalue, lon_distributedValue );
-           gridTransformConnector->transfer(domain->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
-           int localHash = 0;
-           for (int iloc=0; iloc<localSize ; iloc++ ) localHash+=globalIndex(iloc)*lon_distributedValue(iloc)*lat_distributedValue(iloc);
-           int globalHash(0);
-           MPI_Allreduce( &localHash, &globalHash, 1, MPI_INT, MPI_SUM, comm_file  );
-         
+           // The hash of the element will be associated to the default element name (= map key), and to the name really written
+           int globalHash = domain->computeAttributesHash( comm_file ); // Need a MPI_Comm to distribute without redundancy some attributs (value)
+
            StdString defaultNameKey = domain->getDomainOutputName();
            if ( !relDomains_.count ( defaultNameKey ) )
@@ -1056 +1041 @@
         if (!axis->value.isEmpty() )
         {
-           int globalSize = axis->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
-           axis->getLocalView(CElementView::FULL)->createWithoutRedundancyFullViewConnector( globalSize, comm_file, gridTransformConnector, globalIndex );
-           int localSize = globalIndex.numElements();
-
-          CArray<double,1> distributedValue ;
-          gridTransformConnector->transfer(axis->value, 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
-          int localHash = 0;
-          for (int iloc=0; iloc<localSize ; iloc++ ) localHash+=globalIndex(iloc)*distributedValue(iloc);
-          int globalHash(0);
-          MPI_Allreduce( &localHash, &globalHash, 1, MPI_INT, MPI_SUM, comm_file  );
+          // The hash of the element will be associated to the default element name (= map key), and to the name really written
+          int globalHash = axis->computeAttributesHash( comm_file ); // Need a MPI_Comm to distribute without redundancy some attributs (value)
 
           StdString defaultNameKey = axis->getAxisOutputName();