Hi again, Ok. Seems sensible to me. Now, I'm a nuclear physicist so FEM is not quite what I do for a living, but presently I'm giving it a try to see if it can improve some detector simulations we are doing, presently done with finite difference solvers, so bare with me even if I'm not 100% on target all the time. I`ve already tried with dealII and that worked ok, but for good or bad reasons (I'm not sure;) I want simple tets with first order shape functions, hence libmesh.
As for getting inverted element, this means I've moved a node to the other side of the opposite face of my tet? This could very well be a bug in my code. I can, and will, check that I don't move nodes in an unreasnoble way. Also, I'll get back with a stack/call trace and the nodes of the failing element etc. But Paris time is 20:13 and I got other duties calling for me so I guess that will be tomorrows project. cheers Joa On Tue, Oct 13, 2009 at 12:47:26PM -0500, Roy Stogner wrote: > > On Tue, 13 Oct 2009, Joa Ljungvall wrote: > > >If I change the locations of the new nodes when they are made, i.e. somewhere > >deep down in the MeshRefinement class, can I do what I would like to do? > > I don't think that would make any difference. If moving nodes > eventually gives you an inverted element, it probably won't matter > exactly when the nodes are moved. > > I do think it would be a big mistake to start planning major low-level > changes before we fully understand why the less intrusive algorithm is > failing. Maybe get at that negative Jacobian failure in the debugger: > is the inverted element active? A parent? If the latter, what code > is calculating the Jacobian, and can it be worked around somehow? > --- > Roy ------------------------------------------------------------------------------ Come build with us! The BlackBerry(R) Developer Conference in SF, CA is the only developer event you need to attend this year. Jumpstart your developing skills, take BlackBerry mobile applications to market and stay ahead of the curve. Join us from November 9 - 12, 2009. Register now! http://p.sf.net/sfu/devconference _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
