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 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))
101 const polyMesh&
mesh,
105 const scalar innerRadius
108 topoSetPointSource(
mesh),
112 innerRadius_(innerRadius)
125 dict.getCompat<
point>(
"point1", {{
"p1", -2112}}),
135 const polyMesh&
mesh,
137 const bool mandatoryInnerRadius
140 topoSetPointSource(
mesh),
141 point1_(checkIs(is)),
142 point2_(checkIs(is)),
143 radius_(readScalar(checkIs(is))),
146 if (mandatoryInnerRadius)
148 innerRadius_ = readScalar(
checkIs(is));
155 const polyMesh&
mesh,
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_;
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)
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.
addToRunTimeSelectionTable(decompositionMethod, kahipDecomp, dictionary)
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.
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)