On 1/8/21 10:19 AM, Alberto Salvadori wrote:
thank you for the hint. I looked at setp 47 sparsity pattern. If I understood
it right,
the main difference stands in the locally_relevant_dofs, which are no longer
extracted and used.
Rather, the sparsity pattern is based on the whole dof_handler, as per the
instruction
DoFTools::make_flux_sparsity_pattern( this->dof_handler, dsp,
this->hanging_node_constraints, true );
No, you mistook what the difference is. The difference is that
make_sparsity_pattern() only adds entries for the matrix that correspond to
DoFs i and j that are located on the same cell.
make_flux_sparsity_pattern() also adds entries for DoFs i,j that are located
on cells that are face neighbors of each other.
In both cases, you can use the same locally_relevant_set when you set up the
sparsity pattern object.
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/a7c6b38b-8537-26e1-5ea6-72b06cc5b8f3%40colostate.edu.
