Hello all, I am so sorry, Wolfgang, if I failed to express my thoughts in a clear and proper manner. Let me restate my reply to J-P's comment, Hope it is better this time !
> I've attached a slightly amended version of your code that I think will > highlight what the perceived problem is. When you use GridGenerator > ::hyper_rectangle you start out with a single cell that is refined globally > 3 times, which in dim=3 results in 1*(2^3)^3 = 512 cells. For your other > grid you start off with 30*30*1 = 900 cells on the coarsest level, which > upon refinement becomes 900*(2^3)^3 = 460 800 cells. I think that you might > find that its taking a little while to produce the final refined grid. > when using subdivided_hyper_rectangle and starting with 30*30*1 = 900 cells I do not need to refine any more. Triangulation.refine_global () is required for the case of hyper_rectangle. However, I think it is not the producing of the final refined mesh that hangs the code but the add_periodicity when repetition parameter is relatively large(100 for instance). Since, it has been mentioned in the documentation that refinement of the course mesh should be applied after adding periodicity, I think the problem is that using repetition in subdivided_hyper_rectangle, I am refining the course mesh before add_periodicity. > Can you confirm that the code I've attached reaches the end, i.e. that you > see this result: > > $ mpirun -np 2 ./distib_tria_periodicity >> Adding periodicity... done. >> Refinement iteration: 1 >> Refinement iteration: 2 >> Refinement iteration: 3 >> Finished. > > The code reaches the end for small repetition parameter and hangs for large ones. I am wondering if I am right in thinking about the source of the problem. Thank you all for your time and prompt attention. Regards, Hamed -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
