Method generate_joints

tools/Wires/Inflator/SimpleInflator.cpp:153–190

}

void SimpleInflator::generate_joints() {
    const size_t dim = m_wire_network->get_dim();
    const size_t num_vertices = m_wire_network->get_num_vertices();
    const size_t num_edges = m_wire_network->get_num_edges();
    const MatrixFr& vertices = m_wire_network->get_vertices();
    const MatrixIr& edges = m_wire_network->get_edges();
    const size_t loop_size = m_profile->size();

    std::vector<std::vector<const MatrixFr*> > loops(num_vertices);
    std::vector<std::vector<size_t> > source_edge_indices(num_vertices);
    for (size_t i=0; i<num_edges; i++) {
        const auto& end_loops = m_end_loops[i];
        const auto& edge = edges.row(i);
        loops[edge[0]].push_back(&end_loops.first);
        loops[edge[1]].push_back(&end_loops.second);
        source_edge_indices[edge[0]].push_back(i);
        source_edge_indices[edge[1]].push_back(i);
    }

    for (size_t i=0; i<num_vertices; i++) {
        const auto& incident_loops = loops[i];
        const auto& edge_ids = source_edge_indices[i];
        assert(incident_loops.size() == edge_ids.size());
        size_t valance = incident_loops.size();
        MatrixFr pts(valance * loop_size + 1, dim);
        VectorI  source_ids(pts.rows());
        pts.row(0) = vertices.row(i);
        source_ids[0] = -1;
        size_t count = 0;
        for (auto loop : incident_loops) {
            pts.block(count*loop_size+1, 0, loop_size, dim) = *loop;
            source_ids.segment(count*loop_size+1, loop_size) =
                VectorI::Ones(loop_size) * edge_ids[count];
            count++;
        }
        generate_joint(pts, source_ids, i);
    }
}

void SimpleInflator::connect_end_loops() {
    const size_t dim = m_wire_network->get_dim();