Method refine_bd_loops

tools/Wires/Inflator/PeriodicBoundaryRemesher.cpp:378–402

}

void PeriodicBoundaryRemesher::refine_bd_loops(Float ave_edge_len) {
    // Only refine bd edges lying on the edge of the bbox.
    // i.e. These bd edges and all their adjacent faces are all on the boundary.

    const size_t num_vertices = m_vertices.rows();
    assert(m_vertices.cols() == 3);
    std::vector<Float> refined_vertices(num_vertices * 3);
    std::copy(m_vertices.data(), m_vertices.data() + num_vertices * 3,
            refined_vertices.begin());

    refine_bd_loops(min_axis_marker[X], ave_edge_len, refined_vertices);
    refine_bd_loops(max_axis_marker[X], ave_edge_len, refined_vertices);
    refine_bd_loops(min_axis_marker[Y], ave_edge_len, refined_vertices);
    refine_bd_loops(max_axis_marker[Y], ave_edge_len, refined_vertices);
    refine_bd_loops(min_axis_marker[Z], ave_edge_len, refined_vertices);
    refine_bd_loops(max_axis_marker[Z], ave_edge_len, refined_vertices);

    assert(refined_vertices.size() % 3 == 0);
    const size_t num_refined_vertices = refined_vertices.size() / 3;
    if (num_vertices != num_refined_vertices) {
        m_vertices.resize(num_refined_vertices, 3);
        std::copy(refined_vertices.begin(), refined_vertices.end(),
                m_vertices.data());
    }
}

void PeriodicBoundaryRemesher::refine_bd_loops(short label, Float ave_edge_len,
        std::vector<Float>& refined_vertices) {