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) 2004-2011 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
  13     the Free Software Foundation, either version 3 of the License, or
  14     (at your 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, see <http://www.gnu.org/licenses/>.
  23
  24 \*---------------------------------------------------------------------------*/
  25
  26 #include "procLduInterface.H"
  27 #include "lduInterfaceField.H"
  28 #include "cyclicLduInterface.H"
  29 #include "processorLduInterface.H"
  30
  31 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  32
  33 Foam::procLduInterface::procLduInterface
  34 (
  35     const lduInterfaceField& interface,
  36     const scalarField& coeffs
  37 )
  38 :
  39     faceCells_(interface.interface().faceCells()),
  40     coeffs_(coeffs),
  41     myProcNo_(-1),
  42     neighbProcNo_(-1),
  43     tag_(-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         tag_ = pldui.tag();
  53     }
  54     else if (isA<cyclicLduInterface>(interface.interface()))
  55     {
  56     }
  57     else
  58     {
  59         FatalErrorIn
  60         (
  61             "procLduInterface::procLduInterface"
  62             "(const lduInterfaceField&, const scalarField&"
  63         )   << "Unknown lduInterface type "
  64             << interface.interface().type()
  65             << exit(FatalError);
  66     }
  67 }
  68
  69
  70 Foam::procLduInterface::procLduInterface(Istream& is)
  71 :
  72     faceCells_(is),
  73     coeffs_(is),
  74     myProcNo_(readLabel(is)),
  75     neighbProcNo_(readLabel(is)),
  76     tag_(readLabel(is))
  77 {}
  78
  79
  80 // * * * * * * * * * * * * * * * IOstream Operators  * * * * * * * * * * * * //
  81
  82 Foam::Ostream& Foam::operator<<(Ostream& os, const procLduInterface& cldui)
  83 {
  84     os  << cldui.faceCells_
  85         << cldui.coeffs_
  86         << cldui.myProcNo_
  87         << cldui.neighbProcNo_
  88         << cldui.tag_;
  89
  90     return os;
  91 }
  92
  93
  94 // ************************************************************************* //