40 lagrangianFieldDecomposer::lagrangianFieldDecomposer
43 const polyMesh& procMesh,
46 const Cloud<indexedParticle>& lagrangianPositions,
47 const List<SLList<indexedParticle*>*>& cellParticles
51 positions_(procMesh, cloudName, false),
52 particleIndices_(lagrangianPositions.size())
56 forAll(cellProcAddressing, procCelli)
58 label celli = cellProcAddressing[procCelli];
60 if (cellParticles[celli])
62 SLList<indexedParticle*>& particlePtrs = *cellParticles[celli];
64 forAllIter(SLList<indexedParticle*>, particlePtrs, iter)
66 const indexedParticle& ppi = *iter();
67 particleIndices_[pi++] = ppi.index();
82 particleIndices_.setSize(pi);
84 IOPosition<passiveParticle>(positions_).write();