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 iteration using a solution transfer to interpolate the vertices to the refined mesh. After a lot of digging I found the method communicate_locally_moved_vertices with a note that I should convert the triangulation back to its original shape, then apply refinement, and then move the vertices to the desired position. This simple trick works great so far. Kudos to whoever came up with that !! Giorgos On Thursday, February 15, 2018 at 3:54:03 PM UTC-8, Wolfgang Bangerth wrote: > > 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 yet, although it doesnt seem as > > an artifact. I would say that it looks like that the one processor > > doesnt know that there is a boundary in this ghost cell. > > At the moment I set boundary ids only for the locally owned cells. > > Not good enough. Both the IndexSet you pass to ConstrainedMatrix and > VectorTools::interpolate_boundary_values() needs to have correct > information about ghost cells. So the index set needs to contain all > locally relevant DoFs, and you need to set boundary_ids for all ghost > cells as well. > > > > I'd like to know if this is somehow expected when more than one > > processor is used. > > No -- the solution you compute needs to be exactly the same (up to > round-off and solver tolerance, of course) whether you compute with one > or multiple processors. > > Best > W. > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: [email protected] > <javascript:> > www: http://www.math.colostate.edu/~bangerth/ > -- 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 [email protected]. For more options, visit https://groups.google.com/d/optout.
