svd_with_rank

[Jorge-C]

In several cases, we follow the pattern usv = svd(a); r = rank(usv): write util function that does both things.
We don't want to use numpy.linalg.matrix_rank because it recomputes the svd (expensive). It'd probably be helpful to add an optional argument like return_rank to numpy's svd.