Results
Simple code, runs in IDL on Windows laptops and Alpha Station
With 128*128 modes, ~20 iterations in 1 second.
Convergence is very fast in a/?Də limit. About 1 second.
Convergence becomes very slow when a/?D>ɭ.
- Adaptive “time” step becomes necessary. Start with ?t=1, increase gradually to ?t~1010 or more. Convergence in ,000 iterations.