Method gridFSInoDT

SRC/element/PFEMElement/BackgroundMesh.cpp:2057–2850

FSI without Delauney Triangulation

// FSI without Delauney Triangulation
int BackgroundMesh::gridFSInoDT() {
    Domain* domain = OPS_GetDomain();
    if (domain == 0) return 0;
    int ndm = OPS_GetNDM();

    // store cells in a vector
    std::vector<BCell*> cells;
    for (auto& item : bcells) {
        auto& cell = item.second;
        if (cell.getType() == BACKGROUND_STRUCTURE) {
            cells.push_back(&cell);
        }
    }

    // create elements in each cell
    int numele = 0;
    if (ndm == 2) {
        // numele = 2;
        numele = 4;
    } else if (ndm == 3) {
        // numele = 6;
        numele = 12;
    }
    VVInt elends(numele * cells.size());
    VInt gtags(numele * cells.size());
    VInt contact3Ddir(numele * cells.size(), -1);
    int ndtag = Mesh::nextNodeTag();
// #pragma omp parallel for
    for (int j = 0; j < (int)cells.size(); ++j) {
        // get indices
        auto& cindices = cells[j]->getIndices();

        // get bnodes
        auto& cbnodes = cells[j]->getNodes();

        // get types
        VInt types(cbnodes.size());
        for (int i = 0; i < (int)cbnodes.size(); ++i) {
            types[i] = cbnodes[i]->getType();
        }

        // get fluid and structural local indices
        VInt flocal, slocal, fixlocal, fslocal;
        for (int i = 0; i < (int)types.size(); ++i) {
            int type = types[i];
            if (type == BACKGROUND_STRUCTURE) {
                slocal.push_back(i);
            } else if (type == BACKGROUND_FLUID_STRUCTURE) {
                fslocal.push_back(i);
            } else if (type == BACKGROUND_FLUID) {
                flocal.push_back(i);
            } else {
                fixlocal.push_back(i);
            }
        }

        // FSI or SSI contact
        if ((ndm == 2 && slocal.size() + fslocal.size() == 4) ||