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 \*---------------------------------------------------------------------------*/ 00025 00026 #include <OpenFOAM/diagonalPreconditioner.H> 00027 00028 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * // 00029 00030 namespace Foam 00031 { 00032 defineTypeNameAndDebug(diagonalPreconditioner, 0); 00033 00034 lduMatrix::preconditioner:: 00035 addsymMatrixConstructorToTable<diagonalPreconditioner> 00036 adddiagonalPreconditionerSymMatrixConstructorToTable_; 00037 00038 lduMatrix::preconditioner:: 00039 addasymMatrixConstructorToTable<diagonalPreconditioner> 00040 adddiagonalPreconditionerAsymMatrixConstructorToTable_; 00041 } 00042 00043 00044 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * // 00045 00046 Foam::diagonalPreconditioner::diagonalPreconditioner 00047 ( 00048 const lduMatrix::solver& sol, 00049 const dictionary& 00050 ) 00051 : 00052 lduMatrix::preconditioner(sol), 00053 rD(sol.matrix().diag().size()) 00054 { 00055 scalar* __restrict__ rDPtr = rD.begin(); 00056 const scalar* __restrict__ DPtr = solver_.matrix().diag().begin(); 00057 00058 register label nCells = rD.size(); 00059 00060 // Generate reciprocal diagonal 00061 for (register label cell=0; cell<nCells; cell++) 00062 { 00063 rDPtr[cell] = 1.0/DPtr[cell]; 00064 } 00065 } 00066 00067 00068 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * // 00069 00070 void Foam::diagonalPreconditioner::precondition 00071 ( 00072 scalarField& wA, 00073 const scalarField& rA, 00074 const direction 00075 ) const 00076 { 00077 scalar* __restrict__ wAPtr = wA.begin(); 00078 const scalar* __restrict__ rAPtr = rA.begin(); 00079 const scalar* __restrict__ rDPtr = rD.begin(); 00080 00081 register label nCells = wA.size(); 00082 00083 for (register label cell=0; cell<nCells; cell++) 00084 { 00085 wAPtr[cell] = rDPtr[cell]*rAPtr[cell]; 00086 } 00087 } 00088 00089 00090 // ************************ vim: set sw=4 sts=4 et: ************************ //
