Split Operator Method

This is to propagate a 2D Schrodinger equation using Split Operator method. We have used Strang Splitting to approximate the unitary operator.

The code is extremely faster compared with the matrix method (creation and annihilation operator method). For a higher grid 512X512, the matrix based method will take weeks to finish the problem, here

you can have the result within few seconds !!!

This figure is just a snapshot of the propagating wavepacket.