<div dir="ltr">Hello Christophe,<div><br></div><div>First, thanks for a great piece of software!</div><div><br></div><div>I am having troubles creating a 3D mesh that results in a single element type. I want a mesh that contains only tetrahedra, but somehow I end up with some pyramids that cause troubles downstream (e.g. the FENICS/DOLFIN library supports meshes with a single element type). The main problem is that, depending on the mesh size I set using the Mesh.MeshSizeFactor, the mesh ends up having only tets or mostly tets with a few pyramids.</div><div><br></div><div>The domain is a box with a few cylinders inside. Attached is a sketch of the geometry. The cylinders are created by extruding a circle, with no recombination. Attached is a sample GEO file, which I generate from pygmsh. Attached is also a sample geo file, but I am afraid is a bit hard to follow, because it is generated automatically by pygmsh, based on a python program I wrote to automate the mesh generation for different cases. The geo file also takes some time to run.</div><div><br></div><div>What algorithm type should I use (e.g. in the command line)?</div><div>What options should I set for Mesh.Algorithm and Mesh.Algorithm3D to make sure I get only tets?</div><div><br></div><div>Thanks,</div><div>David</div><div><br></div><div><br></div></div>