U, s, V = np.linalg.svd(M, full_matrices=False)