00001 //=========================================================================== 00002 // SINTEF LSMG library - version 1.1 00003 // 00004 // Copyright (C) 2000-2005 SINTEF ICT, Applied Mathematics, Norway. 00005 // 00006 // This program is free software; you can redistribute it and/or 00007 // modify it under the terms of the GNU General Public License 00008 // as published by the Free Software Foundation version 2 of the License. 00009 // 00010 // This program is distributed in the hope that it will be useful, 00011 // but WITHOUT ANY WARRANTY; without even the implied warranty of 00012 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 00013 // GNU General Public License for more details. 00014 // 00015 // You should have received a copy of the GNU General Public License 00016 // along with this program; if not, write to the Free Software 00017 // Foundation, Inc., 00018 // 59 Temple Place - Suite 330, 00019 // Boston, MA 02111-1307, USA. 00020 // 00021 // Contact information: e-mail: tor.dokken@sintef.no 00022 // SINTEF ICT, Department of Applied Mathematics, 00023 // P.O. Box 124 Blindern, 00024 // 0314 Oslo, Norway. 00025 // 00026 // Other licenses are also available for this software, notably licenses 00027 // for: 00028 // - Building commercial software. 00029 // - Building software whose source code you wish to keep private. 00030 //=========================================================================== 00031 #ifndef _SMOOTHMATRIX_H 00032 #define _SMOOTHMATRIX_H 00033 00038 class SmoothMatrix { 00039 int K_; // order 00040 int n1_, n2_; // ??? NOTE: This is not the matrix dimension, but 00041 double alphaU_; // "alpha" , "dilation" factor (1/scale) 00042 double alphaU3_; // = alphaU_^3 00043 double alphaV_; 00044 double alphaV3_; // = alphaV_^3 00045 static double innerProduct_[3][4]; // [der][trans] 00046 00047 static double innerProductQd0_[3][3]; // [offset-1][trans] 00048 static double innerProductQd1_[3][3]; // [offset-1][trans] 00049 static double innerProductQd2_[3][3]; // [offset-1][trans] 00050 int cutOff(int ij, int rs, int mn) const; 00051 00052 public: 00053 SmoothMatrix(int m, int n, double dU, double dV); 00054 ~SmoothMatrix(){} 00055 00056 // indices runs from 0 00057 double operator()(int i, int j, int r, int s) const; 00058 double norm_l2() const; 00059 }; 00060 00061 #endif
1.5.1