Method reflect

tools/Wires/Inflator/IsotropicPeriodicInflator.cpp:185–217

}

void IsotropicPeriodicInflator::reflect() {
    const Float tol = get_distance_threshold();
    typedef std::function<bool(const VectorF&)> IndicatorFunc;
    std::vector<IndicatorFunc> min_bd_indicators;
    min_bd_indicators.push_back([=](const VectorF& v)
            { return fabs(v[0] - m_octa_cell_bbox_min[0]) < tol; });
    min_bd_indicators.push_back([=](const VectorF& v)
            { return fabs(v[1] - m_octa_cell_bbox_min[1]) < tol; });
    min_bd_indicators.push_back([=](const VectorF& v)
            { return fabs(v[2] - m_octa_cell_bbox_min[2]) < tol; });

    for (auto f : min_bd_indicators) {
        SubMesh::Ptr submesh = SubMesh::create_raw(m_vertices, m_faces);
        submesh->filter_vertex_with_custom_function(f);
        submesh->finalize();
        m_vertices = submesh->get_unselected_vertices();
        m_faces = submesh->get_unselected_faces();
    }
    
    for (size_t i=0; i<3; i++) {
        size_t num_vertices = m_vertices.rows();
        MatrixFr r_vts = m_vertices;
        for (size_t j=0; j<num_vertices; j++) {
            r_vts(j, i) = -r_vts(j, i) + 2 * m_octa_cell_bbox_min[i];
        }

        MatrixIr r_faces = m_faces.array() + num_vertices;
        r_faces.col(0).swap(r_faces.col(1));

        m_vertices = MatrixUtils::vstack<MatrixFr>({m_vertices, r_vts});
        m_faces = MatrixUtils::vstack<MatrixIr>({m_faces, r_faces});
    }
}

void IsotropicPeriodicInflator::merge_boundary() {
    auto bd = Boundary::extract_surface_boundary_raw(m_vertices, m_faces);