Hmmm, why does my code work in 2D but fail in 4D (and not always, only for some inputs)? Is it the curse of dimensionality... :woozy_baa:
Update: my code has nothing wrong it seems. In the unit test used a wrong boundary condition. I need unit tests for my unit tests... :woozy_baa:
