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 74 const scalar orad2 =
sqr(radius_);
75 const scalar irad2 = innerRadius_ > 0 ?
sqr(innerRadius_) : -1;
81 const scalar d2 =
magSqr(ctrs[elemi] - origin_);
83 if ((d2 < orad2) && (d2 > irad2))
98 const scalar innerRadius
101 topoSetCellSource(
mesh),
104 innerRadius_(innerRadius)
115 origin_(
dict.getCompat<
vector>(
"origin", {{
"centre", -1806}})),
131 origin_(checkIs(is)),
132 radius_(readScalar(checkIs(is))),
149 Info<<
" Adding cells within sphere," 150 <<
" origin = " << origin_ <<
", radius = " << radius_;
152 if (innerRadius_ > 0)
154 Info<<
", innerRadius = " << innerRadius_;
166 Info<<
" Removing cells within sphere," 167 <<
" origin = " << origin_ <<
", radius = " << radius_;
169 if (innerRadius_ > 0)
171 Info<<
", innerRadius = " << innerRadius_;
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.
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)
virtual void applyToSet(const topoSetSource::setAction action, topoSet &set) const
Apply specified action to the topoSet.
addToRunTimeSelectionTable(functionObject, pointHistory, dictionary)