00001 //=========================================================================== 00002 // GoTools - SINTEF Geometry Tools version 1.1 00003 // 00004 // GoTools module: CORE 00005 // 00006 // Copyright (C) 2000-2007 SINTEF ICT, Applied Mathematics, Norway. 00007 // 00008 // This program is free software; you can redistribute it and/or 00009 // modify it under the terms of the GNU General Public License 00010 // as published by the Free Software Foundation version 2 of the License. 00011 // 00012 // This program is distributed in the hope that it will be useful, 00013 // but WITHOUT ANY WARRANTY; without even the implied warranty of 00014 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 00015 // GNU General Public License for more details. 00016 // 00017 // You should have received a copy of the GNU General Public License 00018 // along with this program; if not, write to the Free Software 00019 // Foundation, Inc., 00020 // 59 Temple Place - Suite 330, 00021 // Boston, MA 02111-1307, USA. 00022 // 00023 // Contact information: E-mail: tor.dokken@sintef.no 00024 // SINTEF ICT, Department of Applied Mathematics, 00025 // P.O. Box 124 Blindern, 00026 // 0314 Oslo, Norway. 00027 // 00028 // Other licenses are also available for this software, notably licenses 00029 // for: 00030 // - Building commercial software. 00031 // - Building software whose source code you wish to keep private. 00032 //=========================================================================== 00033 #ifndef _LUDECOMP_NEW_H 00034 #define _LUDECOMP_NEW_H 00035 00036 #include <vector> 00037 #include <stdexcept> 00038 00039 namespace Go 00040 { 00043 00044 00053 //=========================================================================== 00054 template<typename SquareMatrix> 00055 void LUDecomp(SquareMatrix& mat, int num_rows, int* perm, bool& parity); 00056 //=========================================================================== 00057 00069 //=========================================================================== 00070 template<typename SquareMatrix, typename T> 00071 void LUsolveSystem(SquareMatrix& A, int num_unknowns, T* vec); 00072 //=========================================================================== 00073 00083 //=========================================================================== 00084 template<typename SquareMatrix, typename T> 00085 void forwardSubstitution(const SquareMatrix& L, T* x, int num_unknowns); 00086 //=========================================================================== 00087 00097 //=========================================================================== 00098 template<typename SquareMatrix, typename T> 00099 void backwardSubstitution(const SquareMatrix& U, T* x, int num_unknowns); 00100 //=========================================================================== 00101 00111 //=========================================================================== 00112 template<typename SquareMatrix> 00113 void forwardSubstitution(const SquareMatrix& L, std::vector<double>* x, int num_unknowns); 00114 //=========================================================================== 00115 00125 //=========================================================================== 00126 template<typename SquareMatrix> 00127 void backwardSubstitution(const SquareMatrix& U, std::vector<double>* x, int num_unknowns); 00128 //=========================================================================== 00129 00131 }; // end namespace Go 00132 00133 #include "LUDecomp_implementation.h" 00134 00135 #endif // _LUDECOMP_NEW_H 00136 00137