Method gridFluid

SRC/element/PFEMElement/BackgroundMesh.cpp:1786–2054

}

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

    // 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));
    }

    // create elements in each cell
    bool use_center_node = true;
    int numele = 0;
    int numelenodes = 0;
    if (ndm == 2) {
        if (use_center_node) {
            numele = 4;
        } else {
            numele = 2;
        }
        numelenodes = 3;
    } else if (ndm == 3) {
        if (use_center_node) {
            numele = 12;
        } else {
            numele = 6;
        }
        numelenodes = 4;
    }
    VVInt elends(numele * cells.size());
    VInt gtags(numele * cells.size());
    int ndtag = Mesh::nextNodeTag();
// #pragma omp parallel for
// TODO: setCenterNode will add nodes to domain. Can't be in parallel.
    for (int j = 0; j < (int)cells.size(); ++j) {
        // structural cell
        if (cells[j]->getType() == BACKGROUND_STRUCTURE) continue;

        // find the group of this mesh
        std::map<int, int> numpts;
        auto& pts = cells[j]->getPts();
        for (int i = 0; i < (int)pts.size(); ++i) {
            numpts[pts[i]->getGroupTag()] += 1;
        }
        int num = 0;
        int gtag = 0;
        for (std::map<int, int>::iterator it = numpts.begin();
             it != numpts.end(); ++it) {
            if (num < it->second) {
                num = it->second;
                gtag = it->first;
            }