Changeset 1037 for XIOS/dev/branch_yushan/src/transformation

Timestamp:

01/25/17 16:25:17 (7 years ago)

Author:

yushan

Message:

initialize the branch

Location:

XIOS/dev/branch_yushan/src/transformation

Files:

• Functions/average_reduction.cpp (modified) (2 diffs)
• Functions/average_reduction.hpp (modified) (1 diff)
• Functions/extract.cpp (modified) (1 diff)
• Functions/extract.hpp (modified) (1 diff)
• Functions/max_reduction.cpp (modified) (2 diffs)
• Functions/max_reduction.hpp (modified) (1 diff)
• Functions/min_reduction.cpp (modified) (2 diffs)
• Functions/min_reduction.hpp (modified) (1 diff)
• Functions/reduction.hpp (modified) (1 diff)
• Functions/sum_reduction.cpp (modified) (2 diffs)
• Functions/sum_reduction.hpp (modified) (1 diff)
• axis_algorithm_extract_domain.cpp (modified) (1 diff)
• axis_algorithm_inverse.cpp (modified) (2 diffs)
• axis_algorithm_reduce_domain.cpp (modified) (1 diff)
• domain_algorithm_interpolate.cpp (modified) (9 diffs)
• domain_algorithm_interpolate.hpp (modified) (2 diffs)
• generic_algorithm_transformation.cpp (modified) (2 diffs)
• grid_transformation.cpp (modified) (2 diffs)
• scalar_algorithm_extract_axis.cpp (modified) (1 diff)
• scalar_algorithm_reduce_axis.cpp (modified) (1 diff)
• scalar_algorithm_reduce_domain.cpp (modified) (1 diff)

XIOS/dev/branch_yushan/src/transformation/Functions/average_reduction.cpp

@@ -3 +3 @@
    \author Ha NGUYEN
    \since 8 Sep 2016
-   \date 9 Jan 2017
+   \date 8 Sep 2016
 
    \brief average reduction
  */
 #include "average_reduction.hpp"
-#include "utils.hpp"
 
 namespace xios {
@@ -30 +29 @@
                                        const double* dataInput,
                                        CArray<double,1>& dataOut,
-                                       std::vector<bool>& flagInitial,
-                                       const double& defaultValue)
+                                       std::vector<bool>& flagInitial)
 {
   if (resetWeight_) { weights_.resize(flagInitial.size()); weights_ = 1.0; resetWeight_ = false; }
-  bool hasMissingValue = NumTraits<double>::isnan(defaultValue);
 
-  if (hasMissingValue)
+  int nbLocalIndex = localIndex.size();
+  int currentlocalIndex = 0;
+  double currentWeight  = 0.0;
+  for (int idx = 0; idx < nbLocalIndex; ++idx)
   {
-    int nbLocalIndex = localIndex.size();
-    int currentlocalIndex = 0;
-    double currentWeight  = 0.0;
-    for (int idx = 0; idx < nbLocalIndex; ++idx)
+    currentlocalIndex = localIndex[idx].first;
+    currentWeight     = localIndex[idx].second;
+
+    if (flagInitial[currentlocalIndex])
     {
-      currentlocalIndex = localIndex[idx].first;
-      currentWeight     = localIndex[idx].second;
-      if (!NumTraits<double>::isnan(*(dataInput + idx)))
-      {
-        if (flagInitial[currentlocalIndex])
-        {
-          dataOut(currentlocalIndex) = *(dataInput + idx);
-          flagInitial[currentlocalIndex] = false;
-        }
-        else
-        {
-          dataOut(currentlocalIndex)  += *(dataInput + idx);
-          weights_(currentlocalIndex) += 1.0;
-        }
-      }
+      dataOut(currentlocalIndex) = *(dataInput + idx);
+      flagInitial[currentlocalIndex] = false;
     }
-  }
-  else
-  {
-    int nbLocalIndex = localIndex.size();
-    int currentlocalIndex = 0;
-    double currentWeight  = 0.0;
-    for (int idx = 0; idx < nbLocalIndex; ++idx)
+    else
     {
-      currentlocalIndex = localIndex[idx].first;
-      currentWeight     = localIndex[idx].second;
-
-      if (flagInitial[currentlocalIndex])
-      {
-        dataOut(currentlocalIndex) = *(dataInput + idx);
-        flagInitial[currentlocalIndex] = false;
-      }
-      else
-      {
-        dataOut(currentlocalIndex)  += *(dataInput + idx);
-        weights_(currentlocalIndex) += 1.0;
-      }
+      dataOut(currentlocalIndex)  += *(dataInput + idx);
+      weights_(currentlocalIndex) += 1.0;
     }
   }

XIOS/dev/branch_yushan/src/transformation/Functions/average_reduction.hpp

@@ -26 +26 @@
                      const double* dataInput,
                      CArray<double,1>& dataOut,
-                     std::vector<bool>& flagInitial,
-                     const double& defaultValue);
+                     std::vector<bool>& flagInitial);
 
   virtual void updateData(CArray<double,1>& dataOut);

