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" 69 void Foam::cylinderToFace::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 topoSetFaceSource(
mesh),
112 innerRadius_(innerRadius)
114 if (innerRadius_ > radius_)
117 <<
"inner radius = " << innerRadius_ <<
" cannot be larger than " 118 <<
"outer radius = " << radius_
126 const polyMesh&
mesh,
127 const dictionary&
dict 133 dict.getCompat<
point>(
"point1", {{
"p1", -2112}}),
143 const polyMesh&
mesh,
145 const bool mandatoryInnerRadius
148 topoSetFaceSource(
mesh),
149 point1_(checkIs(is)),
150 point2_(checkIs(is)),
151 radius_(readScalar(checkIs(is))),
154 if (mandatoryInnerRadius)
156 innerRadius_ = readScalar(
checkIs(is));
163 const polyMesh&
mesh,
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_;
virtual void applyToSet(const topoSetSource::setAction action, topoSet &set) const
Apply specified action to the topoSet.
errorManipArg< error, int > exit(error &err, const int errNo=1)
error FatalError
Error stream (stdout output on all processes), with additional 'FOAM FATAL ERROR' header text and sta...
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
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.
A topoSetFaceSource to select all faces whose face centre inside a given bounding cylinder or cylinde...
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)
#define forAll(list, i)
Loop across all elements in list.
cylinderToFace(const polyMesh &mesh, Istream &is, const bool mandatoryInnerRadius)
Construct from Istream with mandatory inner radius.
AccessType combine(const UList< T > &lists, AccessOp aop=accessOp< T >())
Combines sub-lists into a single list.
vectorField pointField
pointField is a vectorField.
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)
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)
const vectorField & faceCentres() const
Subtract elements from current set.
vector point
Point is a vector.
Class with constructor to add usage string to table.
messageStream Info
Information stream (stdout output on master, null elsewhere)
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)