On Tue, 12 Sep 2017, Cyrill Vonplanta wrote:
I am relatively new to libmesh and I would like to create a geometry in parallel with 2 cubes for benchmarking purposes. What would be the easiest way to go about this? I cannot call MeshTools::Generation::build_cube() 2 times with the same mesh, right?
I'm afraid not. Not only does build_cube() clear any existing mesh objects (and number newly created objects based on the assumption that the clear occurred), but it does its new object creation in a replicated fashion, relying on later prepare_for_use() to trim away remote elements in a distributed mesh. Fixing the first problem probably wouldn't be too hard for you (it would just require making use of a couple id offsets everywhere), but fixing the second would basically require manually partitioning and ghosting the mesh on the fly and would be quite tricky. Usually when we want a cube-type geometry in parallel we create a coarse version first (so the replicated code is cheaper) and then uniformly refine it (which can be done distributed-parallel). --- Roy ------------------------------------------------------------------------------ Check out the vibrant tech community on one of the world's most engaging tech sites, Slashdot.org! http://sdm.link/slashdot _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
