61 cylinderToPoint::typeName,
62 "\n Usage: cylinderToPoint (p1X p1Y p1Z) (p2X p2Y p2Z) radius\n\n"
63 " Select points within bounding cylinder\n\n"
73 const vector axis = (point2_ - point1_);
74 const scalar magAxis2 =
magSqr(axis);
75 const scalar orad2 =
sqr(radius_);
76 const scalar irad2 = innerRadius_ > 0 ?
sqr(innerRadius_) : -1;
82 const vector d = ctrs[elemi] - point1_;
83 const scalar magD = d & axis;
85 if ((magD > 0) && (magD < magAxis2))
87 const scalar d2 = (d & d) -
sqr(magD)/magAxis2;
88 if ((d2 < orad2) && (d2 > irad2))
105 const scalar innerRadius
112 innerRadius_(innerRadius)
125 dict.getCompat<
point>(
"point1", {{
"p1", -2112}}),
137 const bool mandatoryInnerRadius
141 point1_(checkIs(is)),
142 point2_(checkIs(is)),
143 radius_(readScalar(checkIs(is))),
146 if (mandatoryInnerRadius)
148 innerRadius_ = readScalar(
checkIs(is));
175 Info<<
" Adding points within cylinder"
176 <<
", with point1 = " << point1_
177 <<
", point2 = " << point2_
178 <<
", radius = " << radius_;
180 if (innerRadius_ > 0)
182 Info<<
", inner radius = " << innerRadius_;
194 Info<<
" Removing points within cylinder"
195 <<
", with point1 = " << point1_
196 <<
", point2 = " << point2_
197 <<
", radius = " << radius_;
199 if (innerRadius_ > 0)
201 Info<<
", inner radius = " << innerRadius_;
Macros for easy insertion into run-time selection tables.
#define addNamedToRunTimeSelectionTable(baseType, thisType, argNames, lookupName)
Add to construction table with 'lookupName' as the key.
#define addToRunTimeSelectionTable(baseType, thisType, argNames)
Add to construction table with typeName as the key.
An Istream is an abstract base class for all input systems (streams, files, token lists etc)....
static MinMax< scalar > ge(const scalar &minVal)
A semi-infinite range from minVal to the type max.
A topoSetPointSource to select all points which are inside a given bounding cylinder or cylinder annu...
virtual void applyToSet(const topoSetSource::setAction action, topoSet &set) const
Apply specified action to the topoSet.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
T getCompat(const word &keyword, std::initializer_list< std::pair< const char *, int > > compat, enum keyType::option matchOpt=keyType::REGEX) const
T getCheckOrDefault(const word &keyword, const T &deflt, const Predicate &pred, enum keyType::option matchOpt=keyType::REGEX) const
Mesh consisting of general polyhedral cells.
virtual const pointField & points() const
Return raw points.
The topoSetPointSource is a intermediate class for handling topoSet sources for selecting points.
Class with constructor to add usage string to table.
Base class of a source for a topoSet.
void addOrDelete(topoSet &set, const label id, const bool add) const
Add or delete id from set. Add when 'add' is true.
setAction
Enumeration defining various actions.
@ SUBTRACT
Subtract elements from current set.
@ ADD
Add elements to current set.
@ NEW
Create a new set and ADD elements to it.
const polyMesh & mesh_
Reference to the mesh.
static Istream & checkIs(Istream &is)
Check state of stream.
General set of labels of mesh quantity (points, cells, faces).
A Vector of values with scalar precision, where scalar is float/double depending on the compilation f...
A class for handling words, derived from Foam::string.
#define defineTypeNameAndDebug(Type, DebugSwitch)
Define the typeName and debug information.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
messageStream Info
Information stream (stdout output on master, null elsewhere)
vectorField pointField
pointField is a vectorField.
Ostream & endl(Ostream &os)
Add newline and flush stream.
void add(FieldField< Field1, typename typeOfSum< Type1, Type2 >::type > &f, const FieldField< Field1, Type1 > &f1, const FieldField< Field2, Type2 > &f2)
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
dict add("bounds", meshBb)
#define forAll(list, i)
Loop across all elements in list.