XIOS/dev/branch_yushan/src/transformation/Functions/extract.cpp

@@ -29 +29 @@
                                        const double* dataInput,
                                        CArray<double,1>& dataOut,
-                                       std::vector<bool>& flagInitial,
-                                       const double& defaultValue)
+                                       std::vector<bool>& flagInitial)
 {
   int nbLocalIndex = localIndex.size();

XIOS/dev/branch_yushan/src/transformation/Functions/extract.hpp

@@ -26 +26 @@
                      const double* dataInput,
                      CArray<double,1>& dataOut,
-                     std::vector<bool>& flagInitial,
-                     const double& defaultValue);
+                     std::vector<bool>& flagInitial);
 
   virtual ~CExtractReductionAlgorithm() {}

XIOS/dev/branch_yushan/src/transformation/Functions/max_reduction.cpp

@@ -8 +8 @@
  */
 #include "max_reduction.hpp"
-#include "utils.hpp"
 
 namespace xios {
@@ -30 +29 @@
                                    const double* dataInput,
                                    CArray<double,1>& dataOut,
-                                   std::vector<bool>& flagInitial,
-                                   const double& defaultValue)
+                                   std::vector<bool>& flagInitial)
 {
-  bool hasMissingValue = NumTraits<double>::isnan(defaultValue);
-
-  if (hasMissingValue)
+  int nbLocalIndex = localIndex.size();
+  int currentlocalIndex = 0;
+  double currentWeight  = 0.0;
+  for (int idx = 0; idx < nbLocalIndex; ++idx)
   {
-    int nbLocalIndex = localIndex.size();
-    int currentlocalIndex = 0;    
-    for (int idx = 0; idx < nbLocalIndex; ++idx)
+    currentlocalIndex = localIndex[idx].first;
+    currentWeight     = localIndex[idx].second;
+    if (flagInitial[currentlocalIndex])
     {
-      currentlocalIndex = localIndex[idx].first;      
-      if (!NumTraits<double>::isnan(*(dataInput + idx)))
-      {
-        if (flagInitial[currentlocalIndex])
-        {
-          dataOut(currentlocalIndex) = *(dataInput + idx);
-          flagInitial[currentlocalIndex] = false;
-        }
-        else
-        {
-          dataOut(currentlocalIndex) = std::max(*(dataInput + idx), dataOut(currentlocalIndex));
-        }
-      }
+      dataOut(currentlocalIndex) = *(dataInput + idx);
+      flagInitial[currentlocalIndex] = false;
     }
-  }
-  else
-  {
-    int nbLocalIndex = localIndex.size();
-    int currentlocalIndex = 0;    
-    for (int idx = 0; idx < nbLocalIndex; ++idx)
+    else
     {
-      currentlocalIndex = localIndex[idx].first;      
-      if (flagInitial[currentlocalIndex])
-      {
-        dataOut(currentlocalIndex) = *(dataInput + idx);
-        flagInitial[currentlocalIndex] = false;
-      }
-      else
-      {
-        dataOut(currentlocalIndex) = std::max(*(dataInput + idx), dataOut(currentlocalIndex));
-      }
+      dataOut(currentlocalIndex) = std::max(*(dataInput + idx), dataOut(currentlocalIndex));
     }
   }

XIOS/dev/branch_yushan/src/transformation/Functions/max_reduction.hpp

@@ -26 +26 @@
                      const double* dataInput,
                      CArray<double,1>& dataOut,
-                     std::vector<bool>& flagInitial,
-                     const double& defaultValue);
+                     std::vector<bool>& flagInitial);
 
   virtual ~CMaxReductionAlgorithm() {}

