[dune-fem] assert failed in adaptiveleafindexset.hh

[Dedner, Andreas]

Hi Thomas.

Tut mr leid, dass ich so spaet erst reagiere und ich kann auch im Moment noch nichts
wirklich sagen...
Gibts irgenteine Art das einfach zu rekonstroieren. Was ist denn der GridPart?

Gruss
Andreas
________________________________________
From: dune-fem-bounces+a.s.dedner=warwick.ac.uk at dune-project.org [dune-fem-bounces+a.s.dedner=warwick.ac.uk at dune-project.org] on behalf of Thomas Albrecht [Thomas.Albrecht at unibas.ch]
Sent: Tuesday, May 22, 2012 2:41 PM
To: Tobias Malkmus
Cc: dune-fem at dune-project.org
Subject: Re: [dune-fem] assert failed in adaptiveleafindexset.hh

Hi,

it is a parallel ALUCubeGrid<3,3>. I did do a load balancing after
creation and none of the grids is empty.

I also looked at the leafIndex_ vector, from which the index is taken.
There are many numbers in there, quite a few are -1, but that is what
the vector is initialized to. So the question is, how it can happen that

return leafIndex_[ entity ].first;

can return one of the -1 values.

Thomas

On 05/22/2012 03:35 PM, Tobias Malkmus wrote:
> Hi Thomas
>
> Which grid type are you using ?
>
> Maybe one process has no elements ?
> Then a load balancing after creation of the grid might help.
>
>
> Best Tobias
>
> On 05/22/2012 03:21 PM, Thomas Albrecht wrote:
>> Hi everybody,
>>
>> I am still working on an image registration algorithm, which is an
>> alteration of the dgtest tutorial program. I'm using dune-fem-1.2.1
>>
>> I am getting a segmentation fault when I run the program with more than
>> two processes. After compiling in Debug mode I found out that only after
>> the first adaptive grid refinement I get an assert failed error in line
>> 439:
>>
>> 438: const int idx = codimSet.index( entity );
>> 439: assert( (idx >= 0) && (idx < codimSet.size()) );
>>
>> This is because the line the idx is -1 . It is returned from the
>> function:
>>
>> //! return leaf index for given entity
>> template <class EntityType>
>> int index ( const EntityType& entity ) const
>> {
>> assert( myCodim_ == EntityType :: codimension );
>> return leafIndex_[ entity ].first;
>> }
>>
>> in codimindexset.hh. By the way: codim is 0.
>>
>> So obviously I have no idea what's going on here, but it must have
>> something to do with the grid adaptation with several processes, which
>> provokes an invalid index to be returned.
>>
>> Do you have any idea?
>>
>> Thanks
>>
>> Thomas
>>
>
>


--
Thomas Albrecht, PhD
Postdoctoral Researcher
Department of Mathematics and Computer Science
University of Basel
Bernoullistrasse 16
4056 Basel
Switzerland

tel: +41 61 267 0442
email: thomas.albrecht at unibas.ch


_______________________________________________
dune-fem mailing list
dune-fem at dune-project.org
http://lists.dune-project.org/mailman/listinfo/dune-fem