Re: [deal.II] Non smooth solution when the BC coincide with the processors boundaries?

2018-02-16 Thread Giorgos Kourakos
Thank you for you reply, The issue with the boundaries was indeed solved when I assigned the boundary ids on the ghost cells as well. As for the other problem with the refinement around the processor borders, I was using a relatively complex way to move the triangulation vertices each iteratio

Re: [deal.II] Non smooth solution when the BC coincide with the processors boundaries?

2018-02-15 Thread Wolfgang Bangerth
On 02/15/2018 04:48 PM, Giorgos Kourakos wrote: I haven't forgotten the ConstraintMatrix::distribute. The code works with 1 processor I do set manually boundary ids but I reset them after every refinement. Isn't this enough to get the correct boundaries? I haven't cross out the other two cases

Re: [deal.II] Non smooth solution when the BC coincide with the processors boundaries?

2018-02-15 Thread Timo Heister
Giorgos, there can be several reasons for this: 1. You forgot a ConstraintMatrix::distribute() call (after solving) 2. Your ConstraintMatrix is inconsistent (for example if you manually set boundary ids or values in the ConstraintMatrix) 3. If you do some kind of nonlinear iteration, you might not

[deal.II] Non smooth solution when the BC coincide with the processors boundaries?

2018-02-15 Thread Giorgos Kourakos
Hi deal.ii experts, I have the following issue and I'd like if you could give me directions what to look for. I'm attaching a figure from two consecutive iterations: The black line shows how the triangulation is split to different processors. The red line indicates a Dirichlet boundary conditi