102 return st.
ii()*st.
ii()*st.
ii();
108 inline SphericalTensor<Cmpt>
inv(
const SphericalTensor<Cmpt>& st)
111 if (
mag(st.ii()) < VSMALL)
114 <<
"SphericalTensor is not invertible due to the zero determinant:" 115 <<
"det(SphericalTensor) = " <<
det(st)
126 inline Cmpt
magSqr(
const SphericalTensor<Cmpt>& st)
128 return Cmpt(3*
mag(st.ii()*st.ii()));
134 inline Cmpt
cmptMax(
const SphericalTensor<Cmpt>& st)
158 inline Cmpt
cmptAv(
const SphericalTensor<Cmpt>& st)
168 inline SphericalTensor<Cmpt>
172 if (
mag(st.
ii()) < VSMALL)
176 <<
" is not divisible due to a zero element in SphericalTensor:" 177 <<
"SphericalTensor = " << st
188 inline SphericalTensor<Cmpt>
189 operator/(
const SphericalTensor<Cmpt>& st,
const Cmpt
s)
195 <<
"SphericalTensor = " << st
196 <<
" is not divisible due to a zero value in Cmpt:" 202 return SphericalTensor<Cmpt>(st.ii()/
s);
218 operator&(
const SphericalTensor<Cmpt>& st,
const Vector<Cmpt>& v)
232 operator&(
const Vector<Cmpt>& v,
const SphericalTensor<Cmpt>& st)
246 operator&&(
const SphericalTensor<Cmpt>& st1,
const SphericalTensor<Cmpt>& st2)
248 return 3*st1.ii()*st2.ii();
263 class outerProduct<Cmpt, SphericalTensor<Cmpt>>
267 typedef SphericalTensor<Cmpt>
type;
void cmptMax(FieldField< Field, typename FieldField< Field, Type >::cmptType > &cf, const FieldField< Field, Type > &f)
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)...
dimensionSet operator &&(const dimensionSet &ds1, const dimensionSet &ds2)
dimensionedScalar operator/(const scalar s1, const dimensionedScalar &ds2)
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.
A templated (3 x 3) diagonal tensor of objects of <T>, effectively containing 1 element, derived from VectorSpace.
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.
tmp< GeometricField< Type, faPatchField, areaMesh > > operator &(const faMatrix< Type > &, const DimensionedField< Type, areaMesh > &)
const SphericalTensor< Cmpt > & T() const
Return non-Hermitian transpose (no-op)
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
tmp< DimensionedField< typename DimensionedField< Type, GeoMesh >::cmptType, GeoMesh >> cmptAv(const DimensionedField< Type, GeoMesh > &df)
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
SphericalTensor()=default
Default construct.
static constexpr const zero Zero
Global zero (0)