00001
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043 #ifndef INTERPOLATION_HPP
00044 #define INTERPOLATION_HPP 1
00045
00046 #include <cstddef>
00047 #include <vector>
00048
00049
00055 class Interpolation2D
00056 {
00057 protected:
00058
00059 size_t _n;
00060 size_t _m;
00061 std::vector<double> _f;
00068 Interpolation2D( size_t n, size_t m, const std::vector<double> &f );
00069
00070 const double &__f( int i, int j ) const;
00071 double &__f( int i, int j );
00072
00073 public:
00074
00077 virtual ~Interpolation2D() {}
00078
00084 virtual double operator()( double x, double y ) const = 0;
00085 };
00086
00087
00088
00089
00090
00096 class ClosestInterpolation2D : public Interpolation2D
00097 {
00098
00099 public:
00100
00103 ClosestInterpolation2D( size_t n, size_t m, const std::vector<double> &f );
00104
00107 virtual ~ClosestInterpolation2D() {}
00108
00114 virtual double operator()( double x, double y ) const;
00115 };
00116
00117
00120 class BiLinearInterpolation2D : public Interpolation2D
00121 {
00122
00123 public:
00124
00127 BiLinearInterpolation2D( size_t n, size_t m, const std::vector<double> &f );
00128
00131 virtual ~BiLinearInterpolation2D() {}
00132
00138 virtual double operator()( double x, double y ) const;
00139 };
00140
00141
00142
00143
00144
00150 class BiCubicInterpolation2D : public Interpolation2D
00151 {
00152
00153 std::vector<double> _fx;
00154 std::vector<double> _fy;
00155 std::vector<double> _fxy;
00156 std::vector<double> _c;
00161 const double &__fx( int i, int j ) const;
00162 const double &__fy( int i, int j ) const;
00163 const double &__fxy( int i, int j ) const;
00164
00165 double &__fx( int i, int j );
00166 double &__fy( int i, int j );
00167 double &__fxy( int i, int j );
00168
00169 static void calc_coefs( double *c, double *x );
00170 static const double wt[16][16];
00171
00172 public:
00173
00176 BiCubicInterpolation2D( size_t n, size_t m, const std::vector<double> &f );
00177
00180 virtual ~BiCubicInterpolation2D() {}
00181
00187 virtual double operator()( double x, double y ) const;
00188 };
00189
00190
00191 #endif
00192
00193
00194
00195
00196
00197
00198
00199
00200
00201
00202
00203
00204
00205
00206
00207
00208
00209
00210