A new algorithm for computing the Geronimus transformation with large shifts