We present a detailed comparison of the efficiency and accuracy of the second- and third-order split operator methods, a time dependent modified Cayley method, and the Chebychev polynomial expansion method for solving the time dependent Schrodinger equation in the one-dimensional double well potential energy function. We also examine the efficiency and accuracy of the split operator and modified Cayley methods for the imaginary time propagation.