The assert at mesh/mesh_communication.C:1672
libmesh_assert (el_id[e] < mesh.n_elem());
has a chance of failing in parallel, given the algorithm:
if proc == 0
read mesh
broadcast data
prepare_for_use.
The element ids of the mesh may not be contiguous when the inequality is
asserted. In my example, the mesh has 1326 elements. In serial, a particular
element e has el_id[e] = 540. On 2 processes, it is el_id[e] = 1532.
Thanks,
Bill
------------------------------------------------------------------------------
Download Intel® Parallel Studio Eval
Try the new software tools for yourself. Speed compiling, find bugs
proactively, and fine-tune applications for parallel performance.
See why Intel Parallel Studio got high marks during beta.
http://p.sf.net/sfu/intel-sw-dev
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
