src/OpenFOAM/matrices/LUscalarMatrix/procLduInterface.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 "procLduInterface.H"
  28 #include "lduInterfaceField.H"
  29 #include "cyclicLduInterface.H"
  30 #include "processorLduInterface.H"
  31
  32 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  33
  34 Foam::procLduInterface::procLduInterface
  35 (
  36     const lduInterfaceField& interface,
  37     const scalarField& coeffs
  38 )
  39 :
  40     faceCells_(interface.interface().faceCells()),
  41     coeffs_(coeffs),
  42     myProcNo_(-1),
  43     neighbProcNo_(-1)
  44 {
  45     if (isA<processorLduInterface>(interface.interface()))
  46     {
  47         const processorLduInterface& pldui =
  48             refCast<const processorLduInterface>(interface.interface());
  49
  50         myProcNo_ = pldui.myProcNo();
  51         neighbProcNo_ = pldui.neighbProcNo();
  52     }
  53     else if (isA<cyclicLduInterface>(interface.interface()))
  54     {
  55     }
  56     else
  57     {
  58         FatalErrorIn
  59         (
  60             "procLduInterface::procLduInterface"
  61             "(const lduInterfaceField&, const scalarField&"
  62         )   << "unknown lduInterface type " << interface.interface().type()
  63             << exit(FatalError);
  64     }
  65 }
  66
  67
  68 Foam::procLduInterface::procLduInterface(Istream& is)
  69 :
  70     faceCells_(is),
  71     coeffs_(is),
  72     myProcNo_(readLabel(is)),
  73     neighbProcNo_(readLabel(is))
  74 {}
  75
  76
  77 // * * * * * * * * * * * * * * * IOstream Operators  * * * * * * * * * * * * //
  78
  79 Foam::Ostream& Foam::operator<<(Ostream& os, const procLduInterface& cldui)
  80 {
  81     os  << cldui.faceCells_
  82         << cldui.coeffs_
  83         << cldui.myProcNo_
  84         << cldui.neighbProcNo_;
  85
  86     return os;
  87 }
  88
  89
  90 // ************************************************************************* //