35 Foam::lagrangianFieldDecomposer::lagrangianFieldDecomposer
48 particleIndices_(lagrangianPositions.
size())
52 labelList decodedProcFaceAddressing(faceProcAddressing.
size());
54 forAll(faceProcAddressing, i)
56 decodedProcFaceAddressing[i] =
mag(faceProcAddressing[i]) - 1;
59 forAll(cellProcAddressing, procCelli)
61 label celli = cellProcAddressing[procCelli];
63 if (cellParticles[celli])
65 SLList<indexedParticle*>& particlePtrs = *cellParticles[celli];
69 const indexedParticle& ppi = *iter();
70 particleIndices_[
pi++] = ppi.index();
72 const label mappedTetFace =
73 decodedProcFaceAddressing.find(ppi.tetFace());
75 if (mappedTetFace == -1)
78 <<
"Face lookup failure." <<
nl 90 ppi.procTetPt(procMesh, procCelli, mappedTetFace)
97 particleIndices_.resize(
pi);
99 IOPosition<Cloud<passiveParticle>>(positions_).
write();
void size(const label n)
Older name for setAddressableSize.
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...
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
Template class for non-intrusive linked lists.
constexpr char nl
The newline '\n' character (0x0a)
#define forAll(list, i)
Loop across all elements in list.
label size() const noexcept
The number of elements in table.
void write(vtk::formatter &fmt, const Type &val, const label n=1)
Component-wise write of a value (N times)
const word cloudName(propsDict.get< word >("cloud"))
A class for handling words, derived from Foam::string.
constexpr scalar pi(M_PI)
errorManip< error > abort(error &err)
A class representing the concept of 0 (zero) that can be used to avoid manipulating objects known to ...
Mesh consisting of general polyhedral cells.
forAllConstIters(mixture.phases(), phase)