Domain Decomposition preconditioning for high-frequency Helmholtz problems with absorption

In this paper we give new results on domain decomposition preconditioners for GMRES when computing piecewise-linear finite-element approximations of the Helmholtz equation $-\Delta u - (k^2+ {\rm i} \varepsilon)u = f$, with absorption parameter $\varepsilon \in \mathbb{R}$. Multigrid approximations of this equation with $\varepsilon \not= 0$ are commonly used as preconditioners for the pure Helmholtz case ($\varepsilon = 0$). However a rigorous theory for such (so-called "shifted Laplace") preconditioners, either for the pure Helmholtz equation, or even the absorptive equation ($\varepsilon \not=0$), is still missing. We present a new theory for the absorptive equation that provides rates of convergence for (left- or right-) preconditioned GMRES, via estimates of the norm and field of values of the preconditioned matrix. This theory uses a $k$- and $\varepsilon$-explicit coercivity result for the underlying sesquilinear form and shows, for example, that if $|\varepsilon|\sim k^2$, then classical overlapping additive Schwarz will perform optimally for the absorptive problem, provided the subdomain and coarse mesh diameters are carefully chosen. Extensive numerical experiments are given that support the theoretical results. The theory for the absorptive case gives insight into how its domain decomposition approximations perform as preconditioners for the pure Helmholtz case $\varepsilon = 0$. At the end of the paper we propose a (scalable) multilevel preconditioner for the pure Helmholtz problem that has an empirical computation time complexity of about $\mathcal{O}(n^{4/3})$ for solving finite element systems of size $n=\mathcal{O}(k^3)$, where we have chosen the mesh diameter $h \sim k^{-3/2}$ to avoid the pollution effect. Experiments on problems with $h\sim k^{-1}$, i.e. a fixed number of grid points per wavelength, are also given

Product quasi-interpolation in logarithmically singular integral equations

A discrete high order method is constructed and justified for a class of Fredholm integral equations of the second kind with kernels that may have boundary and logarithmic diagonal singularities. The method is based on the improving the boundary behaviour of the kernel with the help of a change of variables, and on the product integration using quasi-interpolation by smooth splines of order m. Properties of different proposed calculation schemes are compared through numerical experiments using, in particular, variable precision interval arithmetics

About the convergence type of improper integrals defining fractional derivatives

This article continues the analysis of the class of fractionally differentiable functions. We complete the main result of [4] that characterises the class of fractionally differentiable functions in terms of the pointwise convergence of certain improper integrals containing these functions. Our aim is to present an example, which shows that in order to obtain all fractionally differentiable functions, one may not replace the conditional convergence of those integrals by their absolute convergence

Generalized Qualification and Qualification Levels for Spectral Regularization Methods

The concept of qualification for spectral regularization methods for inverse ill-posed problems is strongly associated to the optimal order of convergence of the regularization error. In this article, the definition of qualification is extended and three different levels are introduced: weak, strong and optimal. It is shown that the weak qualification extends the definition introduced by Mathe and Pereverzev in 2003, mainly in the sense that the functions associated to orders of convergence and source sets need not be the same. It is shown that certain methods possessing infinite classical qualification, e.g. truncated singular value decomposition (TSVD), Landweber's method and Showalter's method, also have generalized qualification leading to an optimal order of convergence of the regularization error. Sufficient conditions for a SRM to have weak qualification are provided and necessary and sufficient conditions for a given order of convergence to be strong or optimal qualification are found. Examples of all three qualification levels are provided and the relationships between them as well as with the classical concept of qualification and the qualification introduced by Mathe and Perevezev are shown. In particular, spectral regularization methods having extended qualification in each one of the three levels and having zero or infinite classical qualification are presented. Finally several implications of this theory in the context of orders of convergence, converse results and maximal source sets for inverse ill-posed problems, are shown.Comment: 20 pages, 1 figur

Методы решения линейных некорректно поставленных задач в гильбертовых пространствах : учебное пособие для студентов математического факультета

• Предисловие • Оглавление • Основные понятия теории некорректно поставленных задач • § 1. Построение класса методов решения некорректно поставленных задач • § 2. Примеры методов • § 3. Оценки разности степеней операторов • § 4. Априорное задание параметра регуляризации • § 5. Принцип невязки: самосопряженная задача • § 6. принцип невязки: несамосопряжеая задача • § 7. Критический уровень невязки • Библиографические замечания • Литература • Предметный указательhttp://tartu.ester.ee/record=b1172450~S1*es

Twitteri korpus ajas muutuva keelekasutuse uurimiseks

https://www.ester.ee/record=b5506700*es

Regularisierung nichtkorrekter Aufgaben

Das sind die Texte der Vorlesungen, die ich im Dezember 1988 - März 1989 an der Universität Kaiserslautern hielt. Die Sektionen 1-4 enthalten Materialien, die in Russisch im Buch [33] und in früheren Arbeiten [27,28] [30-33] publiziert sind. Sektion 5 enthält neue Ergebnisse, die wir während meines Aufenthaltes in Kaiserslautern in Zusammenarbeit mit Herrn Robert Plato (TU Berlin) ausarbeiteten (siehe [21,22]). Sektion 6 ist eine Erweiterung der Arbeit [31]