If we use the old velocity, it can happen that we hit the kernel of
the singular norm component

This will allow us to compute the initial acceleration

As a consequence, the nonlinearity is only assembled once

Kill tests that use the old interface

Also fix the order a-mu-V0-stress-b-state-L everywhere