src/coupledMatrix/coupledLduMatrix/coupledLduMatrix.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2006-7 H. Jasak All rights reserved
   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 Class
  26     coupledLduMatrix
  27
  28 Description
  29     Collection of lduMatrices solved together as a block system
  30
  31 Author
  32     Hrvoje Jasak, Wikki Ltd.  All rights reserved
  33
  34 \*----------------------------------------------------------------------------*/
  35
  36 #include "coupledLduMatrix.H"
  37
  38 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  39
  40 namespace Foam
  41 {
  42     defineTypeNameAndDebug(coupledLduMatrix, 1);
  43 }
  44
  45
  46 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  47
  48
  49 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  50
  51 // Construct given size
  52 Foam::coupledLduMatrix::coupledLduMatrix(const label size)
  53 :
  54     PtrList<lduMatrix>(size)
  55 {}
  56
  57
  58 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
  59
  60 Foam::coupledLduMatrix::~coupledLduMatrix()
  61 {}
  62
  63
  64 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  65
  66 bool Foam::coupledLduMatrix::diagonal() const
  67 {
  68     const PtrList<lduMatrix>& matrices = *this;
  69
  70     bool diag = true;
  71
  72     forAll (matrices, matrixI)
  73     {
  74         diag = diag && matrices[matrixI].diagonal();
  75     }
  76
  77     return diag;
  78 }
  79
  80
  81 bool Foam::coupledLduMatrix::symmetric() const
  82 {
  83     const PtrList<lduMatrix>& matrices = *this;
  84
  85     bool sym = true;
  86
  87     forAll (matrices, matrixI)
  88     {
  89         sym =
  90             (sym && matrices[matrixI].diagonal())
  91          || (sym && matrices[matrixI].symmetric());
  92     }
  93
  94     return sym;
  95 }
  96
  97
  98 bool Foam::coupledLduMatrix::asymmetric() const
  99 {
 100     const PtrList<lduMatrix>& matrices = *this;
 101
 102     bool asym = false;
 103
 104     forAll (matrices, matrixI)
 105     {
 106         asym = (asym || matrices[matrixI].asymmetric());
 107     }
 108
 109     return asym;
 110 }
 111
 112
 113 void Foam::coupledLduMatrix::Amul
 114 (
 115     FieldField<Field, scalar>& result,
 116     const FieldField<Field, scalar>& x,
 117     const PtrList<FieldField<Field, scalar> >& bouCoeffs,
 118     const lduInterfaceFieldPtrsListList& interfaces,
 119     const direction cmpt
 120 ) const
 121 {
 122     const PtrList<lduMatrix>& matrices = *this;
 123
 124     // Reset product to zero
 125     result = 0;
 126
 127     // Initialise the update of coupled interfaces
 128     initMatrixInterfaces
 129     (
 130         bouCoeffs,
 131         interfaces,
 132         x,
 133         result,
 134         cmpt
 135     );
 136
 137     forAll (matrices, rowI)
 138     {
 139         matrices[rowI].AmulCore(result[rowI], x[rowI]);
 140     }
 141
 142     // Update couple interfaces
 143     updateMatrixInterfaces
 144     (
 145         bouCoeffs,
 146         interfaces,
 147         x,
 148         result,
 149         cmpt
 150     );
 151 }
 152
 153
 154 void Foam::coupledLduMatrix::Tmul
 155 (
 156     FieldField<Field, scalar>& result,
 157     const FieldField<Field, scalar>& x,
 158     const PtrList<FieldField<Field, scalar> >& intCoeffs,
 159     const lduInterfaceFieldPtrsListList& interfaces,
 160     const direction cmpt
 161 ) const
 162 {
 163     const PtrList<lduMatrix>& matrices = *this;
 164
 165     // Reset product to zero
 166     result = 0;
 167
 168     // Initialise the update of coupled interfaces
 169     initMatrixInterfaces
 170     (
 171         intCoeffs,
 172         interfaces,
 173         x,
 174         result,
 175         cmpt
 176     );
 177
 178     forAll (matrices, rowI)
 179     {
 180         matrices[rowI].TmulCore(result[rowI], x[rowI]);
 181     }
 182
 183     // Update couple interfaces
 184     updateMatrixInterfaces
 185     (
 186         intCoeffs,
 187         interfaces,
 188         x,
 189         result,
 190         cmpt
 191     );
 192 }
 193
 194
 195 void Foam::coupledLduMatrix::initMatrixInterfaces
 196 (
 197     const PtrList<FieldField<Field, scalar> >& coupleCoeffs,
 198     const lduInterfaceFieldPtrsListList& interfaces,
 199     const FieldField<Field, scalar>& x,
 200     FieldField<Field, scalar>& result,
 201     const direction cmpt
 202 ) const
 203 {
 204     const PtrList<lduMatrix>& matrices = *this;
 205
 206     forAll (matrices, rowI)
 207     {
 208         matrices[rowI].initMatrixInterfaces
 209         (
 210             coupleCoeffs[rowI],
 211             interfaces[rowI],
 212             x[rowI],
 213             result[rowI],
 214             cmpt
 215         );
 216     }
 217 }
 218
 219
 220 void Foam::coupledLduMatrix::updateMatrixInterfaces
 221 (
 222     const PtrList<FieldField<Field, scalar> >& coupleCoeffs,
 223     const lduInterfaceFieldPtrsListList& interfaces,
 224     const FieldField<Field, scalar>& x,
 225     FieldField<Field, scalar>& result,
 226     const direction cmpt
 227 ) const
 228 {
 229     const PtrList<lduMatrix>& matrices = *this;
 230
 231     forAll (matrices, rowI)
 232     {
 233         matrices[rowI].updateMatrixInterfaces
 234         (
 235             coupleCoeffs[rowI],
 236             interfaces[rowI],
 237             x[rowI],
 238             result[rowI],
 239             cmpt
 240         );
 241     }
 242 }
 243
 244
 245 // * * * * * * * * * * * * * * * Member Operators  * * * * * * * * * * * * * //
 246
 247
 248 // * * * * * * * * * * * * * * * Friend Functions  * * * * * * * * * * * * * //
 249
 250
 251 // * * * * * * * * * * * * * * * Friend Operators  * * * * * * * * * * * * * //
 252
 253
 254 // ************************************************************************* //