IBSimu 1.0.4
|
00001 00005 /* Copyright (c) 2005-2010 Taneli Kalvas. All rights reserved. 00006 * 00007 * You can redistribute this software and/or modify it under the terms 00008 * of the GNU General Public License as published by the Free Software 00009 * Foundation; either version 2 of the License, or (at your option) 00010 * any later version. 00011 * 00012 * This library is distributed in the hope that it will be useful, but 00013 * WITHOUT ANY WARRANTY; without even the implied warranty of 00014 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00015 * General Public License for more details. 00016 * 00017 * You should have received a copy of the GNU General Public License 00018 * along with this library (file "COPYING" included in the package); 00019 * if not, write to the Free Software Foundation, Inc., 51 Franklin 00020 * Street, Fifth Floor, Boston, MA 02110-1301 USA 00021 * 00022 * If you have questions about your rights to use or distribute this 00023 * software, please contact Berkeley Lab's Technology Transfer 00024 * Department at TTD@lbl.gov. Other questions, comments and bug 00025 * reports should be sent directly to the author via email at 00026 * taneli.kalvas@jyu.fi. 00027 * 00028 * NOTICE. This software was developed under partial funding from the 00029 * U.S. Department of Energy. As such, the U.S. Government has been 00030 * granted for itself and others acting on its behalf a paid-up, 00031 * nonexclusive, irrevocable, worldwide license in the Software to 00032 * reproduce, prepare derivative works, and perform publicly and 00033 * display publicly. Beginning five (5) years after the date 00034 * permission to assert copyright is obtained from the U.S. Department 00035 * of Energy, and subject to any subsequent five (5) year renewals, 00036 * the U.S. Government is granted for itself and others acting on its 00037 * behalf a paid-up, nonexclusive, irrevocable, worldwide license in 00038 * the Software to reproduce, prepare derivative works, distribute 00039 * copies to the public, perform publicly and display publicly, and to 00040 * permit others to do so. 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