Method moveFixedParticles

SRC/element/PFEMElement/BackgroundMesh.cpp:1514–1784

}

int BackgroundMesh::moveFixedParticles() {
    int ndm = OPS_GetNDM();

    // check each cell
    for (auto it = bcells.begin(); it != bcells.end(); ++it) {
        // get cell
        const VInt& index = it->first;
        BCell& cell = it->second;

        // empty cell
        if (cell.getPts().empty()) {
            continue;
        }

        // check if BACKGROUND_STRUCTURE cell
        if (cell.getType() != BACKGROUND_STRUCTURE) {
            continue;
        }

        // check neighbor cells
        VInt ind = index;
        VVInt indices;
        ind -= 1;
        getCorners(ind, 2, indices);

        // give each cell a score
        VInt scores(indices.size());
        for (int i = 0; i < (int)indices.size(); ++i) {
            auto cellit = bcells.find(indices[i]);
            if (cellit == bcells.end()) {
                // empty cell
                scores[i] = -1;
            } else if (cellit->second.getType() ==
                       BACKGROUND_STRUCTURE) {
                scores[i] = -1;
            } else if (cellit->second.getPts().empty()) {
                scores[i] = -1;
            } else {
                // easier to get into cell already having particles
                scores[i] = 1;
            }
        }
        VVInt cellmap;
        if (ndm == 2) {
            cellmap.resize(9);
            for (int i = 0; i < (int)cellmap.size(); ++i) {
                cellmap[i].resize(3);
                cellmap[i][0] = i;
            }
            cellmap[0][1] = 1;
            cellmap[0][2] = 3;
            cellmap[1][1] = 0;
            cellmap[1][2] = 2;
            cellmap[2][1] = 1;
            cellmap[2][2] = 5;
            cellmap[3][1] = 0;
            cellmap[3][2] = 6;
            cellmap[4][1] = 4;