Add partial pivoting runtime option to LU.

Note: in fact, inverse() always uses partial pivoting because the algo currently used doesn't make sense with complete pivoting.
No num stability issue so far even with size 200x200.

If there is any problem we can of course reimplement inverse on top of LU.