“…The splitting methods form an important group of methods which are quite accurate and efficient [57]. Actually, they have been widely applied for dealing with highly oscillatory systems such as the Schrödinger/nonlinear Schrödinger equations [1,8,9,22,23,55,67], the Dirac/nonlinear Dirac equations [5,6,14,54], the Maxwell-Dirac system [10,49], the Zakharov system [12,13,41,50], the Gross-Pitaevskii equation for Bose-Einstein condensation (BEC) [11], the Stokes equation [21], and the Enrenfest dynamics [32], etc.…”