English

Rank structured approximation method for quasi--periodic elliptic problems

Numerical Analysis 2017-01-03 v1

Abstract

We consider an iteration method for solving an elliptic type boundary value problem Au=f\mathcal{A} u=f, where a positive definite operator A\mathcal{A} is generated by a quasi--periodic structure with rapidly changing coefficients (typical period is characterized by a small parameter ϵ\epsilon) . The method is based on using a simpler operator A0\mathcal{A}_0 (inversion of A0\mathcal{A}_0 is much simpler than inversion of A\mathcal{A}), which can be viewed as a preconditioner for A\mathcal{A}. We prove contraction of the iteration method and establish explicit estimates of the contraction factor qq. Certainly the value of qq depends on the difference between A\mathcal{A} and A0\mathcal{A}_0. For typical quasi--periodic structures, we establish simple relations that suggest an optimal A0\mathcal{A}_0 (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\mathcal{A} admit low rank representations and algebraic operations are performed in tensor structured formats. Under moderate assumptions the storage and solution complexity of our approach depends only weakly (merely linear-logarithmically) on the frequency parameter 1/ϵ1/\epsilon, providing the FEM approximation of the order of O(ϵ1+p)O(\epsilon^{1+p}), p>0p>0.

Keywords

Cite

@article{arxiv.1701.00039,
  title  = {Rank structured approximation method for quasi--periodic elliptic problems},
  author = {B. Khoromskij and S. Repin},
  journal= {arXiv preprint arXiv:1701.00039},
  year   = {2017}
}

Comments

23 pages, 15 figures

R2 v1 2026-06-22T17:38:09.833Z