Method gridNodes

SRC/element/PFEMElement/BackgroundMesh.cpp:1269–1512

}

int BackgroundMesh::gridNodes() {
    // get domain
    int ndm = OPS_GetNDM();
    Domain* domain = OPS_GetDomain();
    if (domain == 0) return 0;

    // vector of iterators
    std::vector<std::map<VInt, BNode>::iterator> iters;
    iters.reserve(bnodes.size());
    for (std::map<VInt, BNode>::iterator it = bnodes.begin();
         it != bnodes.end(); ++it) {
        iters.push_back(it);
    }

    // vector of new objects
    int ndtag = Mesh::nextNodeTag();
    std::vector<Node*> newnodes(iters.size(), 0),
        newpnodes(iters.size(), 0);
    std::vector<Pressure_Constraint*> newpcs(iters.size(), 0);

    int res = 0;

#pragma omp parallel for
    for (int j = 0; j < (int)iters.size(); ++j) {
        // get iterator
        std::map<VInt, BNode>::iterator it = iters[j];

        // get cell
        const VInt& index = it->first;
        BNode& bnode = it->second;
        if (bnode.getType() == BACKGROUND_FIXED) {
            continue;
        }

        // coordinates
        VDouble crds;
        getCrds(index, crds);

        // get particles
        VParticle pts;
        VInt minind = index;
        VInt maxind = index;
        minind -= numave;
        maxind += numave;
        gatherParticles(minind, maxind, pts);

        // get information
        double wt = 0.0, pre = 0.0, pdot = 0.0;
        VDouble vel(ndm), accel(ndm);
        for (int i = 0; i < (int)pts.size(); ++i) {
            // get particle
            if (pts[i] == 0) {
                continue;
            }

            // particle coordinates
            const VDouble& pcrds = pts[i]->getCrds();