On Fri, 26 Mar 2010, Cochran, William K. wrote: > 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.
Can you send us an example of an exodus file for which the parallel read fails? --- Roy ------------------------------------------------------------------------------ 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