XIOS/dev/branch_yushan/src/transformation/Functions/min_reduction.cpp

@@ -3 +3 @@
    \author Ha NGUYEN
    \since 27 June 2016
-   \date 9 Jan 2017
+   \date 27 June 2016
 
    \brief min reduction
  */
 #include "min_reduction.hpp"
-#include "utils.hpp"
 
 namespace xios {
@@ -30 +29 @@
                                    const double* dataInput,
                                    CArray<double,1>& dataOut,
-                                   std::vector<bool>& flagInitial,
-                                   const double& defaultValue)
+                                   std::vector<bool>& flagInitial)
 {
-  bool hasMissingValue = NumTraits<double>::isnan(defaultValue);
-
-  if (hasMissingValue)
+  int nbLocalIndex = localIndex.size();
+  int currentlocalIndex = 0;
+  for (int idx = 0; idx < nbLocalIndex; ++idx)
   {
-    int nbLocalIndex = localIndex.size();
-    int currentlocalIndex = 0;
-    for (int idx = 0; idx < nbLocalIndex; ++idx)
+    currentlocalIndex = localIndex[idx].first;
+    if (flagInitial[currentlocalIndex])
     {
-      currentlocalIndex = localIndex[idx].first;
-      if (!NumTraits<double>::isnan(*(dataInput + idx)))
-      {
-        if (flagInitial[currentlocalIndex])
-        {
-          dataOut(currentlocalIndex) = *(dataInput + idx);
-          flagInitial[currentlocalIndex] = false;
-        }
-        else
-        {
-          dataOut(currentlocalIndex) = std::min(*(dataInput + idx), dataOut(currentlocalIndex));
-        }
-      }
+      dataOut(currentlocalIndex) = *(dataInput + idx);
+      flagInitial[currentlocalIndex] = false;
     }
-  }
-  else
-  {
-    int nbLocalIndex = localIndex.size();
-    int currentlocalIndex = 0;
-    for (int idx = 0; idx < nbLocalIndex; ++idx)
+    else
     {
-      currentlocalIndex = localIndex[idx].first;
-      if (flagInitial[currentlocalIndex])
-      {
-        dataOut(currentlocalIndex) = *(dataInput + idx);
-        flagInitial[currentlocalIndex] = false;
-      }
-      else
-      {
-        dataOut(currentlocalIndex) = std::min(*(dataInput + idx), dataOut(currentlocalIndex));
-      }
+      dataOut(currentlocalIndex) = std::min(*(dataInput + idx), dataOut(currentlocalIndex));
     }
   }

XIOS/dev/branch_yushan/src/transformation/Functions/min_reduction.hpp

@@ -26 +26 @@
                      const double* dataInput,
                      CArray<double,1>& dataOut,
-                     std::vector<bool>& flagInitial,
-                     const double& defaultValue);
+                     std::vector<bool>& flagInitial);
 
   virtual ~CMinReductionAlgorithm() {}

XIOS/dev/branch_yushan/src/transformation/Functions/reduction.hpp

@@ -45 +45 @@
                      const double* dataInput,
                      CArray<double,1>& dataOut,
-                     std::vector<bool>& flagInitial,
-                     const double& defaultValue) = 0;
+                     std::vector<bool>& flagInitial) = 0;
   /*!
     Update local data 

XIOS/dev/branch_yushan/src/transformation/Functions/sum_reduction.cpp

@@ -3 +3 @@
    \author Ha NGUYEN
    \since 27 June 2016
-   \date 9 Jan 2017
+   \date 27 June 2016
 
    \brief sum reduction
  */
 #include "sum_reduction.hpp"
