On 9/17/20 9:07 AM, Gabriel Stankiewicz wrote:
I created IndexSets by using the function
DoFTools::locally_owned_dofs_per_component() and then gathering all indices
corresponding to displacement DoFs (three instances of IndexSet) into one
IndexSet using IndexSet::add_indices() and the fourth instance correponded to
electric potential DoFs. However, the function
DoFTools::locally_owned_dofs_per_component() wasn't returning only locally
owned DoFs but also the locally not owned ones.
No, that isn't right. This is a function that is so widely used that it is
hard to believe that it would be wrong. If you think that it returns something
wrong, I believe that your *understanding* of what this function returns is
wrong, and that might lead to further confusion downstream when you are
thinking how to initialize matrices.
There are quite a number of tutorial programs that build parallel block
matrices. I know that step-32 does this for Trilinos, but I think that the
parallel Navier-Stokes solver does it for PETSc. It would be useful to just
compare your code with how these programs do it.
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].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/7456c0c4-d042-2455-50bc-659ddaf9de16%40colostate.edu.
