61 sphereToCell::typeName,
62 "\n Usage: sphereToCell (centreX centreY centreZ) radius\n\n" 63 " Select all cells with cellCentre within bounding sphere\n\n" 69 void Foam::sphereToCell::combine(topoSet&
set,
const bool add)
const 73 const scalar orad2 =
sqr(radius_);
74 const scalar irad2 = innerRadius_ > 0 ?
sqr(innerRadius_) : -1;
80 const scalar d2 =
magSqr(ctrs[elemi] - origin_);
82 if ((d2 < orad2) && (d2 > irad2))
97 const scalar innerRadius
100 topoSetCellSource(
mesh),
103 innerRadius_(innerRadius)
116 dict.getCompat<
vector>(
"origin", {{
"centre", -1806}}),
130 origin_(checkIs(is)),
131 radius_(readScalar(checkIs(is))),
148 Info<<
" Adding cells within sphere," 149 <<
" origin = " << origin_ <<
", radius = " << radius_;
151 if (innerRadius_ > 0)
153 Info<<
", innerRadius = " << innerRadius_;
165 Info<<
" Removing cells within sphere," 166 <<
" origin = " << origin_ <<
", radius = " << radius_;
168 if (innerRadius_ > 0)
170 Info<<
", innerRadius = " << innerRadius_;
A topoSetCellSource to select cells based on cell centres inside a given bounding sphere...
A list of keyword definitions, which are a keyword followed by a number of values (eg...
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)...
T getCheck(const word &keyword, const Predicate &pred, enum keyType::option matchOpt=keyType::REGEX) const
Find and return a T with additional checking FatalIOError if not found, or if the number of tokens is...
Ostream & endl(Ostream &os)
Add newline and flush 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.
sphereToCell(const polyMesh &mesh, const point &origin, const scalar radius, const scalar innerRadius=0)
Construct from components.
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)
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).
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)
Mesh consisting of general polyhedral cells.
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
virtual void applyToSet(const topoSetSource::setAction action, topoSet &set) const
Apply specified action to the topoSet.
addToRunTimeSelectionTable(functionObject, pointHistory, dictionary)