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]
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.
