Changeset 2294

Timestamp:

02/03/22 09:48:58 (2 years ago)

Author:

jderouillat

Message:

Fix for the average reduction in the new implementation of the GatherConnector::transfer

File:

• XIOS/dev/dev_ym/XIOS_COUPLING/src/distribution/gatherer_connector.hpp (modified) (1 diff)

XIOS/dev/dev_ym/XIOS_COUPLING/src/distribution/gatherer_connector.hpp

@@ -218 +218 @@
                     {
                       if (*(count+connector[j])==0) *(output+connector[j]) = *(input + i) ;
-                      else *(output+connector[j]) += *(input + i) ;
+                      else *(output+connector[j]) = *(output+connector[j])* (*(count+connector[j])) + *(input + i) ;
                       (*(count+connector[j]))++ ;
                     } 
                     j++ ;
                   }
+                for(int i=0, j=0; i<srcSize; i++)
+                    if (mask[i]) 
+                    {
+                      if (!std::isnan(*(input + i)))
+                        *(output+connector[j]) /= (*(count+connector[j]));
+                      j++ ;
+                    }
                 break ;
               default :