Method moveParticles

SRC/element/PFEMElement/BackgroundMesh.cpp:3456–3529

}

int BackgroundMesh::moveParticles() {
    Domain* domain = OPS_GetDomain();
    if (domain == 0) return 0;
    int ndm = OPS_GetNDM();
    double dt = domain->getCurrentTime() - currentTime;

    // get current disp and velocity
    for (std::map<VInt, BNode>::iterator it = bnodes.begin();
         it != bnodes.end(); ++it) {
        BNode& bnode = it->second;
        VInt& tags = bnode.getTags();

        for (int i = 0; i < (int)bnode.size(); ++i) {
            Node* nd = domain->getNode(tags[i]);
            Pressure_Constraint* pc =
                domain->getPressure_Constraint(tags[i]);

            if (pc != 0) {
                auto& pn = bnode.getPressure();
                auto& dpn = bnode.getPdot();
                pn[i] = pc->getPressure();
                dpn[i] = pc->getPdot();
            }
            if (nd != 0) {
                const Vector& vel = nd->getTrialVel();
                const Vector& accel = nd->getTrialAccel();
                auto& vn = bnode.getVel();
                auto& dvn = bnode.getAccel();
                for (int j = 0; j < ndm; ++j) {
                    vn[i][j] = vel(j);
                    dvn[i][j] = accel(j);
                }
            }
        }
    }

    // store cells in a vector
    std::vector<BCell*> cells;
    VVInt indices;
    cells.reserve(bcells.size());
    indices.reserve(bcells.size());
    for (std::map<VInt, BCell>::iterator it = bcells.begin();
         it != bcells.end(); ++it) {
        indices.push_back(it->first);
        cells.push_back(&(it->second));
    }

    // move particles in each cell
    int res = 0;
#pragma omp parallel for
    for (int j = 0; j < (int)cells.size(); ++j) {
        // get particles in cell
        const VParticle& pts = cells[j]->getPts();

        // move the particle
        for (int i = 0; i < (int)pts.size(); ++i) {
            // set update state
            if (pts[i] == 0) continue;