An efficient ADI-solver for scattered data problems with global smoothing (Q1387878)

A minimization problem of the energy functional \[ I(u)\equiv \int_\Omega[(u_{xx})^2+ (u_{yy})^2]dx dy\qquad (\Omega\equiv (0,1)^2) \] is considered with given values of \(u\) at a subset \(D\) of grid points. The central result is connected with an iterative solution of the difference system for this problem via alternating directions implicit (ADI) methods. In the model case, \(D\) corresponds to the points on two boundary grid layers. Then, in an Euclidean space \(E\), the system can be written as \(\Lambda_1u+ \Lambda_2u= f\) with symmetric and positive operators \(\Lambda_1\) and \(\Lambda_2\); moreover, \(\Lambda_1\Lambda_2= \Lambda_2\Lambda_1\). Exactly, these systems were considered by the reviewer [Dokl. Akad. Nauk SSSR 138, 271-274 (1961; Zbl 0108.12805)] as a basis for the inner iterations; later this approach was used with success for various plate and shell problems. In this respect, the results of the paper confirm a superiority of ADI methods (for the commutative case) with respect to many other methods. Two-stage iterative methods can be applied for the general case.