[dune-fem] assert failed in adaptiveleafindexset.hh

[Tobias Malkmus]

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
>


-- 
Tobias Malkmus                 <tomalk at mathematik.uni-freiburg.de>

Mathematisches Institut               Tel: +49 (0) 761 203 5651
Abt. für Angewandte Mathematik        Universität Freiburg
Hermann-Herder-Str. 10
79104 Freiburg