Re: [deal.II] Re: Problem with parallelization when using hyper_cube_slit

2018-04-17 Thread Roberto Porcù
Dear Timo. thank you very much. I removed the check on the cell when setting the boundary indicators and now it is working correctly in parallel. Thank you. Best, Roberto On Monday, April 16, 2018 at 10:35:05 AM UTC-4, Timo Heister wrote: > > Roberto, > > you are setting boundary ids only

Re: [deal.II] Re: Problem with parallelization when using hyper_cube_slit

2018-04-16 Thread Timo Heister
Roberto, you are setting boundary ids only for locally owned cells. I assume this means your constraint matrix will be inconsistent because of that. Can you check what happens if you use the correct boundary indices on all processors? On Fri, Apr 13, 2018 at 1:42 PM, Roberto Porcù wrote: > Dea

[deal.II] Re: Problem with parallelization when using hyper_cube_slit

2018-04-13 Thread Roberto Porcù
Dear Lucas, thanks very much for replying. I tried to solve the same problem without rotating the domain. I rotated the problem. Unfortunately the issue is still there. I attach the new images. I also attach the source code again because I realized that the previous one had a misplaced bracke

[deal.II] Re: Problem with parallelization when using hyper_cube_slit

2018-04-13 Thread Lucas Campos
Dear Roberto, I think the problem is that you are using GridTools::transform (via GridTools::rotate), in a distributed triangulation. This should not be done, as noted here: https://www.dealii.org/8.5.0/doxygen/deal.II/namespaceGridTools.html#a212e99cf0d923cebfa04f1d23fa60b04 If this is indee

[deal.II] Re: Problem with parallelization when using hyper_cube_slit

2018-04-11 Thread Roberto Porcù
Dear all, in the reply to my first post I forgot the code. I attach it to this reply. Best Roberto P. On Tuesday, February 20, 2018 at 7:27:37 PM UTC-5, Roberto Porcù wrote: > > Dear all, > > I'm solving a linear elasticity problem on a cracked domain created by > means > of the GridGenera

[deal.II] Re: Problem with parallelization when using hyper_cube_slit

2018-04-11 Thread Roberto Porcù
Dear all, since I am still experiencing the same issue, I tried to solve a very simple and basic problem on a domain created with the GridGenerator::hyper_cube_slit and then rotated of -90 degrees. As you can see in the attached figures I am experiencing the same issue: when I run the prog