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?

link|improve this question

1  
I haven't analyzed or experimented with this specific case, but certainly the use of ADI can in general affect the stable timestep. However, in the case of the implicit trapezoidal rule, I wouldn't expect any change as it is A-stable. – David Ketcheson Feb 13 at 6:40
@David Ketcheson: You may want to write that up as an answer. – Dan Feb 16 at 17:49
feedback

Know someone who can answer? Share a link to this question via email, Google+, Twitter, or Facebook.

Your Answer

 
or
required, but never shown

Browse other questions tagged or ask your own question.