src/OpenFOAM/matrices/lduMatrix/solvers/diagonalSolver/diagonalSolver.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2009 OpenCFD Ltd.
   6      \\/     M anipulation  |
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of OpenFOAM.
  10
  11     OpenFOAM is free software; you can redistribute it and/or modify it
  12     under the terms of the GNU General Public License as published by the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     option) any later version.
  15
  16     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
  17     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
  18     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
  19     for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with OpenFOAM; if not, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 \*---------------------------------------------------------------------------*/
  26
  27 #include "diagonalSolver.H"
  28
  29 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  30
  31 namespace Foam
  32 {
  33     defineTypeNameAndDebug(diagonalSolver, 0);
  34 }
  35
  36
  37 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  38
  39 Foam::diagonalSolver::diagonalSolver
  40 (
  41     const word& fieldName,
  42     const lduMatrix& matrix,
  43     const FieldField<Field, scalar>& interfaceBouCoeffs,
  44     const FieldField<Field, scalar>& interfaceIntCoeffs,
  45     const lduInterfaceFieldPtrsList& interfaces,
  46     const dictionary& solverControls
  47 )
  48 :
  49     lduMatrix::solver
  50     (
  51         fieldName,
  52         matrix,
  53         interfaceBouCoeffs,
  54         interfaceIntCoeffs,
  55         interfaces,
  56         solverControls
  57     )
  58 {}
  59
  60
  61 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  62
  63 Foam::lduMatrix::solverPerformance Foam::diagonalSolver::solve
  64 (
  65     scalarField& psi,
  66     const scalarField& source,
  67     const direction cmpt
  68 ) const
  69 {
  70     psi = source/matrix_.diag();
  71
  72     return lduMatrix::solverPerformance
  73     (
  74         typeName,
  75         fieldName_,
  76         0,
  77         0,
  78         0,
  79         true,
  80         false
  81     );
  82 }
  83
  84
  85 // ************************************************************************* //