cyclicAMIGAMGInterfaceField.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration | Website: https://openfoam.org
5  \\ / A nd | Copyright (C) 2011-2018 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(cyclicAMIGAMGInterfaceField, 0);
36  (
37  GAMGInterfaceField,
38  cyclicAMIGAMGInterfaceField,
39  lduInterface
40  );
42  (
43  GAMGInterfaceField,
44  cyclicAMIGAMGInterfaceField,
45  lduInterfaceField
46  );
47 }
48 
49 
50 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
51 
52 Foam::cyclicAMIGAMGInterfaceField::cyclicAMIGAMGInterfaceField
53 (
54  const GAMGInterface& GAMGCp,
55  const lduInterfaceField& fineInterface
56 )
57 :
58  GAMGInterfaceField(GAMGCp, fineInterface),
59  cyclicAMIInterface_(refCast<const cyclicAMIGAMGInterface>(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::cyclicAMIGAMGInterfaceField::cyclicAMIGAMGInterfaceField
72 (
73  const GAMGInterface& GAMGCp,
74  const bool doTransform,
75  const int rank
76 )
77 :
78  GAMGInterfaceField(GAMGCp, doTransform, rank),
79  cyclicAMIInterface_(refCast<const cyclicAMIGAMGInterface>(GAMGCp)),
80  doTransform_(doTransform),
81  rank_(rank)
82 {}
83 
84 
85 // * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
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  const cyclicAMIGAMGInterface& thisInterface = cyclicAMIInterface_;
103  const cyclicAMIGAMGInterface& neiInterface = thisInterface.neighbPatch();
104 
105  // Get neighbouring field
106  scalarField pnf(neiInterface.interfaceInternalField(psiInternal));
107 
108  // Transform according to the transformation tensors
109  transformCoupleField(pnf, cmpt);
110 
111  // Transform and interpolate
112  scalarField pf(size(), Zero);
113  if (thisInterface.owner())
114  {
115  forAll(thisInterface.AMIs(), i)
116  {
117  const scalar r =
118  pow
119  (
120  inv(thisInterface.AMITransforms()[i]).R()(cmpt, cmpt),
121  rank()
122  );
123 
124  pf += thisInterface.AMIs()[i].interpolateToSource(r*pnf);
125  }
126  }
127  else
128  {
129  forAll(neiInterface.AMIs(), i)
130  {
131  const scalar r =
132  pow
133  (
134  neiInterface.AMITransforms()[i].R()(cmpt, cmpt),
135  rank()
136  );
137 
138  pf += neiInterface.AMIs()[i].interpolateToTarget(r*pnf);
139  }
140  }
141 
142  // Multiply the field by coefficients and add into the result
143  const labelUList& faceCells = cyclicAMIInterface_.faceCells();
144  forAll(faceCells, elemI)
145  {
146  result[faceCells[elemI]] -= coeffs[elemI]*pf[elemI];
147  }
148 }
149 
150 
151 // ************************************************************************* //
tmp< Field< Type > > interfaceInternalField(const UList< Type > &internalData) const
Return the interface internal field of the given field.
#define forAll(list, i)
Loop across all elements in list.
Definition: UList.H:428
GAMG agglomerated cyclic AMI interface.
commsTypes
Types of communications.
Definition: UPstream.H:64
dimensionedSphericalTensor inv(const dimensionedSphericalTensor &dt)
uint8_t direction
Definition: direction.H:45
An abstract base class for implicitly-coupled interface fields e.g. processor and cyclic patch fields...
virtual const cyclicAMIGAMGInterface & neighbPatch() const
Return neighbour patch.
Macros for easy insertion into run-time selection tables.
Abstract base class for cyclic AMI coupled interfaces.
virtual bool doTransform() const =0
Is the transform required.
static const zero Zero
Definition: zero.H:97
virtual const PtrList< AMIPatchToPatchInterpolation > & AMIs() const
Return a reference to the AMI interpolators.
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:61
addToRunTimeSelectionTable(ensightPart, ensightPartCells, istream)
defineTypeNameAndDebug(combustionModel, 0)
virtual int rank() const =0
Return rank of component for transform.
Abstract base class for GAMG agglomerated interfaces.
Definition: GAMGInterface.H:51
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
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.
#define R(A, B, C, D, E, F, K, M)
Abstract base class for GAMG agglomerated interface fields.
volScalarField & p
Namespace for OpenFOAM.
virtual bool owner() const
Does this side own the interface?
virtual const List< vectorTensorTransform > & AMITransforms() const