We present a family of bidirectional beam propagation methods based on the
application of Fade approximants to the square-root expressions appearing i
n the transition operator. We show that the numerical instability associate
d with the use of real Pade primes can be removed by using branch-cut rotat
ion of the square-root operator.