Multigrid with adaptively refined grids and inner face batches in matrix-free (Duplicate)

2 views
Skip to first unread message

Davit Gyulamiryan

unread,
Jan 19, 2026, 1:46:20 PM (16 hours ago) Jan 19
to deal.II User Group
Dear deal.ii team,

I am developing some code that uses multigrid with adaptively refined triangulations (the multigrid levels are built by passing `construct_multigrid_hierarchy` to the triangulation). I am also using matrix-free operators some of which require iterating over the inner cell batches. However, I have encountered the following message from the code raised from 
the `MatrixFere::reinint()` function when initialising the matrix-free storage for the multigrid 
levels with inner batches:
if (additional_data.mapping_update_flags_inner_faces != update_default)
Assert(face_setup.refinement_edge_faces.empty(),
    ExcNotImplemented("Setting up data structures on MG levels with "
      "hanging nodes is currently not supported."));

It seems that this is not yet implemented. I was wondering if this is something that is planned to be released soon, or if I could get some advice on how to implement it myself.

Looking forward to hearing from you.

Best wishes,
Davit Gyulamiryan

Reply all
Reply to author
Forward
0 new messages