Method addParticles

SRC/element/PFEMElement/BackgroundMesh.cpp:1187–1267

for each particle find lower grid check if out of range, if out, remove it get the containing cell set bnodes of the cell add the particle to the cell add bnodes to the cell

// add the particle to the cell
// add bnodes to the cell
int BackgroundMesh::addParticles() {
    // for all particles
    TaggedObjectIter& meshes = OPS_getAllMesh();
    Mesh* mesh = 0;
    while ((mesh = dynamic_cast<Mesh*>(meshes())) != 0) {
        ParticleGroup* group = dynamic_cast<ParticleGroup*>(mesh);
        if (group == 0) {
            continue;
        }

        // check group tag
        if (group->getTag() == contact_tag) {
            opserr << "WARNING: the particle group tag "
                   << contact_tag;
            opserr
                << " is reserved for internal use. Please select a "
                   "different one\n";
            return -1;
        }

        // remove particles
        VInt rm(group->numParticles(), 0);

        // for all particles
        for (int j = 0; j < group->numParticles(); j++) {
            // get particle
            Particle* p = group->getParticle(j);
            if (p == 0) continue;

            // get particle coordinates
            const VDouble& crds = p->getCrds();

            // get index
            VInt index;
            lowerIndex(crds, index);

            // if out of range
            for (int i = 0; i < (int)index.size(); ++i) {
                if (index[i] < lower[i] || index[i] >= upper[i]) {
                    rm[j] = 1;
                    break;
                }
            }
            if (rm[j] == 1) continue;

            // get bcell
            BCell& bcell = bcells[index];

            // add particles
            bcell.add(p);

            // if initial check structure cell
            if (bcell.getType() == BACKGROUND_STRUCTURE) {
                continue;
            }

            // add bnodes of the cell
            if (bcell.getNodes().empty()) {