 Strang splitting

Strang splitting approximates $e^{h(A+B)}$ by $e^{hA/2}{e}^{hB}{e}^{hA/2}$. Its time symmetry cancels the quadratic local commutator defect, giving local error $O(h^3)$ and global order two.