LU with partial pivoting should work (for nice mats) in c++_reduction; full pivoting is still in progress