61 cylinderToCell::typeName,
62 "\n Usage: cylinderToCell (p1X p1Y p1Z) (p2X p2Y p2Z) radius\n\n" 63 " Select all cells with cell centre within bounding cylinder\n\n" 69 void Foam::cylinderToCell::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 topoSetCellSource(
mesh),
113 innerRadius_(innerRadius)
115 if (innerRadius_ > radius_)
118 <<
"inner radius = " << innerRadius_ <<
" cannot be larger than " 119 <<
"outer radius = " << radius_
127 const polyMesh&
mesh,
128 const dictionary&
dict 132 point1_(
dict.getCompat<
point>(
"point1", {{
"p1", -2112}})),
134 radius_(
dict.
getCompat<scalar>(
"radius", {{
"outerRadius", -2112}})),
146 const bool mandatoryInnerRadius
150 point1_(checkIs(is)),
151 point2_(checkIs(is)),
152 radius_(readScalar(checkIs(is))),
155 if (mandatoryInnerRadius)
157 innerRadius_ = readScalar(
checkIs(is));
164 const polyMesh&
mesh,
184 Info<<
" Adding cells with centres within cylinder" 185 <<
", with point1 = " << point1_
186 <<
", point2 = " << point2_
187 <<
", radius = " << radius_;
189 if (innerRadius_ > 0)
191 Info<<
", inner radius = " << innerRadius_;
203 Info<<
" Removing cells with centres within cylinder" 204 <<
", with point1 = " << point1_
205 <<
", point2 = " << point2_
206 <<
", radius = " << radius_;
208 if (innerRadius_ > 0)
210 Info<<
", inner radius = " << innerRadius_;
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.
cylinderToCell(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)...
A topoSetCellSource to select all cells whose cell centre inside a given bounding cylinder or cylinde...
Ostream & endl(Ostream &os)
Add newline and flush stream.
static Istream & checkIs(Istream &is)
Check state of stream.
The topoSetCellSource is a intermediate class for handling topoSet sources for selecting cells...
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.
AccessType combine(const UList< T > &lists, AccessOp aop=accessOp< T >())
Combines sub-lists into a single list.
vectorField pointField
pointField is a vectorField.
virtual void applyToSet(const topoSetSource::setAction action, topoSet &set) const
Apply specified action to the topoSet.
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 vectorField & cellCentres() const
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.
Subtract elements from current set.
vector point
Point is a vector.
Class with constructor to add usage string to table.
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)