Go to the documentation of this file.
61 cylinderToFace::typeName,
62 "\n Usage: cylinderToFace (p1X p1Y p1Z) (p2X p2Y p2Z) radius\n\n"
63 " Select all faces with face centre 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))
105 const scalar innerRadius
112 innerRadius_(innerRadius)
114 if (innerRadius_ > radius_)
117 <<
"inner radius = " << innerRadius_ <<
"cannot be larger than "
118 <<
"outer radius = " << radius_
148 point1_(checkIs(is)),
149 point2_(checkIs(is)),
150 radius_(readScalar(checkIs(is))),
167 Info<<
" Adding faces with centre within cylinder,"
168 <<
" with p1 = " << point1_ <<
", p2 = " << point2_
169 <<
", radius = " << radius_;
171 if (innerRadius_ > 0)
173 Info<<
", innerRadius = " << innerRadius_;
185 Info<<
" Removing faces with centre within cylinder,"
186 <<
" with p1 = " << point1_ <<
", p2 = " << point2_
187 <<
", radius = " << radius_;
189 if (innerRadius_ > 0)
191 Info<<
", innerRadius = " << innerRadius_;
addToRunTimeSelectionTable(decompositionMethod, kahipDecomp, dictionary)
vectorField pointField
pointField is a vectorField.
A topoSetFaceSource to select all faces whose face centre inside a given bounding cylinder or cylinde...
Class with constructor to add usage string to table.
Ostream & endl(Ostream &os)
Add newline and flush stream.
setAction
Enumeration defining the valid actions.
The topoSetFaceSource is a intermediate class for handling topoSet sources for selecting faces.
T get(const word &keyword, enum keyType::option matchOpt=keyType::REGEX) const
Create a new set and ADD elements to it.
cylinderToFace(const polyMesh &mesh, const point &point1, const point &point2, const scalar radius, const scalar innerRadius=0)
Construct from components.
Mesh consisting of general polyhedral cells.
#define forAll(list, i)
Loop across all elements in list.
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
An Istream is an abstract base class for all input systems (streams, files, token lists etc)....
messageStream Info
Information stream (uses stdout - output is on the master only)
AccessType combine(const UList< T > &lists, AccessOp aop=accessOp< T >())
Combines sub-lists into a single list.
virtual void applyToSet(const topoSetSource::setAction action, topoSet &set) const
Apply specified action to the topoSet.
addNamedToRunTimeSelectionTable(topoSetCellSource, badQualityToCell, word, badQuality)
General set of labels of mesh quantity (points, cells, faces).
T getCheckOrDefault(const word &keyword, const T &deflt, const Predicate &pred, enum keyType::option matchOpt=keyType::REGEX) const
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
void add(FieldField< Field1, typename typeOfSum< Type1, Type2 >::type > &f, const FieldField< Field1, Type1 > &f1, const FieldField< Field2, Type2 > &f2)
Macros for easy insertion into run-time selection tables.
Subtract elements from the set.
Vector< scalar > vector
A scalar version of the templated Vector.
errorManipArg< error, int > exit(error &err, const int errNo=1)
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
static MinMax< T > ge(const T &minVal)
A semi-infinite range from minVal to the type max.
const vectorField & faceCentres() const
T getCheck(const word &keyword, const Predicate &pred, enum keyType::option matchOpt=keyType::REGEX) const
void addOrDelete(topoSet &set, const label id, const bool add) const
Add or delete id from set. Add when 'add' is true.
const polyMesh & mesh_
Reference to the mesh.
defineTypeNameAndDebug(combustionModel, 0)