On Thu, 9 Dec 2010, Michael Povolotskyi wrote:
> Dear libmesh developers, > I'd like to use libmesh for solving the Poisson equation on a large grid. > I'd like to have a parallel 3D mesh. > I'd like to create the mesh myself rather than reading it from a file. > > In my program I have distributed information about nodes and elements > that I want to have in my mesh. > My question: how can I create the mesh? > > I imagine that I have to > a) reserve memory for nodes and elements > b) add points and elements > c) prepare the mesh for use > > I have done it already in for a serial case and it worked. > But can I do the same operations in parallel? ParallelMesh (as opposed to the default SerialMesh) is designed to handle the mesh storage, degree of freedom numbering, adaptive refinement, and repartitioning for this. It's currently got a serious bug with adaptive coarsening, though, and that means that it doesn't get enough usage or testing and it might have more subtle bugs elsewhere. If your coarse mesh will fit on one CPU, then the same mesh generation code can work in both serial and parallel; the prepare_for_use will just discard non-semilocal coarse elements after they've been created, and the fine elements created by refinement will only be created on the processors that need them. --- Roy ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
