The computation of tunneling probabilities in three dimensions is a numerical challenge, because the small transition probabilities associated with the overlap of exponentially vanishing wave function-tails require large computational accuracy. In scattering situations arising, e.g., in electron tunneling in metal-molecule-metal junctions, this is compounded by the need to provide a proper trun...