| 6 | using namespace PyMesh; |
| 7 | |
| 8 | void FaceCentroidAttribute::compute_from_mesh(Mesh& mesh) { |
| 9 | const size_t dim = mesh.get_dim(); |
| 10 | const size_t num_faces = mesh.get_num_faces(); |
| 11 | const size_t vertex_per_face = mesh.get_vertex_per_face(); |
| 12 | |
| 13 | VectorF& centroids = m_values; |
| 14 | centroids.resize(num_faces * dim); |
| 15 | |
| 16 | for (size_t i=0; i<num_faces; i++) { |
| 17 | VectorI face = mesh.get_face(i); |
| 18 | VectorF centroid = VectorF::Zero(dim); |
| 19 | for (size_t j=0; j<vertex_per_face; j++) { |
| 20 | centroid += mesh.get_vertex(face[j]); |
| 21 | } |
| 22 | centroid /= vertex_per_face; |
| 23 | |
| 24 | centroids.segment(i*dim, dim) = centroid; |
| 25 | } |
| 26 | } |
nothing calls this directly
no test coverage detected