Method init_boundary_normals

tools/Assembler/Mesh/TriangleMesh.cpp:267–288

}

void TriangleMesh::init_boundary_normals(const BoundaryPtr& bd) {
    const size_t dim = m_mesh->get_dim();
    const size_t num_edges = m_boundary_edges.rows();
    m_boundary_normals.resize(num_edges, 3);
    for (size_t i=0; i<num_edges; i++) {
        const VectorI& edge = m_boundary_edges.row(i);
        assert(2 == edge.size());
        VectorF v0 = m_mesh->get_vertex(edge[0]);
        VectorF v1 = m_mesh->get_vertex(edge[1]);

        size_t face_idx = bd->get_boundary_element(i);
        VectorF normal = m_mesh->get_attribute("face_normal").segment(
                face_idx * 3, 3);

        Vector3F edge_v, normal_v;
        edge_v.segment(0, dim) = v1 - v0;
        normal_v.segment(0, 3) = normal;

        Vector3F edge_normal = edge_v.cross(normal_v).normalized();
        m_boundary_normals.row(i) = edge_normal;
    }
}

void TriangleMesh::init_boundary_vertex_normals() {
    const size_t dim = getDim();