Multigrid with adaptive refinement and inner face batches in matrix-free

36 views
Skip to first unread message

Davit Gyulamiryan

unread,
Dec 12, 2025, 5:00:46 AM (6 days ago) Dec 12
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 fo ther multigrid 
levels:
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."));

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

Davit Gyulamiryan

unread,
Dec 15, 2025, 6:48:55 AM (3 days ago) Dec 15
to deal.II User Group
Edit: I am also using matrix-free operators some of which require iterating over the inner cell batches  ->  I am also using matrix-free operators some of which require iterating over the inner face batches.
Reply all
Reply to author
Forward
0 new messages