“…Numerical methods involved finite-difference approach [19,21,22], bi-k-Lagrange elements [23], spectral collocation methods with Chebyshev polynomials of the first and second kind [24] as well as basis set expansion technique [25]. Time-dependent equations were solved with implicit and semi-implicit Crank-Nicolson methods [26,27,18,28,29], Euler scheme [22], third and fourth-order adaptive Runge-Kutta methods [30], splitstep finite difference method [22] and time-splitting sine and Fourier pseudospectral methods [31,32]. In the latter case, space was discretized with second-and fourth-order finite differences, exponential splines [29] or with Chebyshev-Tau spectral discretization method [26].…”