-#include "utils.hpp"
 
 namespace xios {
@@ -30 +29 @@
                                    const double* dataInput,
                                    CArray<double,1>& dataOut,
-                                   std::vector<bool>& flagInitial,
-                                   const double& defaultValue)
+                                   std::vector<bool>& flagInitial)
 {
-  bool hasMissingValue = NumTraits<double>::isnan(defaultValue);
-  if (hasMissingValue)
+  int nbLocalIndex = localIndex.size();
+  int currentlocalIndex = 0;
+  double currentWeight  = 0.0;
+  for (int idx = 0; idx < nbLocalIndex; ++idx)
   {
-    int nbLocalIndex = localIndex.size();
-    int currentlocalIndex = 0;    
-    for (int idx = 0; idx < nbLocalIndex; ++idx)
+    currentlocalIndex = localIndex[idx].first;
+    currentWeight     = localIndex[idx].second;
+    if (flagInitial[currentlocalIndex])
     {
-      currentlocalIndex = localIndex[idx].first;   
-      if (!NumTraits<double>::isnan(*(dataInput + idx)))
-      {   
-        if (flagInitial[currentlocalIndex])
-        {
-          dataOut(currentlocalIndex) = *(dataInput + idx);
-          flagInitial[currentlocalIndex] = false;
-        }
-        else
-        {
-          dataOut(currentlocalIndex) += *(dataInput + idx);
-        }
-      }
-    }    
-  }
-  else
-  {
-    int nbLocalIndex = localIndex.size();
-    int currentlocalIndex = 0;    
-    for (int idx = 0; idx < nbLocalIndex; ++idx)
+      dataOut(currentlocalIndex) = *(dataInput + idx);
+      flagInitial[currentlocalIndex] = false;
+    }
+    else
     {
-      currentlocalIndex = localIndex[idx].first;      
-      if (flagInitial[currentlocalIndex])
-      {
-        dataOut(currentlocalIndex) = *(dataInput + idx);
-        flagInitial[currentlocalIndex] = false;
-      }
-      else
-      {
-        dataOut(currentlocalIndex) += *(dataInput + idx);
-      }
+      dataOut(currentlocalIndex) += *(dataInput + idx);
     }
   }

XIOS/dev/branch_yushan/src/transformation/Functions/sum_reduction.hpp

@@ -26 +26 @@
                      const double* dataInput,
                      CArray<double,1>& dataOut,
-                     std::vector<bool>& flagInitial,
-                     const double& defaultValue);
+                     std::vector<bool>& flagInitial);
 
   virtual ~CSumReductionAlgorithm() {}

XIOS/dev/branch_yushan/src/transformation/axis_algorithm_extract_domain.cpp

@@ -71 +71 @@
                                         const double& defaultValue)
 {
-  reduction_->apply(localIndex, dataInput, dataOut, flagInitial, defaultValue);
+  reduction_->apply(localIndex, dataInput, dataOut, flagInitial);
 }
 

XIOS/dev/branch_yushan/src/transformation/axis_algorithm_inverse.cpp

@@ -211 +211 @@
 
   status.resize(requests.size());
+  //printf("(%d) src/transformation/axis_algorithm_inverse 1st waitall\n", clientRank);
   MPI_Waitall(requests.size(), &requests[0], &status[0]);
-
+  //printf("(%d) src/transformation/axis_algorithm_inverse 1st waitall OK\n", clientRank);
 
   std::vector<MPI_Request>().swap(requests);
@@ -245 +246 @@
   }
   status.resize(requests.size());
+  //printf("(%d) src/transformation/axis_algorithm_inverse 2nd waitall\n", clientRank);
+
   MPI_Waitall(requests.size(), &requests[0], &status[0]);
-
+  //printf("(%d) src/transformation/axis_algorithm_inverse 2nd waitall OK\n", clientRank);
 
   size_t nGloAxisDest = axisDest_->n_glo.getValue() - 1;

XIOS/dev/branch_yushan/src/transformation/axis_algorithm_reduce_domain.cpp

@@ -79 +79 @@
                                        const double& defaultValue)
 {
-  reduction_->apply(localIndex, dataInput, dataOut, flagInitial, defaultValue);
+  reduction_->apply(localIndex, dataInput, dataOut, flagInitial);
 }
 

XIOS/dev/branch_yushan/src/transformation/domain_algorithm_interpolate.cpp

@@ -48 +48 @@
 
 CDomainAlgorithmInterpolate::CDomainAlgorithmInterpolate(CDomain* domainDestination, CDomain* domainSource, CInterpolateDomain* interpDomain)
