00001 /***************************************************************************\ 00002 * Copyright (C) by University Paris-Est - MISS team 00003 * SVD.hpp created in 10 2008. 00004 * Mail : biri@univ-mlv.fr 00005 * 00006 * This file is part of the OpenKraken-math. 00007 * 00008 * The OpenKraken-math is free software; you can redistribute it and/or modify 00009 * it under the terms of the GNU Lesser General Public License as published by 00010 * the Free Software Foundation; either version 3 of the License, or 00011 * (at your option) any later version. 00012 * 00013 * The OpenKraken-math is distributed in the hope that it will be useful, 00014 * but WITHOUT ANY WARRANTY; without even the implied warranty of 00015 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 00016 * GNU Lesser General Public License for more details. 00017 * 00018 * You should have received a copy of the GNU Lesser General Public License 00019 * along with this program. If not, see <http://www.gnu.org/licenses/>. 00020 * 00021 \***************************************************************************/ 00022 00023 00024 /* 00025 * Anti-doublon 00026 */ 00027 #ifndef __OPENKN_MATH__SVD_HPP__ 00028 #define __OPENKN_MATH__SVD_HPP__ 00029 00030 #include "Matrix.hpp" 00031 00032 namespace kn{ 00033 00044 void decompositionSVD(Matrix<double>& a, 00045 Vector<double>& w, 00046 Matrix<double>& v); 00047 00048 00057 void sortSingularValuesSVD(Matrix<double>& u, 00058 Vector<double>& w, 00059 Matrix<double>& v); 00060 00071 void solveSVD(const Matrix<double>& u, 00072 const Vector<double>& w, 00073 const Matrix<double>& v, 00074 const Vector<double>& b, 00075 Vector<double>& x); 00076 00077 } 00078 #endif