30 {
32 and !
data.d_modelDeck_p->d_populateElementNodeConnectivity)
33 return;
34
35 for (
auto &p :
data.d_referenceParticles) {
36 auto &particle_mesh_p = p->getMeshP();
37 if (!particle_mesh_p->d_encDataPopulated && particle_mesh_p->d_enc.empty())
38 particle_mesh_p->readElementData(particle_mesh_p->d_filename);
39 }
40
41 size_t totalQuadPoints = 0;
42 const auto quadOrder =
data.d_modelDeck_p->d_quadOrder;
43 for (
auto &p :
data.d_particlesListTypeAll) {
44 const auto &particle_mesh_p = p->getMeshP();
45 const auto nq = particle_mesh_p->getNumElements() *
47
48 p->d_globQuadStart = totalQuadPoints;
49 totalQuadPoints += nq;
50 p->d_globQuadEnd = totalQuadPoints;
51
52 std::cout << std::format("p->id() = {}, "
53 "p->d_globQuadStart = {}, "
54 "totalQuadPoints = {}, "
55 "p->d_globQuadEnd = {}",
56 p->getId(), p->d_globQuadStart, nq, p->d_globQuadEnd)
57 << std::endl;
58 }
59
60 data.d_xQuadCur.resize(totalQuadPoints);
61 data.d_strain.resize(totalQuadPoints);
62 data.d_stress.resize(totalQuadPoints);
63}
size_t numQuadPoints(size_t elemType, size_t quadOrder)
bool isTagInList(const std::string &tag, const std::vector< std::string > &tags)
Returns true if tag is found in the list of tags.