A Nonstandard Cyclic Reduction Method, Its Variants and Stability (Q4243522)

A nonstandard cyclic reduction method is introduced for solving the Poisson equation on rectangular domains. The given system has a block tridiagonal structure \[ \left ( \begin{matrix} A & -T & & & \\ -T & A & -T & & \\ & -T & A & \ddots & \\ & & \ddots & \ddots & -T \\ & & & -T & A \end{matrix} \right) \left ( \begin{matrix} u_1 \\ u_2 \\ \vdots \\ u_{n-1} \\ u_n \end{matrix} \right) = \left ( \begin{matrix} f_1 \\ f_2 \\ \vdots \\ f_{n-1} \\ f_n \end{matrix} \right) \] In the classical formulation lines of the equations are multiplied by \(A\) and \(T\) respectively. After that the equations are summed up. In the alternative version multiplications by \(T A^{-1}\) are performed. The algorithm is compared with the Fourier method, Bunemann's algorithm, the FACR techniques, and other similar methods.