A cyclic block-tridiagonal solver

A simple algorithm for solving a cyclic block-tridiagonal system of equations is presented. Introducing a special form of a new variable, the system is split into two block-tridiagonal systems, which can be solved by known methods. Implementation details of the algorithm are discussed and numerical examples of diagonal and random generated systems are presented.