“…There has been a continuous progress of developing numerical algorithms in the past few decades; see for example [6,7] and references therein. The spatial derivative is typically approximated by the finite element method (FEM) [8,9,10,11,12,13,14,15,16,17,18] and the finite difference method [19,20,21,22,23]. As for the temporal discretization, explicit schemes [15,24], fully implicit schemes [25,26,20], and semi-implicit schemes [19,27,28,29,30,31,32] have been extensively studied.…”