61 if (!isType<cyclicAMIPointPatch>(
p))
64 <<
"patch " << this->
patch().
index() <<
" not cyclicAMI type. "
65 <<
"Patch type = " << p.type()
85 if (!isType<cyclicAMIPointPatch>(this->
patch()))
88 <<
"Field type does not correspond to patch type for patch "
90 <<
"Field type: " << typeName <<
endl
91 <<
"Patch type: " << this->
patch().type()
105 cyclicAMIPatch_(ptf.cyclicAMIPatch_),
116 return cyclicAMIPatch_.coupled();
127 if (cyclicAMIPatch_.cyclicAMIPatch().owner())
138 refCast<const GeometricField<Type, pointPatchField, pointMesh>>
140 this->internalField()
144 refCast<const cyclicAMIPointPatchField<Type>>
146 fld.boundaryField()[nbrPatch.
index()]
150 Field<Type> ptFld(this->patchInternalField(pField));
151 Field<Type> nbrPtFld(nbr.patchInternalField(pField));
156 const tensor& forwardT = this->forwardT()[0];
157 const tensor& reverseT = this->reverseT()[0];
167 Field<Type> nbrFcFld(nbrPpi().pointToFaceInterpolate(nbrPtFld));
170 if (cyclicAMIPatch_.cyclicAMIPatch().applyLowWeightCorrection())
172 Field<Type> fcFld(ppi().pointToFaceInterpolate(ptFld));
175 cyclicAMIPatch_.cyclicAMIPatch().interpolate
184 cyclicAMIPatch_.cyclicAMIPatch().interpolate(nbrFcFld);
188 this->addToInternalField
191 ppi().faceToPointInterpolate(nbrFcFld)()
197 Field<Type> fcFld(ppi().pointToFaceInterpolate(ptFld));
200 if (cyclicAMIPatch_.cyclicAMIPatch().applyLowWeightCorrection())
202 Field<Type> nbrFcFld(nbrPpi().pointToFaceInterpolate(nbrPtFld));
205 cyclicAMIPatch_.cyclicAMIPatch().neighbPatch().interpolate
214 cyclicAMIPatch_.cyclicAMIPatch().neighbPatch().interpolate
221 nbr.addToInternalField
224 nbrPpi().faceToPointInterpolate(fcFld)()
Info<< nl<< "Wrote faMesh in vtk format: "<< writer.output().name()<< nl;}{ vtk::lineWriter writer(aMesh.points(), aMesh.edges(), fileName(aMesh.mesh().time().globalPath()/"finiteArea-edges"));writer.writeGeometry();writer.beginCellData(4);writer.writeProcIDs();{ Field< scalar > fld(faMeshTools::flattenEdgeField(aMesh.magLe(), true))
Field with dimensions and associated with geometry type GeoMesh which is used to size the field and a...
Generic templated field type.
Generic GeometricField class.
commsTypes
Types of communications.
A Coupled boundary condition for pointField.
Cyclic AMI front and back plane patch field.
virtual void swapAddSeparated(const Pstream::commsTypes commsType, Field< Type > &) const
Complete swap of patch point values and add to local values.
virtual bool coupled() const
Cyclic AMI point patch - place holder only.
const cyclicAMIPointPatch & neighbPatch() const
Return neighbour point patch.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
virtual label index() const
Return the index of this patch in the pointBoundaryMesh.
Foam::pointPatchFieldMapper.
const pointPatch & patch() const
Return patch.
Basic pointPatch represents a set of points from the mesh.
virtual label index() const =0
Return the index of this patch in the pointBoundaryMesh.
#define FatalIOErrorInFunction(ios)
Report an error message using Foam::FatalIOError.
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
To & refCast(From &r)
Reference type cast template function.
dimensionSet transform(const dimensionSet &ds)
Return the argument; transformations do not change the dimensions.
Ostream & endl(Ostream &os)
Add newline and flush stream.
errorManipArg< error, int > exit(error &err, const int errNo=1)