Method init_barycentric_solvers

tools/MeshUtils/PointLocator.cpp:97–114

}

void PointLocator::init_barycentric_solvers() {
    const size_t dim = m_mesh->get_dim();

    const size_t num_elements = m_elements.size() / m_vertex_per_element;
    m_barycentric_solvers = MatrixF::Zero(num_elements*dim, dim);
    m_last_vertices = MatrixF::Zero(num_elements, dim);

    for (size_t i=0; i<num_elements; i++) {
        VectorI elem = m_elements.segment(i*m_vertex_per_element, m_vertex_per_element);
        VectorF last_v = m_mesh->get_vertex(elem[dim]);
        MatrixFr M(dim, dim);
        for (size_t j=0; j<dim; j++) {
            M.row(j) = m_mesh->get_vertex(elem[j]) - last_v;
        }
        m_barycentric_solvers.block(i*dim, 0, dim, dim) = M.transpose().inverse();
        m_last_vertices.row(i) = last_v;
    }
}

void PointLocator::init_hash_grid() {
    const size_t dim = m_mesh->get_dim();