“…Applying recursively, we obtain the splitting operator formula for the discrete time t n = nh. Quite interestingly, if, for example, the exact flows ϕ [1] (t) and ϕ [2] (t) are not available, then one can use consistent approximations to these solutions. In this manner, from now on, φ [1] h and φ [2] h represents numerical discretization of the flows on the two sub-problems (P 1 ) and (P 2 ) at time t n = nh, respectively.…”