source: XIOS/dev/XIOS_DEV_CMIP6/src/transformation/Functions/max_reduction.cpp

Last change on this file was 1474, checked in by oabramkina, 3 years ago

DEV_CMIP6: porting changes in r1471.

Now DEV_CMIP6 compiles with PGI.

File size: 2.0 KB
Line 
1/*!
2   \file max.cpp
3   \author Ha NGUYEN
4   \since 27 June 2016
5   \date 27 June 2016
6
7   \brief max reduction
8 */
9#include "max_reduction.hpp"
10#include "utils.hpp"
11
12namespace xios {
13
14CMaxReductionAlgorithm::CMaxReductionAlgorithm()
15  : CReductionAlgorithm()
16{
17}
18
19CReductionAlgorithm* CMaxReductionAlgorithm::create()
20{
21  return (new CMaxReductionAlgorithm());
22}
23
24bool CMaxReductionAlgorithm::registerTrans()
25{
26  return registerOperation(TRANS_REDUCE_MAX, CMaxReductionAlgorithm::create);
27}
28
29void CMaxReductionAlgorithm::apply(const std::vector<std::pair<int,double> >& localIndex,
30                                   const double* dataInput,
31                                   CArray<double,1>& dataOut,
32                                   std::vector<bool>& flagInitial,
33                                   bool ignoreMissingValue, bool firstPass)
34{ 
35  if (ignoreMissingValue)
36  {
37    int nbLocalIndex = localIndex.size();
38    int currentlocalIndex = 0;
39    if (firstPass) dataOut=std::numeric_limits<double>::quiet_NaN();   
40    for (int idx = 0; idx < nbLocalIndex; ++idx)
41    {
42      currentlocalIndex = localIndex[idx].first;     
43      if (!NumTraits<double>::isNan(*(dataInput + idx)))
44      {
45        if (flagInitial[currentlocalIndex])
46        {
47          dataOut(currentlocalIndex) = *(dataInput + idx);
48          flagInitial[currentlocalIndex] = false;
49        }
50        else
51        {
52          dataOut(currentlocalIndex) = std::max(*(dataInput + idx), dataOut(currentlocalIndex));
53        }
54      }
55    }
56  }
57  else
58  {
59    int nbLocalIndex = localIndex.size();
60    int currentlocalIndex = 0;   
61    for (int idx = 0; idx < nbLocalIndex; ++idx)
62    {
63      currentlocalIndex = localIndex[idx].first;     
64      if (flagInitial[currentlocalIndex])
65      {
66        dataOut(currentlocalIndex) = *(dataInput + idx);
67        flagInitial[currentlocalIndex] = false;
68      }
69      else
70      {
71        dataOut(currentlocalIndex) = std::max(*(dataInput + idx), dataOut(currentlocalIndex));
72      }
73    }
74  }
75}
76
77}
Note: See TracBrowser for help on using the repository browser.