Fix VTK output of polyhedrons

By now, the code only manages diamond cells as polyhedrons, but output
is general (up to the simple definition of a polyhedral cell type)