-: CDomainAlgorithmTransformation(domainDestination, domainSource), interpDomain_(interpDomain), writeToFile_(false), readFromFile_(false)
-{
-  CContext* context = CContext::getCurrent();
+: CDomainAlgorithmTransformation(domainDestination, domainSource), interpDomain_(interpDomain), writeToFile_(false)
+{
   interpDomain_->checkValid(domainSource);
-  fileToReadWrite_ = "xios_interpolation_weights_";
-
-  if (interpDomain_->weight_filename.isEmpty())
-  {
-    fileToReadWrite_ += context->getId() + "_" + 
-                    domainSource->getDomainOutputName() + "_" + 
-                    domainDestination->getDomainOutputName() + ".nc";    
-  }
-  else 
-    fileToReadWrite_ = interpDomain_->weight_filename;
-
-  ifstream f(fileToReadWrite_.c_str());  
-  switch (interpDomain_->mode)
-  {
-    case CInterpolateDomain::mode_attr::read:
-      readFromFile_ = true;      
-      break;
-    case CInterpolateDomain::mode_attr::compute:
-      readFromFile_ = false;
-      break;
-    case CInterpolateDomain::mode_attr::read_or_compute:      
-      if (!f.good())
-        readFromFile_ = false;
-      else
-        readFromFile_ = true;
-      break;
-    default:
-      break;
-  } 
-
-  writeToFile_ = interpDomain_->write_weight;  
-    
+  if ((CInterpolateDomain::mode_attr::write == interpDomain_->mode) &&
+      (!interpDomain_->file.isEmpty()))
+    writeToFile_ = true;
 }
 
@@ -384 +354 @@
   }
 
-  if (writeToFile_ && !readFromFile_)
+  if (writeToFile_)
      writeRemapInfo(interpMapValue);
   exchangeRemapInfo(interpMapValue);
