Skip to content

Periodic boundaries broken

ALUGrid has a lot of code related to periodic boundaries but the feature seems to be broken, which I conclude from the minimal example attached. There are no intersections where intersection.neighbor() and intersection.boundary() are true. The same dgf grid file works with dune-spgrid. The result is the same for 2D and 3D.

test-periodic.cc periodic.dgf

Edited by Timo Koch