Go to the documentation of this file.
65 <<
"Span should have two positive and one zero entry. Now:"
68 else if (span[dir] < VSMALL)
86 <<
"Span should have two positive and one zero entry. Now:"
98 const scalar nearestDistSqr
111 info.rawPoint()[normalDir_] = origin_[normalDir_];
116 if (dir != normalDir_)
118 if (info.rawPoint()[dir] < origin_[dir])
120 info.rawPoint()[dir] = origin_[dir];
122 else if (info.rawPoint()[dir] > origin_[dir]+span_[dir])
124 info.rawPoint()[dir] = origin_[dir]+span_[dir];
131 if (
magSqr(info.rawPoint() - sample) > nearestDistSqr)
156 if (
mag(dir[normalDir_]) < VSMALL)
163 scalar t = (origin_[normalDir_]-start[normalDir_]) / dir[normalDir_];
172 info.rawPoint() = start+t*dir;
173 info.rawPoint()[normalDir_] = origin_[normalDir_];
178 if (dir != normalDir_)
180 if (info.rawPoint()[dir] < origin_[dir])
186 else if (info.rawPoint()[dir] > origin_[dir]+span_[dir])
200 treeBoundBox bb(origin_, origin_+span_);
201 bb.min()[normalDir_] -= 1
e-6;
202 bb.max()[normalDir_] += 1
e-6;
204 if (!bb.contains(info.hitPoint()))
207 <<
"bb:" << bb <<
endl
208 <<
"origin_:" << origin_ <<
endl
209 <<
"span_:" << span_ <<
endl
210 <<
"normalDir_:" << normalDir_ <<
endl
211 <<
"hitPoint:" << info.hitPoint()
222 Foam::searchablePlate::searchablePlate
232 normalDir_(calcNormal(span_))
235 <<
" origin:" << origin_
236 <<
" origin+span:" << origin_+span_
239 bounds() =
boundBox(origin_, origin_ + span_);
243 Foam::searchablePlate::searchablePlate
252 normalDir_(calcNormal(span_))
255 <<
" origin:" << origin_
256 <<
" origin+span:" << origin_+span_
259 bounds() =
boundBox(origin_, origin_ + span_);
267 if (regions_.empty())
270 regions_[0] =
"region0";
289 centres[0] = origin_ + 0.5*span_;
291 radiusSqr.setSize(1);
302 auto& pts = tpts.ref();
305 pts[2] = origin_ + span_;
307 if (span_.x() < span_.y() && span_.x() < span_.z())
309 pts[1] = origin_ +
point(0, span_.y(), 0);
310 pts[3] = origin_ +
point(0, 0, span_.z());
312 else if (span_.y() < span_.z())
314 pts[1] = origin_ +
point(span_.x(), 0, 0);
315 pts[3] = origin_ +
point(0, 0, span_.z());
319 pts[1] = origin_ +
point(span_.x(), 0, 0);
320 pts[3] = origin_ +
point(0, span_.y(), 0);
331 (origin_.x() + span_.x()) >= bb.
min().
x()
332 && origin_.x() <= bb.
max().
x()
333 && (origin_.y() + span_.y()) >= bb.
min().
y()
334 && origin_.y() <= bb.
max().
y()
335 && (origin_.z() + span_.z()) >= bb.
min().
z()
336 && origin_.z() <= bb.
max().
z()
341 void Foam::searchablePlate::findNearest
352 info[i] = findNearest(
samples[i], nearestDistSqr[i]);
357 void Foam::searchablePlate::findLine
368 info[i] = findLine(start[i],
end[i]);
380 findLine(start,
end, info);
392 findLine(start,
end, nearestInfo);
397 if (nearestInfo[pointi].hit())
400 info[pointi][0] = nearestInfo[pointi];
404 info[pointi].
clear();
427 normal.setSize(info.size());
431 normal[i][normalDir_] = 1.0;
443 <<
"Volume type not supported for plate."
addToRunTimeSelectionTable(decompositionMethod, kahipDecomp, dictionary)
Defines the attributes of an object for which implicit objectRegistry management is supported,...
const Cmpt & x() const
Access to the vector x component.
static const char *const componentNames[]
virtual void findLineAll(const pointField &start, const pointField &end, List< List< pointIndexHit >> &) const
Get all intersections in order from start to end.
A class for managing temporary objects.
static constexpr const zero Zero
Global zero (0)
Ostream & endl(Ostream &os)
Add newline and flush stream.
const point & max() const
Maximum describing the bounding box.
T get(const word &keyword, enum keyType::option matchOpt=keyType::REGEX) const
virtual bool overlaps(const boundBox &bb) const
Does any part of the surface overlap the supplied bound box?
const Cmpt & z() const
Access to the vector z component.
const point & min() const
Minimum describing the bounding box.
#define forAll(list, i)
Loop across all elements in list.
virtual tmp< pointField > points() const
Get the points that define the surface.
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
This class describes the interaction of (usually) a face and a point. It carries the info of a succes...
#define DebugInFunction
Report an information message using Foam::Info.
Base class of (analytical or triangulated) surface. Encapsulates all the search routines....
addNamedToRunTimeSelectionTable(topoSetCellSource, badQualityToCell, word, badQuality)
scalarField samples(nIntervals, Zero)
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
Macros for easy insertion into run-time selection tables.
constexpr auto end(C &c) -> decltype(c.end())
Return iterator to the end of the container c.
errorManip< error > abort(error &err)
Vector< scalar > vector
A scalar version of the templated Vector.
errorManipArg< error, int > exit(error &err, const int errNo=1)
const Cmpt & y() const
Access to the vector y component.
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
PointIndexHit< point > pointIndexHit
virtual tmp< pointField > coordinates() const
Get representative set of element coordinates.
virtual void findLineAny(const pointField &start, const pointField &end, List< pointIndexHit > &) const
Return any intersection on segment from start to end.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
const dimensionedScalar e
Elementary charge.
void clear()
Clear the list, i.e. set size to zero.
static tmp< T > New(Args &&... args)
Construct tmp of T with forwarding arguments.
A bounding box defined in terms of min/max extrema points.
virtual const wordList & regions() const
Names of regions.
vector point
Point is a vector.
virtual void getNormal(const List< pointIndexHit > &, vectorField &normal) const
From a set of points and indices get the normal.
virtual void boundingSpheres(pointField ¢res, scalarField &radiusSqr) const
Get bounding spheres (centre and radius squared), one per element.
void setSize(const label newSize)
Alias for resize(const label)
defineTypeNameAndDebug(combustionModel, 0)
virtual void getVolumeType(const pointField &, List< volumeType > &) const
Determine type (inside/outside/mixed) for point. unknown if.
static constexpr direction nComponents
Number of components in this vector space.
virtual void getRegion(const List< pointIndexHit > &, labelList ®ion) const
From a set of points and indices get the region.