src/finiteVolume/fields/fvPatchFields/constraint/processor/processorFvPatchScalarField.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 "processorFvPatchScalarField.H"
  28
  29 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  30
  31 namespace Foam
  32 {
  33
  34 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  35
  36 template<>
  37 void processorFvPatchField<scalar>::initInterfaceMatrixUpdate
  38 (
  39     const scalarField& psiInternal,
  40     scalarField&,
  41     const lduMatrix&,
  42     const scalarField&,
  43     const direction,
  44     const Pstream::commsTypes commsType
  45 ) const
  46 {
  47     procPatch_.compressedSend
  48     (
  49         commsType,
  50         patch().patchInternalField(psiInternal)()
  51     );
  52 }
  53
  54
  55 template<>
  56 void processorFvPatchField<scalar>::updateInterfaceMatrix
  57 (
  58     const scalarField&,
  59     scalarField& result,
  60     const lduMatrix&,
  61     const scalarField& coeffs,
  62     const direction,
  63     const Pstream::commsTypes commsType
  64 ) const
  65 {
  66     scalarField pnf
  67     (
  68         procPatch_.compressedReceive<scalar>(commsType, this->size())()
  69     );
  70
  71     const unallocLabelList& faceCells = patch().faceCells();
  72
  73     forAll(faceCells, facei)
  74     {
  75         result[faceCells[facei]] -= coeffs[facei]*pnf[facei];
  76     }
  77 }
  78
  79
  80 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  81
  82 } // End namespace Foam
  83
  84 // ************************************************************************* //