Hi,
SfePy can now uniformly refine tetrahedral meshes using Domain.refine().
It's quite simple to try - look at 'tests/test_domain.py'. It uses an
algorithm from [1], that prevents spoiling element quality in subsequent
refinements.
r.
[1] Juergen Bey: Tetrahedral grid refinement, Computing 55 (1995), no. 4, 355--378