Forwarding to list in case anyone runs into the same problem later ---------- Forwarded message ---------- Date: Sat, 18 Dec 2010 18:18:21 -0500 From: Michael Povolotskyi <[email protected]> To: Roy Stogner <[email protected]> Subject: Re: [Libmesh-users] node id
On 12/18/2010 12:05 AM, Roy Stogner wrote: > > > On Wed, 15 Dec 2010, Michael Povolotskyi wrote: > >> Dear Libmesh developers, >> I have a question about the node id. >> >> Is it assumed in the libmesh code that node id's are not just the unique >> numbers but rather belong to consecutive set from 0 to n? > > I think we don't make that assumption in most code, but we do renumber > nodes and elements to make them contiguous during a typical > prepare_for_use and I think many of the IO methods assume they're > contiguous. > > Also SerialMesh stores things in vectors indexed by ID; *too* much > sparsity would be a memory suck. > --- > Roy Thank you, Roy. I see it now. I was calling my_mesh->prepare_for_use() with no argument. So, the nodes were not renumbered and the output routine crashed. Michael. ------------------------------------------------------------------------------ Learn how Oracle Real Application Clusters (RAC) One Node allows customers to consolidate database storage, standardize their database environment, and, should the need arise, upgrade to a full multi-node Oracle RAC database without downtime or disruption http://p.sf.net/sfu/oracle-sfdevnl _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
