101 return st.
ii()*st.
ii()*st.
ii();
107 inline SphericalTensor<Cmpt>
inv(
const SphericalTensor<Cmpt>& st)
110 if (
mag(st.ii()) < VSMALL)
113 <<
"SphericalTensor not invertible, determinant:" 114 <<
det(st) <<
" tensor:" << st <<
nl 125 inline Foam::scalar
magSqr(
const SphericalTensor<Cmpt>& st)
133 inline Cmpt
cmptMax(
const SphericalTensor<Cmpt>& st)
157 inline Cmpt
cmptAv(
const SphericalTensor<Cmpt>& st)
165 inline SphericalTensor<Cmpt>
lerp 167 const SphericalTensor<Cmpt>& a,
184 if (
mag(st.
ii()) < VSMALL)
188 <<
" is not divisible due to a zero element in SphericalTensor:" 189 <<
"SphericalTensor = " << st
207 <<
"SphericalTensor = " << st
208 <<
" is not divisible due to a zero value in Cmpt:" 230 operator&(
const SphericalTensor<Cmpt>& st,
const Vector<Cmpt>& v)
258 operator&&(
const SphericalTensor<Cmpt>& st1,
const SphericalTensor<Cmpt>& st2)
260 return 3*st1.ii()*st2.ii();
275 class outerProduct<Cmpt, SphericalTensor<Cmpt>>
279 typedef SphericalTensor<Cmpt>
type;
void cmptMax(FieldField< Field, typename FieldField< Field, Type >::cmptType > &cf, const FieldField< Field, Type > &f)
tmp< DimensionedField< typename DimensionedField< Type, GeoMesh >::cmptType, GeoMesh >> cmptAv(const DimensionedField< Type, GeoMesh > &f1)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
error FatalError
Error stream (stdout output on all processes), with additional 'FOAM FATAL ERROR' header text and sta...
Cmpt cmptSum(const SphericalTensor< Cmpt > &st)
Return the sum of components of a SphericalTensor.
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
dimensionedSphericalTensor inv(const dimensionedSphericalTensor &dt)
An Istream is an abstract base class for all input systems (streams, files, token lists etc)...
constexpr char nl
The newline '\n' character (0x0a)
dimensionSet operator &&(const dimensionSet &ds1, const dimensionSet &ds2)
dimensionedScalar operator/(const scalar s1, const dimensionedScalar &ds2)
const Cmpt & y() const noexcept
Access to the vector y component.
typeOfRank< typename pTraits< arg1 >::cmptType, direction(pTraits< arg1 >::rank)+direction(pTraits< arg2 >::rank) >::type type
dimensionedScalar det(const dimensionedSphericalTensor &dt)
SphericalTensor< Cmpt > sph(const DiagTensor< Cmpt > &dt)
Return the spherical part of a DiagTensor as a SphericalTensor.
scalar diagSqr() const
The L2-norm squared of the diagonal.
A templated (3 x 3) diagonal tensor of objects of <T>, effectively containing 1 element, derived from VectorSpace.
const dimensionedScalar b
Wien displacement law constant: default SI units: [m.K].
const Cmpt & ii() const noexcept
dimensionedScalar tr(const dimensionedSphericalTensor &dt)
errorManip< error > abort(error &err)
Templated 3D Vector derived from VectorSpace adding construction from 3 components, element access using x(), y() and z() member functions and the inner-product (dot-product) and cross-product operators.
const Cmpt & x() const noexcept
Access to the vector x component.
dimensioned< Type > lerp(const dimensioned< Type > &a, const dimensioned< Type > &b, const scalar t)
const Cmpt & z() const noexcept
Access to the vector z component.
tmp< GeometricField< Type, faPatchField, areaMesh > > operator &(const faMatrix< Type > &, const DimensionedField< Type, areaMesh > &)
A class representing the concept of 0 (zero) that can be used to avoid manipulating objects known to ...
void cmptMin(FieldField< Field, typename FieldField< Field, Type >::cmptType > &cf, const FieldField< Field, Type > &f)
gmvFile<< "tracers "<< particles.size()<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().x()<< " ";}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().y()<< " ";}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().z()<< " ";}gmvFile<< nl;forAll(lagrangianScalarNames, i){ word name=lagrangianScalarNames[i];IOField< scalar > s(IOobject(name, runTime.timeName(), cloud::prefix, mesh, IOobject::MUST_READ, IOobject::NO_WRITE))
typeOfRank< typename pTraits< arg1 >::cmptType, direction(pTraits< arg1 >::rank)+direction(pTraits< arg2 >::rank) - 2 >::type type
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
SphericalTensor()=default
Default construct.
static constexpr const zero Zero
Global zero (0)