1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
|
// Copyright (C) 2016 EDF
// All Rights Reserved
// This code is published under the GNU Lesser General Public License (GNU LGPL)
#ifndef SPARSEINTERPOLATORSPECTRAL_H
#define SPARSEINTERPOLATORSPECTRAL_H
#include <Eigen/Dense>
#include "StOpt/core/grids/SparseSpaceGrid.h"
#include "StOpt/core/grids/InterpolatorSpectral.h"
/** \file SparseInterpolatorSpectral.h
* \brief Defines an interpolator for a sparse grid : here is a global interpolator, storing the representation of the function
* to interpolate.
* \author Xavier Warin
*/
namespace StOpt
{
/// \class SparseInterpolatorSpectral SparseInterpolatorSpectral.h
/// Sparse interpolation object with any points
/// Templated are the basis functions
class SparseInterpolatorSpectral : public InterpolatorSpectral
{
const SparseSpaceGrid *m_grid ; //< grid used
Eigen::ArrayXd m_hierar ; //< Hierarchized value of the function to interpolate
public :
/** \brief Constructor taking in values on the grid
* \param p_grid is the sparse grid used to interpolate
* \param p_values Function values on the sparse grid
*/
SparseInterpolatorSpectral(const SparseSpaceGrid *p_grid, const Eigen::ArrayXd &p_values) : m_grid(p_grid), m_hierar(p_values)
{
// store hierarchized values
p_grid->toHierarchize(m_hierar);
}
/** \brief Constructor convenient for deserialization
* \param p_hierar Hierarchical values
*/
SparseInterpolatorSpectral(const Eigen::ArrayXd &p_hierar): m_hierar(p_hierar) {}
/** \brief Affect the grid : use for deserialization
* \param p_grid the grid to affect
*/
void setGrid(const SpaceGrid *p_grid)
{
m_grid = static_cast< const SparseSpaceGrid *>(p_grid);
}
/** \brief Get back grid associated to operator
*/
const StOpt::SpaceGrid *getGrid()
{
return static_cast< const SpaceGrid *>(m_grid);
}
/** \brief interpolate
* \param p_point coordinates of the point for interpolation
* \return interpolated value
*/
inline double apply(const Eigen::ArrayXd &p_point) const
{
return m_grid->createInterpolator(p_point)->apply(m_hierar);
}
/** \brief Get back hierarchical values
*/
const Eigen::ArrayXd &getHierar() const
{
return m_hierar;
}
};
}
#endif /* SPARSEINTERPOLATORSPECTRAL_H */
|
