Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

make_flux_sparsity_pattern() revision #14868

Merged
merged 3 commits into from
Mar 13, 2023
Merged

make_flux_sparsity_pattern() revision #14868

Changes from 1 commit
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
8e02293
face_has_flux_coupling may be a predicate made of multiple, complicat…
vyushut Feb 11, 2023
b016359
a test is added: make_flux_sparsity_pattern() is run with a face filt…
vyushut Mar 10, 2023
44fa21d
both 2D and 3D cases with sparsity pattern printed out; the test desc…
vyushut Mar 13, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
14 changes: 8 additions & 6 deletions source/dofs/dof_tools_sparsity.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -852,9 +852,6 @@ namespace DoFTools
}
else
{
if (!face_has_flux_coupling(cell, face_n))
continue;

typename DoFHandler<dim,
spacedim>::level_cell_iterator
neighbor =
Expand Down Expand Up @@ -888,6 +885,10 @@ namespace DoFTools
neighbor->is_locally_owned())
continue; // (the neighbor is finer)

if (!face_has_flux_coupling(cell, face_n))
continue;


const unsigned int neighbor_face_n =
periodic_neighbor ?
cell->periodic_neighbor_face_no(face_n) :
Expand Down Expand Up @@ -1197,9 +1198,6 @@ namespace DoFTools
level_cell_iterator neighbor =
cell->neighbor_or_periodic_neighbor(face);

if (!face_has_flux_coupling(cell, face))
continue;

// Like the non-hp-case: If the cells are on the same
// level (and both are active, locally-owned cells)
// then only add to the sparsity pattern if the
Expand Down Expand Up @@ -1229,6 +1227,10 @@ namespace DoFTools
neighbor->is_locally_owned())
continue; // (the neighbor is finer)


if (!face_has_flux_coupling(cell, face))
continue;

// In 1d, go straight to the cell behind this
// particular cell's most terminal cell. This makes us
// skip the if (neighbor->has_children()) section
Expand Down