On 2/5/24 12:26, Alex Quinlan wrote:
The dynamic sparsity pattern comes from this approach:
dof_handler_sld.distribute_dofs(fe_sld);
sld_owned_dofs = dof_handler_sld.locally_owned_dofs();
DynamicSparsityPattern debug_dsp(sld_owned_dofs);
DoFTools::make_sparsity_pattern(dof_handler_sld, debug_dsp);
It makes sense that I should have an entry at position (36,48), since
dof 36 and dof 48 share a cell. Therefore, I think my stiffness matrix
process is okay. My problem seems to be with my sparsity pattern.
So, I'm confused as to why Proc0 would not have any Sparsity Pattern
entries for (36,48) when Proc1 does indeed have (48,36).
I don't know whether that's going to fix your problem, but the canonical
place to look up how to generate sparsity patterns in parallel is
step-40, where the corresponding code looks like this:
DynamicSparsityPattern dsp(locally_relevant_dofs);
DoFTools::make_sparsity_pattern(dof_handler, dsp, constraints, false);
SparsityTools::distribute_sparsity_pattern(dsp,
dof_handler.locally_owned_dofs(),
mpi_communicator,
locally_relevant_dofs);
Best
W.
--
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/d39e972f-1975-427e-a9d1-c947e32bc02c%40colostate.edu.
