src/finiteVolume/blockLduSystem/BlockLduSystem.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | foam-extend: Open Source CFD
   4    \\    /   O peration     | Version:     3.2
   5     \\  /    A nd           | Web:         http://www.foam-extend.org
   6      \\/     M anipulation  | For copyright notice see file Copyright
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of foam-extend.
  10
  11     foam-extend 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 3 of the License, or (at your
  14     option) any later version.
  15
  16     foam-extend is distributed in the hope that it will be useful, but
  17     WITHOUT ANY WARRANTY; without even the implied warranty of
  18     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  19     General Public License for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 \*---------------------------------------------------------------------------*/
  25
  26 #include "BlockLduSystem.H"
  27 #include "IOstreams.H"
  28
  29 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  30
  31 template<class blockType, class sourceType>
  32 Foam::BlockLduSystem<blockType, sourceType>::BlockLduSystem
  33 (
  34     const lduMesh& ldu
  35 )
  36 :
  37     BlockLduMatrix<blockType>(ldu),
  38     source_(ldu.lduAddr().size(), pTraits<sourceType>::zero)
  39 {}
  40
  41
  42 template<class blockType, class sourceType>
  43 Foam::BlockLduSystem<blockType, sourceType>::BlockLduSystem
  44 (
  45     const lduMesh& ldu,
  46     const Field<sourceType>& s
  47 )
  48 :
  49     BlockLduMatrix<blockType>(ldu),
  50     source_(s)
  51 {
  52     if (ldu.lduAddr().size() != s.size())
  53     {
  54         FatalErrorIn
  55         (
  56             "BlockLduSystem::BlockLduSystem\n"
  57             "(\n"
  58             "    const lduMesh& ldu,"
  59             "    const Field<sourceType>& s,"
  60             ")\n"
  61         )   << "Sizes of ldu addressing and source field are not the same."
  62             << abort(FatalError);
  63     }
  64 }
  65
  66
  67 template<class blockType, class sourceType>
  68 Foam::BlockLduSystem<blockType, sourceType>::BlockLduSystem
  69 (
  70     const BlockLduMatrix<blockType>& bm,
  71     const Field<sourceType>& s
  72 )
  73 :
  74     BlockLduMatrix<blockType>(bm),
  75     source_(s)
  76 {
  77     if (this->lduAddr().size() != s.size())
  78     {
  79         FatalErrorIn
  80         (
  81             "BlockLduSystem::BlockLduSystem\n"
  82             "(\n"
  83             "    const BlockLduMatrix<blockType>& bm,"
  84             "    const Field<sourceType>& s,"
  85             ")\n"
  86         )   << "Sizes of block matrix and source field are not the same."
  87             << abort(FatalError);
  88     }
  89 }
  90
  91
  92 template<class blockType, class sourceType>
  93 Foam::BlockLduSystem<blockType, sourceType>::BlockLduSystem
  94 (
  95     const BlockLduSystem<blockType, sourceType>& bs
  96 )
  97 :
  98     BlockLduMatrix<blockType>(bs),
  99     source_(bs.source())
 100 {}
 101
 102
 103 // * * * * * * * * * * * * * * * IOstream Operators  * * * * * * * * * * * * //
 104
 105 template<class blockType, class sourceType>
 106 Foam::Ostream& Foam::operator<<
 107 (
 108     Ostream& os,
 109     const BlockLduSystem<blockType, sourceType>& bs
 110 )
 111 {
 112     os  << static_cast<const BlockLduMatrix<blockType>&>(bs) << nl
 113         << bs.source() << endl;
 114
 115     return os;
 116 }
 117
 118
 119 // ************************************************************************* //