I'm using the Crank-Nicholson method to solve the time-dependent Schrödinger equation, doing the implicit part with ADI (called the split-operator method in physics circles). I'm getting some weird results that are probably the result of a bug somewhere in my code.
Just in case, I thought I should probably check to see if the method I'm using is unstable.
Does ADI change the stability properties of the Crank-Nicolson method?
If so, how?