00001 /*---------------------------------------------------------------------------*\ 00002 ========= | 00003 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox 00004 \\ / O peration | 00005 \\ / A nd | Copyright (C) 1991-2010 OpenCFD Ltd. 00006 \\/ M anipulation | 00007 ------------------------------------------------------------------------------- 00008 License 00009 This file is part of OpenFOAM. 00010 00011 OpenFOAM is free software: you can redistribute it and/or modify it 00012 under the terms of the GNU General Public License as published by 00013 the Free Software Foundation, either version 3 of the License, or 00014 (at your option) any later version. 00015 00016 OpenFOAM is distributed in the hope that it will be useful, but WITHOUT 00017 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or 00018 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License 00019 for more details. 00020 00021 You should have received a copy of the GNU General Public License 00022 along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>. 00023 00024 Class 00025 Foam::DICPreconditioner 00026 00027 Description 00028 Simplified diagonal-based incomplete Cholesky preconditioner for symmetric 00029 matrices (symmetric equivalent of DILU). The reciprocal of the 00030 preconditioned diagonal is calculated and stored. 00031 00032 SourceFiles 00033 DICPreconditioner.C 00034 00035 \*---------------------------------------------------------------------------*/ 00036 00037 #ifndef DICPreconditioner_H 00038 #define DICPreconditioner_H 00039 00040 #include <OpenFOAM/lduMatrix.H> 00041 00042 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00043 00044 namespace Foam 00045 { 00046 00047 /*---------------------------------------------------------------------------*\ 00048 Class DICPreconditioner Declaration 00049 \*---------------------------------------------------------------------------*/ 00050 00051 class DICPreconditioner 00052 : 00053 public lduMatrix::preconditioner 00054 { 00055 // Private data 00056 00057 //- The reciprocal preconditioned diagonal 00058 scalarField rD_; 00059 00060 00061 public: 00062 00063 //- Runtime type information 00064 TypeName("DIC"); 00065 00066 00067 // Constructors 00068 00069 //- Construct from matrix components and preconditioner solver controls 00070 DICPreconditioner 00071 ( 00072 const lduMatrix::solver&, 00073 const dictionary& solverControlsUnused 00074 ); 00075 00076 00077 // Destructor 00078 00079 virtual ~DICPreconditioner() 00080 {} 00081 00082 00083 // Member Functions 00084 00085 //- Calculate the reciprocal of the preconditioned diagonal 00086 static void calcReciprocalD(scalarField& rD, const lduMatrix& matrix); 00087 00088 //- Return wA the preconditioned form of residual rA 00089 virtual void precondition 00090 ( 00091 scalarField& wA, 00092 const scalarField& rA, 00093 const direction cmpt=0 00094 ) const; 00095 }; 00096 00097 00098 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00099 00100 } // End namespace Foam 00101 00102 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00103 00104 #endif 00105 00106 // ************************ vim: set sw=4 sts=4 et: ************************ //
