The aim is to compute the analytic singular value decomposition for a few selected singular values/vectors. The proposed method is suitable for large sparse matrix valued functions.