Yes, sfepy cannot do hanging nodes, so i implemented the uniform refinement. The order p has to be the same in all elements of a field. It can be > 2, but then the quadratures and lagrangian basis limit its usefullness for now.
r.
----- Reply message ----- From: "Ondrej Certik" ond...@certik.cz To: sfepy...@googlegroups.com Subject: tetrahedral refinement Date: Fri, Apr 29, 2011 00:29
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
-- You received this message because you are subscribed to the Google Groups "sfepy-devel" group. To post to this group, send email to sfepy...@googlegroups.com. To unsubscribe from this group, send email to sfepy-devel...@googlegroups.com. For more options, visit this group at http://groups.google.com/group/sfepy-devel?hl=en.