Re: [Libmesh-users] node id (fwd)

[Roy S. <roy...@ic...>]

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 <mpo...@pu...>
To: Roy Stogner <roy...@ic...>
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.