Ignore:
Timestamp:
04/13/18 16:25:46 (3 years ago)
Author:
yushan
Message:

Branch EP merged with Dev_cmip6 @r1481

File:
1 edited

Legend:

Unmodified
Added
Removed
  • XIOS/dev/branch_openmp/src/transformation/domain_algorithm_interpolate.cpp

    r1460 r1482  
    916916       renormalizationFactor=1 ; 
    917917     } 
     918      
     919    if (firstPass) 
     920    { 
     921      allMissing.resize(dataOut.numElements()) ; 
     922      allMissing=true ; 
     923    } 
    918924 
    919925    for (int idx = 0; idx < nbLocalIndex; ++idx) 
    920926    { 
    921       if (NumTraits<double>::isnan(*(dataInput + idx))) 
     927      if (NumTraits<double>::isNan(*(dataInput + idx))) 
    922928      { 
    923         flagInitial[localIndex[idx].first] = false; 
     929        allMissing(localIndex[idx].first) = allMissing(localIndex[idx].first) && true; 
    924930        if (renormalize) renormalizationFactor(localIndex[idx].first)-=localIndex[idx].second ; 
    925931      } 
    926932      else 
    927933      { 
    928         dataOut(localIndex[idx].first) += *(dataInput + idx) * localIndex[idx].second; 
    929         flagInitial[localIndex[idx].first] = true; // Reset flag to indicate not all data source are nan 
     934        dataOut(localIndex[idx].first) += *(dataInput + idx) * localIndex[idx].second;    
    930935      } 
    931936    } 
    932937 
    933     // If all data source are nan then data destination must be nan 
     938  } 
     939  else 
     940  { 
    934941    for (int idx = 0; idx < nbLocalIndex; ++idx) 
    935942    { 
    936       if (!flagInitial[localIndex[idx].first]) 
    937         dataOut(localIndex[idx].first) = defaultValue; 
    938     } 
    939   } 
    940   else 
    941   { 
    942     for (int idx = 0; idx < nbLocalIndex; ++idx) 
    943     { 
    944943      dataOut(localIndex[idx].first) += *(dataInput + idx) * localIndex[idx].second; 
    945944    } 
     
    949948void CDomainAlgorithmInterpolate::updateData(CArray<double,1>& dataOut) 
    950949{ 
    951   if (detectMissingValue && renormalize) 
    952   { 
    953     if (renormalizationFactor.numElements()>0) dataOut/=renormalizationFactor ; // In some case, process doesn't received any data for interpolation (mask) 
    954                                                                                 // so renormalizationFactor is not initialized 
    955   } 
    956 } 
    957  
    958 } 
     950  if (detectMissingValue) 
     951  { 
     952    double defaultValue = std::numeric_limits<double>::quiet_NaN(); 
     953    size_t nbIndex=dataOut.numElements() ;  
     954 
     955    for (int idx = 0; idx < nbIndex; ++idx) 
     956    { 
     957      if (allMissing(idx)) dataOut(idx) = defaultValue; // If all data source are nan then data destination must be nan 
     958    } 
     959     
     960    if (renormalize) 
     961    { 
     962      if (renormalizationFactor.numElements()>0) dataOut/=renormalizationFactor ; // In some case, process doesn't received any data for interpolation (mask) 
     963                                                                                 // so renormalizationFactor is not initialized 
     964    }     
     965  } 
     966} 
     967 
     968} 
Note: See TracChangeset for help on using the changeset viewer.