Reduktionsverfahren für Differenzengleichungen bei Randwertaufgaben I

SummaryThis paper describes a fast and numerically stable method for solving the discrete Dirichlet problem for Poisson's equation in case of a rectangle (and mainly, a square). By using a special calculus for difference operators, the system of linear equations is reduced to a block-triangular system such that the diagonal blocks are heavily diagonally dominant. For a standard version of the algorithm, the number of operations and the computing time are proportional toh−2 (h=mesh width). The method is one oftotal reduction compared with the method ofblock-cyclic reduction (odd-even reduction) [2], which we describe as a method ofpartial reduction.—Due to the developed calculus, many generalizations are possible.—In a following part II of the paper, the algorithm and numerical results will be described in detail.