<div dir="ltr">Dear Dune,<div><br></div><div>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:</div><div><br></div><div><div>0.0</div><div>1.04</div><div>11</div></div><div><br></div><div>The .dgf file is also attached.</div><div><br></div><div>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.</div><div><br></div><div>Regards,</div><div>Arya</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jul 29, 2015 at 11:28 AM, Arya Fallahi <span dir="ltr"><<a href="mailto:arya.fallahi@gmail.com" target="_blank">arya.fallahi@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear Dune,<div><br></div><div>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:</div><div><br></div><div>dgtd-max: /afs/<a href="http://desy.de/user/a/afallahi/workspace/dune-alugrid/dune/alugrid/common/defaultindexsets.hh:228" target="_blank">desy.de/user/a/afallahi/workspace/dune-alugrid/dune/alugrid/common/defaultindexsets.hh:228</a>: 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.<br></div><div><br></div><div>There are also some strange lines in the execution of the code:</div><div><br></div><div><div>DefaultIndexSet[ 0 ]: 352631 s | g 314280</div><div>DefaultIndexSet[ 0 ]: 59350 s | g 46407</div><div>DefaultIndexSet[ -1 ]: 59350 s | g 46407<br></div></div><div><br></div><div>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.</div><div><br></div><div>Thanks in advance</div><span class="HOEnZb"><font color="#888888"><div>Arya</div></font></span></div>
</blockquote></div><br></div>