52 const scalar
c = n1 & n2;
54 const scalar magSqrN3 =
magSqr(n3);
64 + (1 -
c)*
sqr(n3)/magSqrN3;
73 const scalar magSqrN4 =
magSqr(n4);
74 return -
I + 2*
sqr(n4)/magSqrN4;
148 a.
y()*a.
x()*(1 -
c) - a.
z()*
s,
149 a.
x()*a.
z()*(1 -
c) + a.
y()*
s,
151 a.
x()*a.
y()*(1 -
c) + a.
z()*
s,
153 a.
y()*a.
z()*(1 -
c) - a.
x()*
s,
155 a.
x()*a.
z()*(1 -
c) - a.
y()*
s,
156 a.
y()*a.
z()*(1 -
c) + a.
x()*
s,
274 template<
class Type1,
class Type2>
312 scalar
cos = vec & e0;
313 scalar
sin = vec & e1;
Templated 3D SphericalTensor derived from VectorSpace adding construction from 1 component,...
Templated 3D symmetric tensor derived from VectorSpace adding construction from 6 components,...
Convenience class to handle the input of constant rotational speed. Reads an omega entry with default...
gmvFile<< "tracers "<< particles.size()<< nl;forAllConstIter(Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().position().x()<< " ";}gmvFile<< nl;forAllConstIter(Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().position().y()<< " ";}gmvFile<< nl;forAllConstIter(Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().position().z()<< " ";}gmvFile<< nl;forAll(lagrangianScalarNames, i){ word name=lagrangianScalarNames[i];IOField< scalar > s(IOobject(name, runTime.name(), cloud::prefix, mesh, IOobject::MUST_READ, IOobject::NO_WRITE))
const scalar piByTwo(0.5 *pi)
const dimensionedScalar c
Speed of light in a vacuum.
dimensionedSymmTensor symm(const dimensionedSymmTensor &dt)
tensor rotationTensor(const vector &n1, const vector &n2)
Rotational transformation tensor from unit vector n1 to n2.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
Tensor< scalar > tensor
Tensor of scalars.
tensor Ra(const vector &a, const scalar omega)
Rotational transformation tensor about axis a by omega radians.
tensor Rx(const scalar &omega)
Rotational transformation tensor about the x-axis by omega radians.
dimensionedScalar sin(const dimensionedScalar &ds)
SphericalTensor< Cmpt > sph(const DiagTensor< Cmpt > &dt)
Return the spherical part of a diagonal tensor.
static const Identity< scalar > I
tensor Ry(const scalar &omega)
Rotational transformation tensor about the y-axis by omega radians.
Type1 transformMask(const Type2 &t)
sphericalTensor transformMask< sphericalTensor >(const symmTensor &st)
dimensionSet transform(const dimensionSet &)
scalar pseudoAngle(const vector &e0, const vector &e1, const vector &vec)
Estimate angle of vec in coordinate system (e0, e1, e0^e1).
tensor Rz(const scalar &omega)
Rotational transformation tensor about the z-axis by omega radians.
dimensionSet perpendicular(const dimensionSet &)
symmTensor transformMask< symmTensor >(const symmTensor &st)
dimensioned< scalar > magSqr(const dimensioned< Type > &)
dimensionedScalar cos(const dimensionedScalar &ds)