Rank structured approximation method for quasi--periodic elliptic problems
Boris N. Khoromskij and Sergey Repin
Contact the author: Please use for correspondence this email.
Submission date: 27. Jan. 2017
published in: Computational methods in applied mathematics, 17 (2017) 3, p. 457-477
DOI number (of the published article): 10.1515/cmam-2017-0014
MSC-Numbers: 65F30, 65F50, 65N35, 65F10
Keywords and phrases: elliptic problems with quasi--periodi, precondition methods, tensor type methods, guaranteed error bounds
We consider an iteration method for solving an elliptic type boundary value problem Au = f, where a positive definite operator A is generated by a quasi–periodic structure with rapidly changing coefficients (typical period is characterized by a small parameter 𝜖) . The method is based on using a simpler operator A0 (inversion of A0 is much simpler than inversion of A), which can be viewed as a preconditioner for A. We prove contraction of the iteration method and establish explicit estimates of the contraction factor q. Certainly the value of q depends on the difference between A and A0. For typical quasi–periodic structures, we establish simple relations that suggest an optimal A0 (in a selected set of ”simple” structures) and compute the corresponding contraction factor. Further, this allows us to deduce fully computable two–sided a posteriori estimates able to control numerical solutions on any iteration. The method is especially efficient if the coefficients of A admit low rank representations and algebraic operations are performed by tensor type methods. Under moderate assumptions the storage and solution complexity of our approach depends only weakly (merely logarithmically) on the frequency parameter 1∕𝜖, providing the FEM approximation of the order of O(𝜖1+q), q > 0.
You can download the paper at arxiv.org/abs/1701.00039