Method node_map_t

SRC/domain/domain/DomainModalProperties.cpp:207–274

        std::map<int, size_t> pos;

        node_map_t(Domain* domain, int ndm, int ndf) {
            size_t n = static_cast<size_t>(domain->getNumNodes());
            if (n > 0) {
                nodes.resize(n);
                node_ids.resize(n);
                node_u_flags.resize(n);
                Node* nodePtr;
                NodeIter& theNodes = domain->getNodes();
                size_t counter = 0;
                while ((nodePtr = theNodes()) != 0) {
                    // node
                    nodes[counter] = nodePtr;
                    // map
                    pos.emplace(std::make_pair(nodePtr->getTag(), counter));
                    // id (original), needed to see if this dof is fixed
                    const ID& id_source = nodePtr->getDOF_GroupPtr()->getID();
                    // id (remapped)
                    auto& id = node_ids[counter];
                    int node_ndf = nodePtr->getNumberDOF();
                    id.resize(node_ndf);
                    int offset = static_cast<int>(counter)* ndf; // each node has ndf DOFs in this context
                    for (int j = 0; j < node_ndf; ++j) {
                        if (j < ndf) {
                            // exclude pressure dofs
                            // we can detect a pressure node (node_ndf=4) in 3d.. but not in 2d!
                            if (node_ndf == 4 && ndf == 6 && j == 3) {
                                id(j) = DOF_EXCLUDED;
                            }
                            else {
                                if (id_source(j) == -1)
                                    id(j) = DOF_FIXED;
                                else
                                    id(j) = offset;
                            }
                            ++offset;
                        }
                        else {
                            // exclude any dof >= ndf
                            id(j) = DOF_EXCLUDED;
                        }
                    }
                    // local dof flags
                    std::vector<int>& flags = node_u_flags[counter];
                    flags.resize(node_ndf);
                    for (int j = 0; j < node_ndf; ++j)
                        flags[static_cast<size_t>(j)] = id(j) == DOF_EXCLUDED ? -1 : j;
                    // go on
                    ++counter;
                }
            }

#ifdef DMP_VERBOSE
            opserr << "NODE RE-MAPPING:\n";
            opserr << "num eq: " << (int)(n * ndf) << "\n";
            for (size_t i = 0; i < n; ++i) {
                auto node = nodes[i];
                const auto& id = node_ids[i];
                const auto& uf = node_u_flags[i];