random question, while catching up with some linalg (matrix decomposition and least squares with reduced rank) scipy has rank-revealing (pivoting) QR for a few versions now. Would it be worth it to add a scipy version of matrix_rank that uses QR, instead of SVD as in numpy (IIRC) What's the cutoff for a QR version? Josef