Triangular splitting implementation of RKN‐type Fourier collocation methods for second‐order differential equations

A triangular splitting implementation of Runge–Kutta–Nyström–type Fourier collocation methods is presented and analyzed in this paper. The proposed implementation relies on a reformulation of the method and on the Crout factorization of a corresponding matrix associated with the method. The excellent behavior of the splitting implementation is confirmed by its performance on a few numerical tests.