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 held by original author
   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 template<>
  81 void processorFvPatchField<scalar>::initInterfaceMatrixUpdate
  82 (
  83     const scalarField& psiInternal,
  84     scalarField&,
  85     const BlockLduMatrix<scalar>&,
  86     const CoeffField<scalar>&,
  87     const Pstream::commsTypes commsType
  88 ) const
  89 {
  90     procPatch_.compressedSend
  91     (
  92         commsType,
  93         patch().patchInternalField(psiInternal)()
  94     );
  95 }
  96
  97
  98 template<>
  99 void processorFvPatchField<scalar>::updateInterfaceMatrix
 100 (
 101     const scalarField&,
 102     scalarField& result,
 103     const BlockLduMatrix<scalar>&,
 104     const CoeffField<scalar>& coeffs,
 105     const Pstream::commsTypes commsType
 106 ) const
 107 {
 108     scalarField pnf
 109     (
 110         procPatch_.compressedReceive<scalar>(commsType, this->size())()
 111     );
 112
 113     const unallocLabelList& faceCells = patch().faceCells();
 114     const scalarField& scalarCoeffs = coeffs.asScalar();
 115
 116     forAll(faceCells, facei)
 117     {
 118         result[faceCells[facei]] -= scalarCoeffs[facei]*pnf[facei];
 119     }
 120 }
 121
 122 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 123
 124 } // End namespace Foam
 125
 126 // ************************************************************************* //