On Tue, Mar 8, 2011 at 10:26 AM, Robert Cimrman <cimr...@ntc.zcu.cz> wrote:

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.

Thanks for the paper. I studied it a bit --- and it seems to me, that the algorithm in there can also refine non-uniformly, without spoiling the element quality, right?

So the version that you implemented is just doing regular refinements, e.g. split each tetra into several (8?) tetras?

I then have another question --- if nonuniform tetras are used, let's say generated by tetgen, or by hand, can you still have abitrary order "p" basis functions on them in sfepy, or only uniform "p" in the whole mesh? Or maybe only second order elements, for now?

Ondrej