Method gridFSI

SRC/element/PFEMElement/BackgroundMesh.cpp:2852–3290

}

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

    TriangleMeshGenerator gen;
    TetMeshGenerator tetgen;

    // gather bnodes
    std::map<VInt, BNode*> fsibnodes;
    for (std::map<VInt, BCell>::iterator it = bcells.begin();
         it != bcells.end(); ++it) {
        // only for structural cells
        BCell& bcell = it->second;
        if (bcell.getType() == BACKGROUND_FLUID) continue;

        // get bnode
        auto& bnodes = bcell.getNodes();
        auto& bindices = bcell.getIndices();
        for (int j = 0; j < (int)bnodes.size(); ++j) {
            BNode* bnode = bnodes[j];
            VInt bindex = bindices[j];
            if (bnode == 0) {
                opserr << "WARNING: failed to get bnode -- gridFSI\n";
                return -1;
            }
            fsibnodes[bindex] = bnode;
        }
    }

    // add points
    VInt ndtags, ndtypes, ndsids;
    VVInt ndindex;
    VDouble min, max;
    for (std::map<VInt, BNode*>::iterator it = fsibnodes.begin();
         it != fsibnodes.end(); ++it) {
        VInt bindex = it->first;
        BNode* bnode = it->second;

        const VInt& tags = bnode->getTags();
        const VVDouble& crdsn = bnode->getCrds();
        int type = bnode->getType();
        const VInt& sid = bnode->getSid();

        for (int i = -1; i < (int)tags.size(); ++i) {
            // -1 is only for fixed bnode
            if (i == -1 && type != BACKGROUND_FIXED) {
                continue;
            }
            if (i == -1) {
                ndtags.push_back(0);
                ndsids.push_back(0);
            } else {
                ndtags.push_back(tags[i]);
                ndsids.push_back(sid[i]);
            }
            ndtypes.push_back(type);
            ndindex.push_back(bindex);