opennl: opennl: nl_preconditioners.h Source File

00001 /*
00002  *  OpenNL: Numerical Library
00003  *  Copyright (C) 2004 Bruno Levy
00004  *
00005  *  This program is free software; you can redistribute it and/or modify
00006  *  it under the terms of the GNU General Public License as published by
00007  *  the Free Software Foundation; either version 2 of the License, or
00008  *  (at your option) any later version.
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., 675 Mass Ave, Cambridge, MA 02139, USA.
00018  *
00019  *  If you modify this software, you should include a notice giving the
00020  *  name of the person performing the modification, the date of modification,
00021  *  and the reason for such modification.
00022  *
00023  *  Contact: Bruno Levy
00024  *
00025  *     levy@loria.fr
00026  *
00027  *     ISA Project
00028  *     LORIA, INRIA Lorraine, 
00029  *     Campus Scientifique, BP 239
00030  *     54506 VANDOEUVRE LES NANCY CEDEX 
00031  *     FRANCE
00032  *
00033  *  Note that the GNU General Public License does not permit incorporating
00034  *  the Software into proprietary programs. 
00035  */
00036 
00037 #ifndef __NL_PRECONDITIONERS__
00038 #define __NL_PRECONDITIONERS__
00039 
00040 #include <NL/nl_private.h>
00041 
00042 /************************************************************************/
00043 /* preconditioners */
00044 void nlPreconditioner_Jacobi(NLdouble* x, NLdouble* y) ;
00045 void nlPreconditioner_SSOR(NLdouble* x, NLdouble* y) ;
00046 
00047 /* Utilities for preconditioners                 */
00048 /* (use the matrix in the current context) */
00049 void nlMultDiagonal(NLdouble* xy, NLdouble omega) ;
00050 void nlMultDiagonalInverse(NLdouble* xy, NLdouble omega) ;
00051 void nlMultLowerInverse(NLdouble* x, NLdouble* y, double omega) ;
00052 void nlMultUpperInverse(NLdouble* x, NLdouble* y, NLdouble omega) ;
00053 
00054 #endif