minor : check that A.lu().solve is give same results than A.fullPivLu().solve