matrix_solve.h

Go to the documentation of this file.
00001 /** @file
00002  *  @author Bram de Greve (bramz@users.sourceforge.net)
00003  *  @author Tom De Muer (tomdemuer@users.sourceforge.net)
00004  *
00005  *  *** BEGIN LICENSE INFORMATION ***
00006  *  
00007  *  The contents of this file are subject to the Common Public Attribution License 
00008  *  Version 1.0 (the "License"); you may not use this file except in compliance with 
00009  *  the License. You may obtain a copy of the License at 
00010  *  http://lass.sourceforge.net/cpal-license. The License is based on the 
00011  *  Mozilla Public License Version 1.1 but Sections 14 and 15 have been added to cover 
00012  *  use of software over a computer network and provide for limited attribution for 
00013  *  the Original Developer. In addition, Exhibit A has been modified to be consistent 
00014  *  with Exhibit B.
00015  *  
00016  *  Software distributed under the License is distributed on an "AS IS" basis, WITHOUT 
00017  *  WARRANTY OF ANY KIND, either express or implied. See the License for the specific 
00018  *  language governing rights and limitations under the License.
00019  *  
00020  *  The Original Code is LASS - Library of Assembled Shared Sources.
00021  *  
00022  *  The Initial Developer of the Original Code is Bram de Greve and Tom De Muer.
00023  *  The Original Developer is the Initial Developer.
00024  *  
00025  *  All portions of the code written by the Initial Developer are:
00026  *  Copyright (C) 2004-2007 the Initial Developer.
00027  *  All Rights Reserved.
00028  *  
00029  *  Contributor(s):
00030  *
00031  *  Alternatively, the contents of this file may be used under the terms of the 
00032  *  GNU General Public License Version 2 or later (the GPL), in which case the 
00033  *  provisions of GPL are applicable instead of those above.  If you wish to allow use
00034  *  of your version of this file only under the terms of the GPL and not to allow 
00035  *  others to use your version of this file under the CPAL, indicate your decision by 
00036  *  deleting the provisions above and replace them with the notice and other 
00037  *  provisions required by the GPL License. If you do not delete the provisions above,
00038  *  a recipient may use your version of this file under either the CPAL or the GPL.
00039  *  
00040  *  *** END LICENSE INFORMATION ***
00041  */
00042 
00043 
00044 
00045 #ifndef LASS_GUARDIAN_OF_INCLUSION_NUM_IMPL_MATRIX_SOLVE_H
00046 #define LASS_GUARDIAN_OF_INCLUSION_NUM_IMPL_MATRIX_SOLVE_H
00047 
00048 #include "../num_common.h"
00049 
00050 namespace lass
00051 {
00052 namespace num
00053 {
00054 namespace impl
00055 {
00056 
00057 template 
00058 <
00059     typename T, 
00060     typename RandomIterator1, 
00061     typename RandomIterator2
00062 >
00063 bool ludecomp(RandomIterator1 ioMatrix,
00064               RandomIterator2 oIndex,
00065               size_t iSize,
00066               int& iD);
00067 
00068 
00069 
00070 template 
00071 <
00072     typename T, 
00073     typename RandomIterator1, 
00074     typename RandomIterator2,
00075     typename RandomIterator3
00076 >
00077 void lusolve(RandomIterator1 iMatrix,
00078              RandomIterator2 iIndex,
00079              RandomIterator3 ioColumn,
00080              size_t iSize);
00081 
00082 
00083 
00084 template
00085 <
00086     typename T,
00087     typename RandomIterator1,
00088     typename RandomIterator2,
00089     typename RandomIterator3,
00090     typename RandomIterator4,
00091     typename RandomIterator5
00092 >
00093 void lumprove(RandomIterator1 iMatrix,
00094               RandomIterator2 iMatrixLU,
00095               RandomIterator3 iIndex,
00096               RandomIterator4 iColumn,
00097               RandomIterator5 ioX,
00098               size_t iSize);
00099 
00100 
00101 
00102 template <typename T, typename RandomIterator1, typename RandomIterator2> 
00103 bool cramer2(RandomIterator1 iMatrixRowMajor, 
00104              RandomIterator2 ioColumnFirst, RandomIterator2 ioColumnLast);
00105 
00106 template <typename T, typename RandomIterator1, typename RandomIterator2> 
00107 bool cramer3(RandomIterator1 iMatrixRowMajor,
00108              RandomIterator2 ioColumnFirst, RandomIterator2 ioColumnLast);
00109 
00110 
00111 
00112 
00113 template <typename T, typename RandomIterator1, typename RandomIterator2, typename RandomIterator3>
00114 bool solveTridiagonal(RandomIterator1 iA_1, RandomIterator1 iB_0, RandomIterator1 iC_0,
00115                       RandomIterator2 ioSolution, RandomIterator3 ioTemp, std::size_t iSize);
00116 
00117 }
00118 
00119 }
00120 
00121 }
00122 
00123 #include "matrix_solve.inl"
00124 
00125 #endif
00126 
00127 // EOF
library of assembled shared sources

matrix_solve.h
