solve_LS_fast(U,S,V,XB) fails unit testing with XB.rows() > 2, see LinearAlgebra.cpp:90