Moving to devel list. On Apr 23, 2010, at 8:46 AM, Lorenzo Botti wrote:
> Some time ago Derek pointed me to the nemesis library in order to > have a > full parallel IO in libMesh. > > To start smoothly with the exodusII format I've implemented output for > discontinuous solution, see the attached patch. > It seems to work with ex21... Did anyone receive this patch? I don't see the attachment. > I would like to remove the limitation of having a single element > type mesh > and then move to the output in the nemesis format and maybe add the > possibility to output h-refined meshes in case on high orders > approximations. Note that the limitation isn't that it's a single element type mesh: it's that one "block" or "subdomain" can only contain one type of element. This is fairly ingrained in exodus (you literally can't have more than one element type in a block with exodus) but we could come up with a scheme for dealing with it (like having "virtual" blocks that are just the regular block number plus 1000 or something... and then we reconstruct the original subdomain ids when we read the mesh). I'm open to ideas. 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 didn't manage to compile without including "equation_system.h" in > exodusII_io.h... Any help is appreciated. I'll take a look when I can see the attachment. Derek ------------------------------------------------------------------------------ _______________________________________________ Libmesh-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
