From: Roy Stogner <roystgnr@ic...> - 2010-12-24 03:48:03
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 <mpovolot@...>
To: Roy Stogner <roystgnr@...>
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
> Also SerialMesh stores things in vectors indexed by ID; *too* much
> sparsity would be a memory suck.
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.