Rank structured approximation method for quasi-periodic elliptic problems

From MaRDI portal




Abstract: We consider an iteration method for solving an elliptic type boundary value problem mathcalAu=f, where a positive definite operator mathcalA 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 mathcalA0 (inversion of mathcalA0 is much simpler than inversion of mathcalA), which can be viewed as a preconditioner for mathcalA. 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 mathcalA and mathcalA0. For typical quasi--periodic structures, we establish simple relations that suggest an optimal mathcalA0 (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 mathcalA 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/epsilon, providing the FEM approximation of the order of O(epsilon1+p), p>0.



Cites work







This page was built for publication: Rank structured approximation method for quasi-periodic elliptic problems

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1692721)