61 cylinderToFace::typeName,
62 "\n Usage: cylinderToFace (p1X p1Y p1Z) (p2X p2Y p2Z) radius\n\n"
63 " Select faces with centres 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)
114 if (innerRadius_ > radius_)
117 <<
"inner radius = " << innerRadius_ <<
" cannot be larger than "
118 <<
"outer radius = " << radius_
133 dict.getCompat<
point>(
"point1", {{
"p1", -2112}}),
145 const bool mandatoryInnerRadius
149 point1_(checkIs(is)),
150 point2_(checkIs(is)),
151 radius_(readScalar(checkIs(is))),
154 if (mandatoryInnerRadius)
156 innerRadius_ = readScalar(
checkIs(is));
183 Info<<
" Adding faces with centres within cylinder"
184 <<
", with point1 = " << point1_
185 <<
", point2 = " << point2_
186 <<
", radius = " << radius_;
188 if (innerRadius_ > 0)
190 Info<<
", inner radius = " << innerRadius_;
202 Info<<
" Removing faces with centres within cylinder"
203 <<
", with point1 = " << point1_
204 <<
", point2 = " << point2_
205 <<
", radius = " << radius_;
207 if (innerRadius_ > 0)
209 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 topoSetFaceSource to select all faces whose face centre inside a given bounding cylinder or cylinde...
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.
const vectorField & faceCentres() const
The topoSetFaceSource is a intermediate class for handling topoSet sources for selecting faces.
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.
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
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)
errorManipArg< error, int > exit(error &err, const int errNo=1)
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
dict add("bounds", meshBb)
#define forAll(list, i)
Loop across all elements in list.