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

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2008 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 "BICCG.H"
  28
  29 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  30
  31 namespace Foam
  32 {
  33     defineTypeNameAndDebug(BICCG, 0);
  34
  35     lduMatrix::solver::addasymMatrixConstructorToTable<BICCG>
  36         addBICCGSymMatrixConstructorToTable_;
  37 }
  38
  39
  40 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  41
  42 Foam::IStringStream Foam::BICCG::solverDataStream
  43 (
  44     const scalar tolerance,
  45     const scalar relTol
  46 ) const
  47 {
  48     return IStringStream
  49     (
  50         "{ preconditioner   DILU;"
  51         "  tolerance " + name(tolerance) + "; relTol " + name(relTol) + "; }"
  52     );
  53 }
  54
  55
  56 Foam::IStringStream Foam::BICCG::solverDataStream
  57 (
  58     const word& solverName,
  59     Istream& solverData
  60 ) const
  61 {
  62     scalar tolerance(readScalar(solverData));
  63     scalar relTol(readScalar(solverData));
  64
  65     return IStringStream
  66     (
  67         solverName + "{ preconditioner   DILU;"
  68         "  tolerance " + name(tolerance) + "; relTol " + name(relTol) + "; }"
  69     );
  70 }
  71
  72
  73 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  74
  75 Foam::BICCG::BICCG
  76 (
  77     const word& fieldName,
  78     const lduMatrix& matrix,
  79     const FieldField<Field, scalar>& interfaceBouCoeffs,
  80     const FieldField<Field, scalar>& interfaceIntCoeffs,
  81     const lduInterfaceFieldPtrsList& interfaces,
  82     const scalar tolerance,
  83     const scalar relTol
  84 )
  85 :
  86     PBiCG
  87     (
  88         fieldName,
  89         matrix,
  90         interfaceBouCoeffs,
  91         interfaceIntCoeffs,
  92         interfaces,
  93         solverDataStream(tolerance, relTol)()
  94     )
  95 {}
  96
  97
  98 Foam::BICCG::BICCG
  99 (
 100     const word& fieldName,
 101     const lduMatrix& matrix,
 102     const FieldField<Field, scalar>& interfaceBouCoeffs,
 103     const FieldField<Field, scalar>& interfaceIntCoeffs,
 104     const lduInterfaceFieldPtrsList& interfaces,
 105     Istream& solverData
 106 )
 107 :
 108     PBiCG
 109     (
 110         fieldName,
 111         matrix,
 112         interfaceBouCoeffs,
 113         interfaceIntCoeffs,
 114         interfaces,
 115         solverDataStream(word::null, solverData)()
 116     )
 117 {}
 118
 119
 120 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 121
 122 void Foam::BICCG::read(Istream& solverData)
 123 {
 124     word solverName(solverData);
 125     PBiCG::read(solverDataStream(solverName, solverData)());
 126 }
 127
 128
 129 // ************************************************************************* //