Skip to content

Segfault in LatticeBoltzmannStep with D3Q19 free-slip

When setting up a simulation with a free-slip boundary condition in LatticeBoltzmannStep with a D3Q19 stencil, there is a segmentation fault during the run command.
Using a D2Q9 or D3Q27 stencil seems to work, as well as when the free-slip boundary condition is used outside of LatticeBoltzmannStep, i.e., by creating manually the data and boundary handling (as done in the free-slip test case).

A non-working example is attached. freeslip_d3q19_bug.py