43 const twoDPointCorrector& correct2D
46 const vector& normal = correct2D.planeNormal();
74 typeIOobject<IOdictionary> motionObj
77 mesh.time().constant(),
79 IOobject::MUST_READ_IF_MODIFIED,
83 if (motionObj.headerOk())
85 Info<<
"Reading " << mesh.time().constant() /
"motionProperties"
88 IOdictionary motionProperties(motionObj);
90 Switch twoDMotion(motionProperties.lookup(
"twoDMotion"));
103 const polyMesh& mesh,
104 const twoDPointCorrector* correct2DPtr
111 ? getNormalDir(*correct2DPtr)
126 scalar maxX = -great;
130 scalar maxY = -great;
134 scalar maxZ = -great;
137 scalar minOther = great;
138 scalar maxOther = -great;
144 const edge&
e = edges[edgeI];
148 scalar eMag =
mag(eVec);
154 minX =
min(minX, eMag);
155 maxX =
max(maxX, eMag);
160 minY =
min(minY, eMag);
161 maxY =
max(maxY, eMag);
166 minZ =
min(minZ, eMag);
167 maxZ =
max(maxZ, eMag);
172 minOther =
min(minOther, eMag);
173 maxOther =
max(maxOther, eMag);
177 os <<
"Mesh bounding box:" << boundBox(mesh_.
points()) <<
nl <<
nl
178 <<
"Mesh edge statistics:" <<
nl
179 <<
" x aligned : number:" << nX <<
"\tminLen:" << minX
180 <<
"\tmaxLen:" << maxX <<
nl
181 <<
" y aligned : number:" << nY <<
"\tminLen:" << minY
182 <<
"\tmaxLen:" << maxY <<
nl
183 <<
" z aligned : number:" << nZ <<
"\tminLen:" << minZ
184 <<
"\tmaxLen:" << maxZ <<
nl
185 <<
" other : number:" << mesh_.
nEdges() - nX - nY - nZ
186 <<
"\tminLen:" << minOther
187 <<
"\tmaxLen:" << maxOther <<
nl <<
endl;
191 return min(minY,
min(minZ, minOther));
193 else if (normalDir_ == 1)
195 return min(minX,
min(minZ, minOther));
197 else if (normalDir_ == 2)
199 return min(minX,
min(minY, minOther));
203 return min(minX,
min(minY,
min(minZ, minOther)));
#define forAll(list, i)
Loop across all elements in list.
static const scalar edgeTol_
edgeStats(const polyMesh &mesh)
Construct from mesh.
scalar minLen(Ostream &os) const
Calculate minimum edge length and print.
virtual const pointField & points() const
Return raw points.
const edgeList & edges() const
Return mesh edges. Uses calcEdges.
autoPtr< CompressibleMomentumTransportModel > New(const volScalarField &rho, const volVectorField &U, const surfaceScalarField &phi, const viscosity &viscosity)
const dimensionedScalar e
Elementary charge.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Ostream & endl(Ostream &os)
Add newline and flush stream.
layerAndWeight min(const layerAndWeight &a, const layerAndWeight &b)
Vector< scalar > vector
A scalar version of the templated Vector.
dimensioned< scalar > mag(const dimensioned< Type > &)
layerAndWeight max(const layerAndWeight &a, const layerAndWeight &b)