Hi again, When I loop over the nodes, find a face that is on the boundary, and then move the nodes the problem is, if I understood it correct, that I might move a hanging node out of the line between the nodes of the not so refined neighbour. But I can check if this is the case, and instead of moving only the hanging node I move all three nodes keeping them on a line to minimize the error vis-a-vis the geometry. If this is a region where the solution changes fast this not so refined element will be refined allowing an improved description of the geometry.
Or have I missed something important? Such as that it will not work in 3D... Joa On Tue, Oct 13, 2009 at 05:23:07PM -0500, Roy Stogner wrote: > > > On Tue, 13 Oct 2009, Roy Stogner wrote: > > >Anyway, the solution is probably to, after you've looped over all > >nodes, loop over all constraint equations in the dof map and apply > >them to nodes' xyz coordinates. > > Thinking about it again, this is an oversimplification. For it to > work, you'd need to have a variable (either in your system, or better > in an ExplicitSystem created for the purpose) which is of type > LAGRANGE and which is of the same FE order as your geometric Elem > order. You'd also need to have some lookup between that variable's > dof indices and your node ids; they won't be the same thing. This is > still the right solution, it's just easier said than done. > --- > 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
