cyclicACMIGAMGInterfaceField.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration |
5  \\ / A nd | Copyright (C) 2013 OpenFOAM Foundation
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 
28 #include "lduMatrix.H"
29 
30 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
31 
32 namespace Foam
33 {
34  defineTypeNameAndDebug(cyclicACMIGAMGInterfaceField, 0);
36  (
37  GAMGInterfaceField,
38  cyclicACMIGAMGInterfaceField,
39  lduInterface
40  );
42  (
43  GAMGInterfaceField,
44  cyclicACMIGAMGInterfaceField,
45  lduInterfaceField
46  );
47 }
48 
49 
50 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
51 
52 Foam::cyclicACMIGAMGInterfaceField::cyclicACMIGAMGInterfaceField
53 (
54  const GAMGInterface& GAMGCp,
55  const lduInterfaceField& fineInterface
56 )
57 :
58  GAMGInterfaceField(GAMGCp, fineInterface),
59  cyclicACMIInterface_(refCast<const cyclicACMIGAMGInterface>(GAMGCp)),
60  doTransform_(false),
61  rank_(0)
62 {
64  refCast<const cyclicAMILduInterfaceField>(fineInterface);
65 
66  doTransform_ = p.doTransform();
67  rank_ = p.rank();
68 }
69 
70 
71 Foam::cyclicACMIGAMGInterfaceField::cyclicACMIGAMGInterfaceField
72 (
73  const GAMGInterface& GAMGCp,
74  const bool doTransform,
75  const int rank
76 )
77 :
78  GAMGInterfaceField(GAMGCp, doTransform, rank),
79  cyclicACMIInterface_(refCast<const cyclicACMIGAMGInterface>(GAMGCp)),
80  doTransform_(doTransform),
81  rank_(rank)
82 {}
83 
84 
85 // * * * * * * * * * * * * * * * * Desstructor * * * * * * * * * * * * * * * //
86 
88 {}
89 
90 
91 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
92 
94 (
95  scalarField& result,
96  const scalarField& psiInternal,
97  const scalarField& coeffs,
98  const direction cmpt,
100 ) const
101 {
102  // Get neighbouring field
103  scalarField pnf
104  (
105  cyclicACMIInterface_.neighbPatch().interfaceInternalField(psiInternal)
106  );
107 
108  // Transform according to the transformation tensors
109  transformCoupleField(pnf, cmpt);
110 
111  if (cyclicACMIInterface_.owner())
112  {
113  pnf = cyclicACMIInterface_.AMI().interpolateToSource(pnf);
114  }
115  else
116  {
117  pnf = cyclicACMIInterface_.neighbPatch().AMI().interpolateToTarget(pnf);
118  }
119 
120  const labelUList& faceCells = cyclicACMIInterface_.faceCells();
121 
122  forAll(faceCells, elemI)
123  {
124  result[faceCells[elemI]] -= coeffs[elemI]*pnf[elemI];
125  }
126 }
127 
128 
129 // ************************************************************************* //
unsigned char direction
Definition: direction.H:43
virtual bool doTransform() const =0
Is the transform required.
Namespace for OpenFOAM.
Abstract base class for GAMG agglomerated interface fields.
volScalarField & p
Definition: createFields.H:51
virtual void updateInterfaceMatrix(scalarField &result, const scalarField &psiInternal, const scalarField &coeffs, const direction cmpt, const Pstream::commsTypes commsType) const
Update result field based on interface functionality.
commsTypes
Types of communications.
Definition: UPstream.H:64
#define forAll(list, i)
Definition: UList.H:421
virtual int rank() const =0
Return rank of component for transform.
Macros for easy insertion into run-time selection tables.
A 1D vector of objects of type <T>, where the size of the vector is known and can be used for subscri...
Definition: HashTable.H:60
Abstract base class for GAMG agglomerated interfaces.
Definition: GAMGInterface.H:51
Abstract base class for cyclic AMI coupled interfaces.
An abstract base class for implicitly-coupled interface fields e.g. processor and cyclic patch fields...
addToRunTimeSelectionTable(ensightPart, ensightPartCells, istream)
defineTypeNameAndDebug(combustionModel, 0)