40 const scalar coupledPolyPatch::defaultMatchTol_ = 1
e-4;
50 { transformType::UNKNOWN,
"unknown" },
51 { transformType::ROTATIONAL,
"rotational" },
52 { transformType::TRANSLATIONAL,
"translational" },
53 { transformType::COINCIDENTFULLMATCH,
"coincidentFullMatch" },
54 { transformType::NOORDERING,
"noOrdering" },
62 os <<
"v " << pt.x() <<
' ' << pt.y() <<
' ' << pt.z() <<
endl;
94 os <<
"l " << vertI-1 <<
' ' << vertI <<
nl;
111 for (
const face&
f : faces)
115 if (foamToObj.insert(
f[fp], vertI))
125 os <<
' ' << foamToObj[
f[fp]]+1;
127 os <<
' ' << foamToObj[
f[0]]+1 <<
nl;
134 const UList<face>& faces,
146 anchors[facei] =
points[faces[facei][0]];
154 const face&
f = faces[facei];
164 for (label fp2 = 0; fp2 <
f.
size(); ++fp2)
166 if (
f[fp1] ==
f[fp2])
191 anchors[facei] =
points[faces[facei][0]];
202 const UList<face>& faces,
212 const point& cc = faceCentres[facei];
214 const face&
f = faces[facei];
218 scalar maxLenSqr = -GREAT;
221 scalar maxCmpt = -GREAT;
226 maxLenSqr =
max(maxLenSqr,
magSqr(pt - cc));
249 scalar minDistSqr = GREAT;
255 if (distSqr < minDistSqr)
257 minDistSqr = distSqr;
262 if (anchorFp == -1 ||
Foam::sqrt(minDistSqr) > tol)
273 if (distSqr == minDistSqr && fp != anchorFp)
276 <<
"Cannot determine unique anchor point on face " 277 << UIndirectList<point>(
points,
f)
279 <<
"Both at index " << anchorFp <<
" and " << fp
280 <<
" the vertices have the same distance " 282 <<
" to the anchor " << anchor
283 <<
". Continuing but results might be wrong." 307 Pout<<
"coupledPolyPatch::calcTransformTensors : " <<
name() <<
endl 308 <<
" transform:" << transformTypeNames[
transform] <<
nl 309 <<
" (half)size:" << Cf.size() <<
nl 310 <<
" absTol:" << absTol <<
nl 311 <<
" smallDist min:" <<
min(smallDist) <<
nl 312 <<
" smallDist max:" <<
max(smallDist) <<
nl 313 <<
" sum(mag(nf & nr)):" <<
sum(
mag(nf & nr)) <<
endl;
334 scalar error = absTol*
Foam::sqrt(1.0*Cf.size());
347 && (
sum(
mag(nf & nr)) < Cf.size() - error)
357 forwardT_.setSize(Cf.size());
358 reverseT_.setSize(Cf.size());
359 collocated_.setSize(Cf.size());
370 Pout<<
" sum(mag(forwardT_ - forwardT_[0])):" 371 <<
sum(
mag(forwardT_ - forwardT_[0]))
375 if (
sum(
mag(forwardT_ - forwardT_[0])) < error)
377 forwardT_.setSize(1);
378 reverseT_.setSize(1);
379 collocated_.setSize(1);
383 Pout<<
" difference in rotation less than" 384 <<
" local tolerance " 385 << error <<
". Assuming uniform rotation." <<
endl;
396 separation_ = Cr - Cf;
398 collocated_.setSize(separation_.size());
406 bool sameSeparation =
true;
407 bool doneWarning =
false;
409 forAll(separation_, facei)
411 scalar smallSqr =
sqr(smallDist[facei]);
413 collocated_[facei] = (
magSqr(separation_[facei]) < smallSqr);
416 if (
magSqr(separation_[facei] - separation_[0]) > smallSqr)
418 sameSeparation =
false;
420 if (!doneWarning &&
debug)
424 Pout<<
" separation " << separation_[facei]
426 <<
" differs from separation[0] " << separation_[0]
427 <<
" by more than local tolerance " 429 <<
". Assuming non-uniform separation." <<
endl;
441 Pout<<
" separation " <<
mag(separation_[0])
442 <<
" less than local tolerance " << smallDist[0]
443 <<
". Assuming zero separation." <<
endl;
453 Pout<<
" separation " <<
mag(separation_[0])
454 <<
" more than local tolerance " << smallDist[0]
455 <<
". Assuming uniform separation." <<
endl;
458 separation_.setSize(1);
467 Pout<<
" separation_:" << separation_.size() <<
nl 468 <<
" forwardT size:" << forwardT_.size() <<
endl;
481 const polyBoundaryMesh& bm,
482 const word& patchType,
487 matchTolerance_(defaultMatchTol_),
498 const word& patchType
502 matchTolerance_(
dict.getOrDefault(
"matchTolerance", defaultMatchTol_)),
505 transformTypeNames.getOrDefault
522 matchTolerance_(
pp.matchTolerance_),
523 transform_(
pp.transform_)
534 matchTolerance_(
pp.matchTolerance_),
535 transform_(
pp.transform_)
549 matchTolerance_(
pp.matchTolerance_),
550 transform_(
pp.transform_)
564 matchTolerance_(
pp.matchTolerance_),
565 transform_(
pp.transform_)
583 os.
writeEntry(
"transform", transformTypeNames[transform_]);
static const Enum< transformType > transformTypeNames
void size(const label n)
Older name for setAddressableSize.
void cmptMax(FieldField< Field, typename FieldField< Field, Type >::cmptType > &cf, const FieldField< Field, Type > &f)
dimensioned< Type > sum(const DimensionedField< Type, GeoMesh > &f1)
virtual void write(Ostream &os) const
Write the polyPatch data as a dictionary.
A class for handling file names.
A face is a list of labels corresponding to mesh vertices.
coupledPolyPatch(const word &name, const label size, const label start, const label index, const polyBoundaryMesh &bm, const word &patchType, const transformType transform)
Construct from components.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
labelList pointLabels(nPoints, -1)
A list of keyword definitions, which are a keyword followed by a number of values (eg...
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
Output to file stream as an OSstream, normally using std::ofstream for the actual output...
constexpr char nl
The newline '\n' character (0x0a)
dimensionedScalar sqrt(const dimensionedScalar &ds)
Ostream & endl(Ostream &os)
Add newline and flush stream.
The coupledPolyPatch is an abstract base class for patches that couple regions of the computational d...
Ostream & writeEntry(const keyType &key, const T &value)
Write a keyword/value entry.
Smooth ATC in cells next to a set of patches supplied by type.
tensor rotationTensor(const vector &n1, const vector &n2)
Rotational transformation tensor from vector n1 to n2.
Various functions to operate on Lists.
#define forAll(list, i)
Loop across all elements in list.
static pointField getAnchorPoints(const UList< face > &, const pointField &, const transformType)
Get a unique anchor point for all faces.
vectorField pointField
pointField is a vectorField.
const dimensionedScalar e
Elementary charge.
word name(const expressions::valueTypeCode typeCode)
A word representation of a valueTypeCode. Empty for expressions::valueTypeCode::INVALID.
virtual void write(Ostream &os) const
Write the polyPatch data as a dictionary.
A class for handling words, derived from Foam::string.
Field< scalar > scalarField
Specialisation of Field<T> for scalar.
static void writeOBJ(Ostream &os, const point &pt)
Write point in OBJ format.
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
A polyBoundaryMesh is a polyPatch list with additional search methods and registered IO...
void cmptMag(FieldField< Field, Type > &cf, const FieldField< Field, Type > &f)
An Ostream is an abstract base class for all output systems (streams, files, token lists...
A Vector of values with scalar precision, where scalar is float/double depending on the compilation f...
int debug
Static debugging option.
OBJstream os(runTime.globalPath()/outputName)
defineTypeNameAndDebug(combustionModel, 0)
static label getRotation(const pointField &points, const face &f, const point &anchor, const scalar tol)
Get the number of vertices face f needs to be rotated such that.
static scalarField calcFaceTol(const UList< face > &faces, const pointField &points, const pointField &faceCentres)
Calculate typical tolerance per face. Is currently max distance.
void calcTransformTensors(const vectorField &Cf, const vectorField &Cr, const vectorField &nf, const vectorField &nr, const scalarField &smallDist, const scalar absTol, const transformType=UNKNOWN) const
Calculate the transformation tensors.
virtual ~coupledPolyPatch()
Destructor.
vector point
Point is a vector.
#define WarningInFunction
Report a warning using Foam::Warning.
Enum is a wrapper around a list of names/values that represent particular enumeration (or int) values...
Field< vector > vectorField
Specialisation of Field<T> for vector.
A patch is a list of labels that address the faces in the global face list.
dimensionSet transform(const dimensionSet &ds)
Return the argument; transformations do not change the dimensions.
List< bool > boolList
A List of bools.
prefixOSstream Pout
OSstream wrapped stdout (std::cout) with parallel prefix.
const volScalarField & p0
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
uindirectPrimitivePatch pp(UIndirectList< face >(mesh.faces(), faceLabels), mesh.points())