[Dune] alugrid problem

Arya Fallahi arya.fallahi at gmail.com
Wed Jul 29 14:16:25 CEST 2015


Dear Dune,

I got to simplify my grid to a small one, which returns the above error.
Note that it always runs fine on one processor. So, there is probably a bug
in parallelization. my alugrid.cfg file looks as the following:

0.0
1.04
11

The .dgf file is also attached.

I will be thankful if you have a look at it and let me know if the bug is
from my side or from alugrid.

Regards,
Arya

On Wed, Jul 29, 2015 at 11:28 AM, Arya Fallahi <arya.fallahi at gmail.com>
wrote:

> Dear Dune,
>
> I am currently stuck with a strange parallelization problem in ALUGrid.
> When I run the code on one processor everything is fine. However, after
> running on just two processors, the following error appears:
>
> dgtd-max: /afs/
> desy.de/user/a/afallahi/workspace/dune-alugrid/dune/alugrid/common/defaultindexsets.hh:228:
> Dune::DefaultIndexSet<GridImp, IteratorImp>::IndexType
> Dune::DefaultIndexSet<GridImp, IteratorImp>::index(const typename
> GridImp::Codim<cd>::Entity&) const [with int cd = 0; GridImp =
> Dune::ALU3dGrid<(Dune::ALU3dGridElementType)4u, Dune::ALUGridMPIComm>;
> IteratorImp = Dune::EntityIterator<0, const
> Dune::ALU3dGrid<(Dune::ALU3dGridElementType)4u, Dune::ALUGridMPIComm>,
> Dune::ALU3dGridLeafIterator<0, (Dune::PartitionIteratorType)4u, const
> Dune::ALU3dGrid<(Dune::ALU3dGridElementType)4u, Dune::ALUGridMPIComm> > >;
> Dune::DefaultIndexSet<GridImp, IteratorImp>::IndexType = unsigned int;
> typename GridImp::Codim<cd>::Entity = Dune::Entity<0, 3, const
> Dune::ALU3dGrid<(Dune::ALU3dGridElementType)4u, Dune::ALUGridMPIComm>,
> Dune::ALU3dGridEntity>]: Assertion `indexContainer( cd )[ en ].index() >=
> 0' failed.
>
> There are also some strange lines in the execution of the code:
>
> DefaultIndexSet[ 0 ]: 352631 s | g 314280
> DefaultIndexSet[ 0 ]: 59350 s | g 46407
> DefaultIndexSet[ -1 ]: 59350 s | g 46407
>
> The mesh creation with gmsh runs pretty fine. There seems to be a problem
> with index sets. I have no idea about the returned error. Therefore, I can
> not begin with debugging my code. I will be thankful if you let me know
> about the possible origins of this error message.
>
> Thanks in advance
> Arya
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20150729/c6901e02/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: grid5p.dgf
Type: application/octet-stream
Size: 23554 bytes
Desc: not available
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20150729/c6901e02/attachment.obj>


More information about the Dune mailing list