> > > As for h-refinement for high order approximations... I like the idea, but I > don't know exactly how you would implement that with libMesh. The trouble I > have is that if we change the mesh and project the solution during output... > haven't we just mucked up our solution vector and mesh for when we go back > to solving the problem? But, I seem to remember that some mesh file formats > (GMV) has similar capability now... so I may just not understand how to do > this with libmesh. > > > >> I implemented it in a home made vtk output and my idea was not to modify the mesh. Just create a single element mesh refine it in place and output it. That approach was implemented in the mesh methods I sent you some time ago.
It is not so cpu efficient but it doesn't need a lot of memory. If you have better ideas I'm open for discussion... Thank you. Lorenzo
------------------------------------------------------------------------------
_______________________________________________ Libmesh-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
