Method pointlist

SRC/element/PFEMElement/ParticleGroup.cpp:442–466

}

int ParticleGroup::pointlist(VDouble &pointdata) {
    int ndm = OPS_GetNDM();
    pointdata.clear();
    pointdata.reserve(particles.size() * (3 * ndm + 1));
    for (auto particle : particles) {
        auto tag = particle->getTag();
        const auto& crds = particle->getCrds(); 
        const auto& vel = particle->getVel(); 
        const auto& accel = particle->getAccel(); 
        double p = particle->getPressure(); 
        pointdata.push_back(tag);
        for (int j = 0; j < ndm; ++j) {
            pointdata.push_back(crds[j]);
        }
        for (int j = 0; j < ndm; ++j) {
            pointdata.push_back(vel[j]);
        }
        for (int j = 0; j < ndm; ++j) {
            pointdata.push_back(accel[j]);
        }
        pointdata.push_back(p);
    }

    return 0;
}

int ParticleGroup::pointlist(const VDouble &pointdata, int ndm) {
    VDouble crds(ndm);