I cleaned up a TODO that was left in the new Shan Chen Tutorials.
However, I found a new problem: the multi-phase version did not consistently take the force-dependent velocity shift into account during initialization. Is this intended?
When I switch to the seemingly 'correct' initialization the test fails.