source: XIOS/dev/dev_olga/src/functor/maximum.cpp @ 1158

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

Two server levels: merging with trunk r1137.
There are bugs.

  • Property copyright set to
    Software name : XIOS (Xml I/O Server)
    http://forge.ipsl.jussieu.fr/ioserver
    Creation date : January 2009
    Licence : CeCCIL version2
    see license file in root directory : Licence_CeCILL_V2-en.txt
    or http://www.cecill.info/licences/Licence_CeCILL_V2-en.html
    Holder : CEA/LSCE (Laboratoire des Sciences du CLimat et de l'Environnement)
    CNRS/IPSL (Institut Pierre Simon Laplace)
    Project Manager : Yann Meurdesoif
    yann.meurdesoif@cea.fr
File size: 1.5 KB
RevLine 
[266]1#include "maximum.hpp"
[369]2#include "array_new.hpp"
[1158]3#include "utils.hpp"
[266]4
[335]5namespace xios
[266]6{
7   namespace func
8   {
9      /// ////////////////////// Définitions ////////////////////// ///
10
[369]11      CMaximum::CMaximum(CArray<double,1>& doutput)
[266]12         : SuperClass(StdString("maximum"), doutput)
13      { /* Ne rien faire de plus */ }
14
[470]15      CMaximum::CMaximum(CArray<double,1>& doutput, double missingValue)
16         : SuperClass(StdString("maximum"), doutput, missingValue)
17      { /* Ne rien faire de plus */ }
18
[266]19      CMaximum::~CMaximum(void)
20      { /* Ne rien faire de plus */ }
21
22      //---------------------------------------------------------------
23
[369]24      void CMaximum::apply(const CArray<double,1>& _dinput,
25                                 CArray<double,1>& _doutput)
[266]26      {
[369]27        const double * it1  = _dinput.dataFirst(),
28                     * end1 = _dinput.dataFirst() + _dinput.numElements();
29              double * it   = _doutput.dataFirst();
30         if (this->nbcall == 1) for (; it1 != end1; it1++, it++) *it = *it1;
[470]31         else 
32         {
33           if (hasMissingValue) 
34           { 
35             for (; it1 != end1; it1++, it++) 
[1158]36               if (!NumTraits<double>::isnan(*it1))
[470]37               {
[1158]38                 if (!NumTraits<double>::isnan(*it)) *it = std::max(*it1, *it);
[470]39                 else *it=*it1 ; 
40               }
41           }
42           else for (; it1 != end1; it1++, it++) *it = std::max(*it1, *it);
43         }
[286]44
[369]45
[266]46      }
47
48      //---------------------------------------------------------------
49
50   } // namespace func
[591]51} // namespace xios
Note: See TracBrowser for help on using the repository browser.