source: XIOS/trunk/src/transformation/axis_algorithm_interpolate.hpp @ 2033

Last change on this file since 2033 was 2033, checked in by ymipsl, 10 months ago

Add extrapolation to vertical interpolation.
YM

File size: 2.9 KB
Line 
1/*!
2   \file axis_algorithm_interpolate.hpp
3   \author Ha NGUYEN
4   \since 23 June 2015
5   \date 23 June 2015
6
7   \brief Algorithm for interpolation on an axis.
8 */
9#ifndef __XIOS_AXIS_ALGORITHM_INTERPOLATE_HPP__
10#define __XIOS_AXIS_ALGORITHM_INTERPOLATE_HPP__
11
12#include "axis_algorithm_transformation.hpp"
13#include "transformation.hpp"
14
15namespace xios {
16
17class CAxis;
18class CGrid;
19class CInterpolateAxis;
20 
21/*!
22  \class CAxisAlgorithmInterpolate
23  Implementing interpolation on axis
24  The values on axis source are assumed monotonic
25*/
26class CAxisAlgorithmInterpolate : public CAxisAlgorithmTransformation
27{
28public:
29  CAxisAlgorithmInterpolate(CAxis* axisDestination, CAxis* axisSource, CInterpolateAxis* interpAxis);
30
31  virtual ~CAxisAlgorithmInterpolate() {}
32
33  static bool registerTrans();
34
35  virtual StdString getName() {return "interpolate_axis";}
36  bool isInversed() {return isInversed_;}
37
38protected:
39  void computeIndexSourceMapping_(const std::vector<CArray<double,1>* >& dataAuxInputs);
40
41private:
42  void retrieveAllAxisValue(const CArray<double,1>& axisValue, const CArray<bool,1>& axisMask,
43                            std::vector<double>& recvBuff, std::vector<int>& indexVec);
44  void computeInterpolantPoint(const std::vector<double>& recvBuff, const std::vector<int>&, int transPos = 0);
45  void computeInterpolantPoint(const std::vector<double>& recvBuff, const std::vector<int>&, 
46                               const std::vector<CArray<double,1>* >& dataAuxInputs, int transPos = 0);
47  void computeWeightedValueAndMapping(CArray<double,1>& axisDestValue, const std::map<int, std::vector<std::pair<int,double> > >& interpolatingIndexValues, int transPos = 0);
48  void fillInAxisValue(std::vector<CArray<double,1> >& vecAxisValue,
49                       const std::vector<CArray<double,1>* >& dataAuxInputs);
50
51private:
52  // Interpolation order
53  int order_;
54  bool extrapolate_ ;
55  StdString coordinate_;
56  StdString coordinateDST_;
57  bool isInversed_;
58  std::vector<std::vector<int> > transPosition_;
59
60private:
61  static CGenericAlgorithmTransformation* create(CGrid* gridDst, CGrid* gridSrc,
62                                                CTransformation<CAxis>* transformation,
63                                                int elementPositionInGrid,
64                                                std::map<int, int>& elementPositionInGridSrc2ScalarPosition,
65                                                std::map<int, int>& elementPositionInGridSrc2AxisPosition,
66                                                std::map<int, int>& elementPositionInGridSrc2DomainPosition,
67                                                std::map<int, int>& elementPositionInGridDst2ScalarPosition,
68                                                std::map<int, int>& elementPositionInGridDst2AxisPosition,
69                                                std::map<int, int>& elementPositionInGridDst2DomainPosition);
70};
71
72}
73
74#endif // __XIOS_AXIS_ALGORITHM_INTERPOLATE_HPP__
Note: See TracBrowser for help on using the repository browser.