src/ParaMEDMEM/LinearTimeInterpolator.cxx

   1 // Copyright (C) 2007-2024  CEA, EDF
   2 //
   3 // This library is free software; you can redistribute it and/or
   4 // modify it under the terms of the GNU Lesser General Public
   5 // License as published by the Free Software Foundation; either
   6 // version 2.1 of the License, or (at your option) any later version.
   7 //
   8 // This library is distributed in the hope that it will be useful,
   9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11 // Lesser General Public License for more details.
  12 //
  13 // You should have received a copy of the GNU Lesser General Public
  14 // License along with this library; if not, write to the Free Software
  15 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  16 //
  17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  18 //
  19
  20 #include "LinearTimeInterpolator.hxx"
  21
  22 using namespace std;
  23
  24 namespace MEDCoupling
  25 {
  26
  27   LinearTimeInterpolator::LinearTimeInterpolator( double InterpPrecision, int nStepBefore,
  28                                                   int nStepAfter ):
  29     TimeInterpolator( InterpPrecision, nStepBefore, nStepAfter )
  30   {
  31   }
  32
  33   LinearTimeInterpolator::~LinearTimeInterpolator()
  34   {
  35   }
  36
  37   void LinearTimeInterpolator::doInterp( double time0, double time1, double time,
  38                                          int recvcount , int nbuff0, int nbuff1,
  39                                          int **recvbuff0, int **recvbuff1, int *result )
  40   {
  41     for(int i = 0 ; i < recvcount ; i++ )
  42         result[i] = (int) ((recvbuff0[0][i]*(time1 - time) + recvbuff1[0][i]*(time - time0))/(time1 - time0) + _interp_precision);
  43   }
  44
  45   void LinearTimeInterpolator::doInterp( double time0, double time1, double time,
  46                                          int recvcount , int nbuff0, int nbuff1,
  47                                          double **recvbuff0, double **recvbuff1,
  48                                          double *result )
  49   {
  50     for(int i = 0 ; i < recvcount ; i++ )
  51       result[i] = (recvbuff0[0][i]*(time1 - time) + recvbuff1[0][i]*(time - time0))/(time1 - time0);
  52   }
  53
  54 }