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