"A new block Krylov subspace framework with applications to functions of matrices acting on multiple vectors"