61 cylinderToPoint::typeName,
62 "\n Usage: cylinderToPoint (p1X p1Y p1Z) (p2X p2Y p2Z) radius\n\n" 63 " Select points within bounding cylinder\n\n" 69 void Foam::cylinderToPoint::combine(topoSet&
set,
const bool add)
const 74 const vector axis = (point2_ - point1_);
75 const scalar magAxis2 =
magSqr(axis);
76 const scalar orad2 =
sqr(radius_);
77 const scalar irad2 = innerRadius_ > 0 ?
sqr(innerRadius_) : -1;
83 const vector d = ctrs[elemi] - point1_;
84 const scalar magD = d & axis;
86 if ((magD > 0) && (magD < magAxis2))
88 const scalar d2 = (d & d) -
sqr(magD)/magAxis2;
89 if ((d2 < orad2) && (d2 > irad2))
102 const polyMesh&
mesh,
106 const scalar innerRadius
109 topoSetPointSource(
mesh),
113 innerRadius_(innerRadius)
124 point1_(
dict.getCompat<
point>(
"point1", {{
"p1", -2112}})),
126 radius_(
dict.
getCompat<scalar>(
"radius", {{
"outerRadius", -2112}})),
138 const bool mandatoryInnerRadius
142 point1_(checkIs(is)),
143 point2_(checkIs(is)),
144 radius_(readScalar(checkIs(is))),
147 if (mandatoryInnerRadius)
149 innerRadius_ = readScalar(
checkIs(is));
156 const polyMesh&
mesh,
176 Info<<
" Adding points within cylinder" 177 <<
", with point1 = " << point1_
178 <<
", point2 = " << point2_
179 <<
", radius = " << radius_;
181 if (innerRadius_ > 0)
183 Info<<
", inner radius = " << innerRadius_;
195 Info<<
" Removing points within cylinder" 196 <<
", with point1 = " << point1_
197 <<
", point2 = " << point2_
198 <<
", radius = " << radius_;
200 if (innerRadius_ > 0)
202 Info<<
", inner radius = " << innerRadius_;
The topoSetPointSource is a intermediate class for handling topoSet sources for selecting points...
A list of keyword definitions, which are a keyword followed by a number of values (eg...
cylinderToPoint(const polyMesh &mesh, Istream &is, const bool mandatoryInnerRadius)
Construct from Istream with mandatory inner radius.
Create a new set and ADD elements to it.
Add elements to current set.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
An Istream is an abstract base class for all input systems (streams, files, token lists etc)...
Ostream & endl(Ostream &os)
Add newline and flush stream.
static Istream & checkIs(Istream &is)
Check state of stream.
void addOrDelete(topoSet &set, const label id, const bool add) const
Add or delete id from set. Add when 'add' is true.
Macros for easy insertion into run-time selection tables.
addNamedToRunTimeSelectionTable(topoSetCellSource, badQualityToCell, word, badQuality)
virtual const pointField & points() const
Return raw points.
#define forAll(list, i)
Loop across all elements in list.
AccessType combine(const UList< T > &lists, AccessOp aop=accessOp< T >())
Combines sub-lists into a single list.
vectorField pointField
pointField is a vectorField.
A topoSetPointSource to select all points which are inside a given bounding cylinder or cylinder annu...
static MinMax< scalar > ge(const scalar &minVal)
A semi-infinite range from minVal to the type max.
T getCheckOrDefault(const word &keyword, const T &deflt, const Predicate &pred, enum keyType::option matchOpt=keyType::REGEX) const
Find and return a T, or return the given default value. FatalIOError if it is found and the number of...
setAction
Enumeration defining various actions.
const polyMesh & mesh_
Reference to the mesh.
void add(FieldField< Field1, typename typeOfSum< Type1, Type2 >::type > &f, const FieldField< Field1, Type1 > &f1, const FieldField< Field2, Type2 > &f2)
A Vector of values with scalar precision, where scalar is float/double depending on the compilation f...
defineTypeNameAndDebug(combustionModel, 0)
General set of labels of mesh quantity (points, cells, faces).
T getCompat(const word &keyword, std::initializer_list< std::pair< const char *, int >> compat, enum keyType::option matchOpt=keyType::REGEX) const
Find and return a T using any compatibility names if needed. FatalIOError if not found, or if there are excess tokens.
Subtract elements from current set.
vector point
Point is a vector.
Class with constructor to add usage string to table.
virtual void applyToSet(const topoSetSource::setAction action, topoSet &set) const
Apply specified action to the topoSet.
tmp< pointField > transform(const pointField &points) const
Coordinate transform (optionally) coordinates. Returns reference to input data if no transform is act...
messageStream Info
Information stream (stdout output on master, null elsewhere)
Mesh consisting of general polyhedral cells.
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
addToRunTimeSelectionTable(functionObject, pointHistory, dictionary)