Source author record

Harry Yserentant

Harry Yserentant appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Numerical Analysis math-ph math.AP math.MP math.RA

Catalog footprint

What is connected

5works

6topics

2close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

A measure concentration effect for matrices of high, higher, and even higher dimension

Let $n>m$, and let $A$ be an $(m\times n)$-matrix of full rank. Then obviously the estimate $\|Ax\|\leq\|A\|\|x\|$ holds for the euclidean norm of $x$ and $Ax$ and the spectral norm as the assigned matrix norm. We study the sets of all $x$ for which, for fixed $δ<1$, conversely $\|Ax\|\geqδ\,\|A\|\|x\|$ holds. It turns out that these sets fill, in the high-dimensional case, almost the complete space once $δ$ falls below a bound that depends on the extremal singular values of $A$ and on the ratio of the dimensions. This effect has much to do with the random projection theorem, which plays an important role in the data sciences. As a byproduct, we calculate the probabilities this theorem deals with exactly.

preprint2022arXiv

The Laplace operator, measure concentration, Gauss functions, and quantum mechanics

We represent in this note the solutions of the electronic Schrödinger equation as traces of higher-dimensional functions. This allows to decouple the electron-electron interaction potential but comes at the price of a degenerate elliptic operator replacing the Laplace operator on the higher-dimensional space. The surprising observation is that this operator can without much loss again be substituted by the Laplace operator, the more successful the larger the system under consideration is. This is due to a concentration of measure effect that has much to do with the random projection theorem known from probability theory. The text is in parts based on the publications [Numer. Math. 146, 219--238 (2020)] and [SIAM J. Matrix Anal. Appl., 43, 464--478 (2022)] of the author and adapts the findings there to the needs of quantum mechanics. Our observations could for example find use in iterative methods that map sums of products of orbitals and geminals onto functions of the same type.

preprint2020arXiv

Numerical Homogenization of Fractal Interface Problems

We consider the numerical homogenization of a class of fractal elliptic interface problems inspired by related mechanical contact problems from the geosciences. A particular feature is that the solution space depends on the actual fractal geometry. Our main results concern the construction of projection operators with suitable stability and approximation properties. The existence of such projections then allows for the application of existing concepts from localized orthogonal decomposition (LOD) and successive subspace correction to construct first multiscale discretizations and iterative algebraic solvers with scale-independent convergence behavior for this class of problems.

preprint2020arXiv

On the expansion of solutions of Laplace-like equations into traces of separable higher dimensional functions

This paper deals with the equation $-Δu+μu=f$ on high-dimensional spaces $\mathbb{R}^m$ where $μ$ is a positive constant. If the right-hand side $f$ is a rapidly converging series of separable functions, the solution $u$ can be represented in the same way. These constructions are based on approximations of the function $1/r$ by sums of exponential functions. The aim of this paper is to prove results of similar kind for more general right-hand sides $f(x)=F(Tx)$ that are composed of a separable function on a space of a dimension $n$ greater than $m$ and a linear mapping given by a matrix $T$ of full rank. These results are based on the observation that in the high-dimensional case, for $ω$ in most of the $\mathbb{R}^n$, the euclidian norm of the vector $T^tω$ in the lower dimensional space $\mathbb{R}^m$ behaves like the euclidian norm of $ω$.

preprint2016arXiv

A Note on Approximate Inverse Iteration

Different variants of approximate inverse iteration like the locally optimal block preconditioned conjugate gradient method became in recent years increasingly popular for the solution of the large matrix eigenvalue problems arising from the discretization of selfadjoint elliptic partial differential equations, in particular for the calculation of the minimum eigenvalue. We extend in this little note the classical convergence theory of D'yakonov and Orekhov [Math. Notes 27 (1980)] to the case of operators with an essential spectrum on infinite dimensional Hilbert spaces and allow for arbitrary, sufficiently small perturbations of the solutions of the equation that links the iterates. The note complements the much more elaborate convergence theory of Neymeyr and Knyazev and Neymeyr for the matrix case (see [Knyazev and Neymeyr, SIAM J. Matrix Anal. Appl. 31 (2009)] and the references therein) and is suitable for classroom presentation.