I've been running in devel mode. Things start to look complicated... What about this as an approach, escentially what was proposed to me in the beginning.
1) Use gmsh to generate a mesh that is fine enough to describe the geometry to the precision it is known 2) Solve the equations 3) Refine and COARSE the mesh using libmesh, nothing fancy 4) -> 2 until I've gotten rid of enough cells... This should work "out of the box", no? cheers 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
