
<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hi guys,<div class=""><br class=""></div><div class="">Sorry to bother you again :-) I’m trying to set up my transient DG solver in parallel on both UG and YaspGrid. The latter works out well. I create an OverlappingEntitySet, and a QkDGFEM with P0ParallelConstraints, and I use the OVLP_AMG_4_DG linear solver.</div><div class=""><br class=""></div><div class="">Now I want to find a setup for the unstructured UG grid with ghosts. I use grids with simplices or rectangles exclusively, which means that I either use a PkFEM or a QkDGFEM. I try to follow the Changelog instructions of PDELab v2.4, and the non-overlapping solver examples provided in the ’dune/pdelab/test’ directory, but it doesn’t work out.</div><div class=""><br class=""></div><div class="">I want to use the NOVLP_BCGS_AMG_SSOR solver. As already noted in [MR !251](<a href="https://gitlab.dune-project.org/pdelab/dune-pdelab/merge_requests/251" class="">https://gitlab.dune-project.org/pdelab/dune-pdelab/merge_requests/251</a>), the OneStepGridOperator does not supply the required makeConsistent() function. So I plug in my first GridOperator instead.</div><div class=""><br class=""></div><div class="">In the Changelog, it’s stated that ‘Nonoverlapping computations should now use `ConformingDirichletConstraints` instead of `NonoverlappingConformingDirichletConstraints`, as there are no ghost DOFs that need to be constrained anymore.’.</div><div class="">Is this also true for a DG setup? The PDELab non-overlapping tests all assemble a GFS with QkFEM and no constraints. So does this FEM handle restricting the DOFs at processor boundaries itself?</div><div class=""><br class=""></div><div class="">I use the NonOverlappingEntitySet, as that's explicitly necessary for non-overlapping computations.</div><div class=""><br class=""></div><div class="">PkFEM works fine when run on one process.</div><div class="">However, the linear solver throws an error if run on more than one process. The backtrace is bt_simplex_novlp_ghostcc.txt.</div><div class=""><br class=""></div><div class="">QkDGFEM throws the same segfault if either executed sequentially or in parallel. The backtrace is bt_cube_novlp_ghostcc.txt.</div><div class=""><br class=""></div><div class="">Both errors are independent of whether I use NoConstraints, P0ParallelGhostConstraints or P0ParallelConstraints.</div><div class=""><br class=""></div><div class="">Can you help me out on this one?</div><div class="">Thanks and cheers,</div><div class="">Lukas</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""></div></body></html>