src/foam/matrices/LUscalarMatrix/procLduMatrix.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 "procLduMatrix.H"
  27 #include "procLduInterface.H"
  28 #include "lduMatrix.H"
  29
  30 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  31
  32 Foam::procLduMatrix::procLduMatrix
  33 (
  34     const lduMatrix& ldum,
  35     const FieldField<Field, scalar>& interfaceCoeffs,
  36     const lduInterfaceFieldPtrsList& interfaces
  37 )
  38 :
  39     upperAddr_(ldum.lduAddr().upperAddr()),
  40     lowerAddr_(ldum.lduAddr().lowerAddr()),
  41     diag_(ldum.diag()),
  42     upper_(ldum.upper()),
  43     lower_(ldum.lower())
  44 {
  45     label nInterfaces = 0;
  46
  47     forAll(interfaces, i)
  48     {
  49         if (interfaces.set(i))
  50         {
  51             nInterfaces++;
  52         }
  53     }
  54
  55     interfaces_.setSize(nInterfaces);
  56
  57     nInterfaces = 0;
  58
  59     forAll(interfaces, i)
  60     {
  61         if (interfaces.set(i))
  62         {
  63             interfaces_.set
  64             (
  65                 nInterfaces++,
  66                 new procLduInterface
  67                 (
  68                     interfaces[i],
  69                     interfaceCoeffs[i]
  70                 )
  71             );
  72         }
  73     }
  74
  75 }
  76
  77
  78 Foam::procLduMatrix::procLduMatrix(Istream& is)
  79 :
  80     upperAddr_(is),
  81     lowerAddr_(is),
  82     diag_(is),
  83     upper_(is),
  84     lower_(is),
  85     interfaces_(is)
  86 {}
  87
  88
  89 // * * * * * * * * * * * * * * * IOstream Operators  * * * * * * * * * * * * //
  90
  91 Foam::Ostream& Foam::operator<<(Ostream& os, const procLduMatrix& cldum)
  92 {
  93     os  << cldum.upperAddr_
  94         << cldum.lowerAddr_
  95         << cldum.diag_
  96         << cldum.upper_
  97         << cldum.lower_
  98         << cldum.interfaces_;
  99
 100     return os;
 101 }
 102
 103
 104 // ************************************************************************* //