ENH: sparse.linalg: take LinearOperator in svds

Should allow solving PCA for sparse matrix.
DOC: advertise svds(LinearOperator) + vanity