Change default solver for the dynamics from `solver.solve` to `solver.pinv_smooth`

Created by: PhilipVinc

As discussed with @gcarleo , it seems that this solver is more stable.