@@ -469 +439 @@
 void CDomainAlgorithmInterpolate::computeIndexSourceMapping_(const std::vector<CArray<double,1>* >& dataAuxInputs)
 {
-  if (readFromFile_ && !writeToFile_)  
+  if (CInterpolateDomain::mode_attr::read == interpDomain_->mode)  
     readRemapInfo();
   else
@@ -478 +448 @@
 
 void CDomainAlgorithmInterpolate::writeRemapInfo(std::map<int,std::vector<std::pair<int,double> > >& interpMapValue)
-{  
-  writeInterpolationInfo(fileToReadWrite_, interpMapValue);
+{
+  std::string filename = interpDomain_->file.getValue();
+  writeInterpolationInfo(filename, interpMapValue);
 }
 
 void CDomainAlgorithmInterpolate::readRemapInfo()
-{  
+{
+  std::string filename = interpDomain_->file.getValue();
   std::map<int,std::vector<std::pair<int,double> > > interpMapValue;
-  readInterpolationInfo(fileToReadWrite_, interpMapValue);
+  readInterpolationInfo(filename, interpMapValue);
 
   exchangeRemapInfo(interpMapValue);
@@ -636 +608 @@
     int countBuff = 0;
     MPI_Get_count(&recvStatus, MPI_INT, &countBuff);
+    #ifdef _usingMPI
     clientSrcRank = recvStatus.MPI_SOURCE;
-
+    #elif _usingEP
+    clientSrcRank = recvStatus.ep_src;
+    #endif
     MPI_Recv((recvIndexSrcBuff + receivedSize),
              recvBuffSize,
@@ -681 +656 @@
                                                                 const StdSize size)
 {
-  return CONetCDF4::addDimension(name, size);  
+  CONetCDF4::addDimension(name, size);  
 }
 
@@ -687 +662 @@
                                                                const std::vector<StdString>& dim)
 {
-  return CONetCDF4::addVariable(name, type, dim);
-}
-
-void CDomainAlgorithmInterpolate::WriteNetCdf::endDefinition()
-{
-  CONetCDF4::definition_end();
+  CONetCDF4::addVariable(name, type, dim);
 }
 
@@ -756 +726 @@
   MPI_Scan(&localNbWeight, &startIndex, 1, MPI_LONG, MPI_SUM, client->intraComm);
   
-  if (0 == globalNbWeight)
-  {
-    info << "There is no interpolation weights calculated between "
-         << "domain source: " << domainSrc_->getDomainOutputName()
-         << " and domain destination: " << domainDest_->getDomainOutputName()
-         << std::endl;
-    return;
-  }
-
   std::vector<StdSize> start(1, startIndex - localNbWeight);
   std::vector<StdSize> count(1, localNbWeight);
-  
-  WriteNetCdf netCdfWriter(filename, client->intraComm);  
+
+  WriteNetCdf netCdfWriter(filename, client->intraComm);
+
+  // netCdfWriter = CONetCDF4(filename, false, false, true, client->intraComm, false);
 
   // Define some dimensions
@@ -782 +745 @@
   netCdfWriter.addVariableWrite("weight", NC_DOUBLE, dims);
 
-  // End of definition
-  netCdfWriter.endDefinition();
-
   // // Write variables
-  if (0 != localNbWeight)
-  {
-    netCdfWriter.writeDataIndex(src_idx, "src_idx", false, 0, &start, &count);
-    netCdfWriter.writeDataIndex(dst_idx, "dst_idx", false, 0, &start, &count);
-    netCdfWriter.writeDataIndex(weights, "weight", false, 0, &start, &count);
-  }
+  netCdfWriter.writeDataIndex(src_idx, "src_idx", true, 0, &start, &count);
+  netCdfWriter.writeDataIndex(dst_idx, "dst_idx", true, 0, &start, &count);
+  netCdfWriter.writeDataIndex(weights, "weight", true, 0, &start, &count);
 
   netCdfWriter.closeFile();

XIOS/dev/branch_yushan/src/transformation/domain_algorithm_interpolate.hpp

@@ -48 +48 @@
   CInterpolateDomain* interpDomain_;
   bool writeToFile_;
-  bool readFromFile_;
-  StdString fileToReadWrite_;
 
   // class WriteNetCdf;
@@ -59 +57 @@
     int addVariableWrite(const StdString& name, nc_type type,
                          const std::vector<StdString>& dim);
-    void endDefinition();
     void writeDataIndex(const CArray<int,1>& data, const StdString& name,
                         bool collective, StdSize record,

XIOS/dev/branch_yushan/src/transformation/generic_algorithm_transformation.cpp

@@ -11 +11 @@
 #include "context_client.hpp"
 #include "client_client_dht_template.hpp"
-#include "utils.hpp"
 
 namespace xios {
@@ -32 +31 @@
                                             const double& defaultValue)
 {
-  int nbLocalIndex = localIndex.size();  
-  bool hasMissingValue = NumTraits<double>::isnan(defaultValue);  
+  int nbLocalIndex = localIndex.size();
+  bool hasMissingValue = (0.0 != defaultValue) ? true : false;
   if (hasMissingValue)
   {
     for (int idx = 0; idx < nbLocalIndex; ++idx)
     {
-      if (NumTraits<double>::isnan(*(dataInput + idx)))
+      if (defaultValue == *(dataInput + idx))
       {
         flagInitial[localIndex[idx].first] = false;

XIOS/dev/branch_yushan/src/transformation/grid_transformation.cpp

@@ -531 +531 @@
 
   status.resize(requests.size());
+  //printf("(%d) src/transformation/grid_transformation 1st waitall\n", clientRank);
   MPI_Waitall(requests.size(), &requests[0], &status[0]);
+  //printf("(%d) src/transformation/grid_transformation 1st waitall OK\n", clientRank);
 
   // Okie, now use the mask to identify which index source we need to send, then also signal the destination which masked index we will return
@@ -583 +585 @@
   }
   status.resize(requests.size());
+  //printf("(%d) src/transformation/grid_transformation 2nd waitall\n", clientRank);
   MPI_Waitall(requests.size(), &requests[0], &status[0]);
+  //printf("(%d) src/transformation/grid_transformation 2nd waitall OK\n", clientRank);
 
   // Cool, now we can fill in local index of grid destination (counted for masked index)

XIOS/dev/branch_yushan/src/transformation/scalar_algorithm_extract_axis.cpp

@@ -58 +58 @@
                                          const double& defaultValue)
 {
-  reduction_->apply(localIndex, dataInput, dataOut, flagInitial, defaultValue);
+  reduction_->apply(localIndex, dataInput, dataOut, flagInitial);
 }
 

XIOS/dev/branch_yushan/src/transformation/scalar_algorithm_reduce_axis.cpp

@@ -90 +90 @@
                                          const double& defaultValue)
 {
-  reduction_->apply(localIndex, dataInput, dataOut, flagInitial, defaultValue);
+  reduction_->apply(localIndex, dataInput, dataOut, flagInitial);
 }
 

XIOS/dev/branch_yushan/src/transformation/scalar_algorithm_reduce_domain.cpp

@@ -84 +84 @@
                                          const double& defaultValue)
 {
-  reduction_->apply(localIndex, dataInput, dataOut, flagInitial, defaultValue);
+  reduction_->apply(localIndex, dataInput, dataOut, flagInitial);
 }