<div dir="ltr"><div style="font-size:12.8000001907349px"><div><div><div><div><div><div>Dear Dune,<br><br></div>Sometime ago I upgraded my alugrid and started using the one integrated into dune library, namely dune-alugrid. I was always suffering from slow computation time that I suddenly realized today that something is going wrong with my load balancing. In the following, I am sending you an example of what I do in my code:<br><br>typedef GridSelector::GridType Grid;<br><br>MPI_Init(&argc,&argv);<br>MPI_Comm_rank(MPI_COMM_WORLD,&rank);<br>MPI_Comm_size(MPI_COMM_WORLD,&size);<br><br>GridPtr<Grid> ag ( "./grid1.dgf" , MPI_COMM_WORLD );<br>Grid& g = *ag;<br>g.loadBalance();<br>typedef LeafMultipleCodimMultipleGeomTypeMapper<Grid,P0Layout> emap;<br>emap em(g);<br>std::cout << em.size() << std::endl;<br><br></div></div>Without load balancing I get the following output on 4 processors:<br><br>0<br>0<br>0<br>15468<br><br></div>meaning that the grid has 15648 elements. However, with load balancing I get the following output:<br><br>11235<br>11869<br>12906<br>11244<br><br></div>It seems that the total number of elements are getting considerably higher. Do you have any guess what I am doing wrong?<br><br></div>Thank you very much for your help,<br></div><span style="font-size:12.8000001907349px">Arya</span><br></div>