29 inline const T Foam::SVD::sign(
const T& a,
const T& b)
31 return b >= 0 ? (a >= 0 ? a : -a) : (a >= 0 ? -a : a);
68 if (s > VSMALL && s < minS) minS =
s;
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
void size(const label)
Override size to be inconsistent with allocated storage.
const scalarRectangularMatrix & VSinvUt() const
Return VSinvUt (the pseudo inverse)
const scalarRectangularMatrix & V() const
Return the square matrix V.
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.timeName(), cloud::prefix, mesh, IOobject::MUST_READ, IOobject::NO_WRITE))
const scalarRectangularMatrix & U() const
Return U.
scalar minNonZeroS() const
Return the minimum non-zero singular value.
const scalarDiagonalMatrix & S() const
Return the singular values.
label nZeros() const
Return the number of zero singular values.