Bug in mesh.num_cells()

[Johann Heller]

I think there is a bug in the mesh.num_cells() function, just try:

mesh = UnitCubeMesh(1,1,1)
mesh.num_cells()

this gives 6, the right answer would be 12

this only occures for (1,1,1), for all other it seems correct.

[Prof Garth Wells]

Why should it be 12?

A cube can be split into 6 tetrahedra.

[Johann Heller]

If you do

plot(mesh)

and hit "c" to show the cell numbering, you can see that the function makes 12 cells

or am i missing something ?

[Prof Garth Wells]

There is no bug in Mesh::num_cells(), but 'c' from the keyboard appears to annotate facet indices, not cells indices, which is inconsistent with the documentation.

I'll register a separate Issue.

[Prof Garth Wells]

Removing milestone: 1.5 (automated comment)