BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Folkmar Bornemann

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

ResearcherUnclaimed source record
math.NAmath-phmath.MPmath.PRmath.CVNumerical Analysiscond-mat.stat-mechmath.CAphysics.comp-ph

Catalog footprint

What is connected

13works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

13 published item(s)

preprint2022arXiv

The Challenge of Sixfold Integrals: The Closed-Form Evaluation of Newton Potentials between Two Cubes

The challenge of explicitly evaluating, in elementary closed form, the weakly singular sixfold integrals for potentials and forces between two cubes has been taken up at various places in the mathematics and physics literature. It created some strikingly specific results, with an aura of arbitrariness, and a single intricate general procedure due to Hackbusch. Those scattered instances were mostly addressing the problem heads on, by successive integration while keeping track of a thicket of primitives generated at intermediate stages. In this paper we present a substantially easier and shorter approach, based on a Laplace transform of the kernel. We clearly exhibit the structure of the results as obtained by an explicit algorithm, just computing with rational polynomials. The method extends, up to the evaluation of single integrals, to higher dimensions. Among other examples, we easily reproduce Fornberg's startling closed form solution of Trefethen's two-cubes problem and Waldvogel's symmetric formula for the Newton potential of a rectangular cuboid.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Efficient numerical evaluation of thermodynamic quantities on infinite (semi-)classical chains

This work presents an efficient numerical method to evaluate the free energy density and associated thermodynamic quantities of (quasi) one-dimensional classical systems, by combining the transfer operator approach with a numerical discretization of integral kernels using quadrature rules. For analytic kernels, the technique exhibits exponential convergence in the number of quadrature points. As demonstration, we apply the method to a classical particle chain, to the semiclassical nonlinear Schrödinger equation and to a classical system on a cylindrical lattice.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A note on the expansion of the smallest eigenvalue distribution of the LUE at the hard edge

In a recent paper, Edelman, Guionnet and Péché conjectured a particular $n^{-1}$ correction term of the smallest eigenvalue distribution of the Laguerre unitary ensemble (LUE) of order $n$ in the hard-edge scaling limit: specifically, the derivative of the limit distribution, that is, the density, shows up in that correction term. We give a short proof by modifying the hard-edge scaling to achieve an optimal $O(n^{-2})$ rate of convergence of the smallest eigenvalue distribution. The appearance of the derivative follows then by a Taylor expansion of the less optimal, standard hard-edge scaling. We relate the $n^{-1}$ correction term further to the logarithmic derivative of the Bessel kernel Fredholm determinant in the work of Tracy and Widom.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the Scaling Limits of Determinantal Point Processes with Kernels Induced by Sturm-Liouville Operators

By applying an idea of Borodin and Olshanski [J. Algebra 313 (2007), 40-60], we study various scaling limits of determinantal point processes with trace class projection kernels given by spectral projections of selfadjoint Sturm-Liouville operators. Instead of studying the convergence of the kernels as functions, the method directly addresses the strong convergence of the induced integral operators. We show that, for this notion of convergence, the Dyson, Airy, and Bessel kernels are universal in the bulk, soft-edge, and hard-edge scaling limits. This result allows us to give a short and unified derivation of the known formulae for the scaling limits of the classical random matrix ensembles with unitary invariance, that is, the Gaussian unitary ensemble (GUE), the Wishart or Laguerre unitary ensemble (LUE), and the MANOVA (multivariate analysis of variance) or Jacobi unitary ensemble (JUE).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Numerical Methods for the Discrete Map $Z^a$

As a basic example in nonlinear theories of discrete complex analysis, we explore various numerical methods for the accurate evaluation of the discrete map $Z^a$ introduced by Agafonov and Bobenko. The methods are based either on a discrete Painlevé equation or on the Riemann-Hilbert method. In the latter case, the underlying structure of a triangular Riemann-Hilbert problem with a non-triangular solution requires special care in the numerical approach. Complexity and numerical stability are discussed, the results are illustrated by numerical examples

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Singular Values and Evenness Symmetry in Random Matrix Theory

Complex Hermitian random matrices with a unitary symmetry can be distinguished by a weight function. When this is even, it is a known result that the distribution of the singular values can be decomposed as the superposition of two independent eigenvalue sequences distributed according to particular matrix ensembles with chiral unitary symmetry. We give decompositions of the distribution of singular values, and the decimation of the singular values --- whereby only even, or odd, labels are observed --- for real symmetric random matrices with an orthogonal symmetry, and even weight. This requires further specifying the functional form of the weight to one of three types --- Gauss, symmetric Jacobi or Cauchy. Inter-relations between gap probabilities with orthogonal and unitary symmetry follow as a corollary. The Gauss case has appeared in a recent work of Bornemann and La Croix. The Cauchy case, when appropriately specialised and upon stereographic projection, gives decompositions for the analogue of the singular values for the circular unitary and circular orthogonal ensembles.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Singular Values of the GOE

As a unifying framework for examining several properties that nominally involve eigenvalues, we present a particular structure of the singular values of the Gaussian orthogonal ensemble (GOE): the even-location singular values are distributed as the positive eigenvalues of a Gaussian ensemble with chiral unitary symmetry (anti-GUE), while the odd-location singular values, conditioned on the even-location ones, can be algebraically transformed into a set of independent $χ$-distributed random variables. We discuss three applications of this structure: first, there is a pair of bidiagonal square matrices, whose singular values are jointly distributed as the even- and odd-location ones of the GOE; second, the magnitude of the determinant of the GOE is distributed as a product of simple independent random variables; third, on symmetric intervals, the gap probabilities of the GOE can be expressed in terms of the Laguerre unitary ensemble (LUE). We work specifically with matrices of finite order, but by passing to a large matrix limit, we also obtain new insight into asymptotic properties such as the central limit theorem of the determinant or the gap probabilities in the bulk-scaling limit. The analysis in this paper avoids much of the technical machinery (e.g. Pfaffians, skew-orthogonal polynomials, martingales, Meijer $G$-function, etc.) that was previously used to analyze some of the applications.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Automatic Deformation of Riemann-Hilbert Problems with Applications to the Painlevé II Transcendents

The stability and convergence rate of Olver's collocation method for the numerical solution of Riemann-Hilbert problems (RHPs) is known to depend very sensitively on the particular choice of contours used as data of the RHP. By manually performing contour deformations that proved to be successful in the asymptotic analysis of RHPs, such as the method of nonlinear steepest descent, the numerical method can basically be preconditioned, making it asymptotically stable. In this paper, however, we will show that most of these preconditioning deformations, including lensing, can be addressed in an automatic, completely algorithmic fashion that would turn the numerical method into a black-box solver. To this end, the preconditioning of RHPs is recast as a discrete, graph-based optimization problem: the deformed contours are obtained as a system of shortest paths within a planar graph weighted by the relative strength of the jump matrices. The algorithm is illustrated for the RHP representing the Painlevé II transcendents.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the convergence rates of Gauss and Clenshaw-Curtis quadrature for functions of limited regularity

We study the optimal general rate of convergence of the n-point quadrature rules of Gauss and Clenshaw-Curtis when applied to functions of limited regularity: if the Chebyshev coefficients decay at a rate O(n^{-s-1}) for some s > 0, Clenshaw-Curtis and Gauss quadrature inherit exactly this rate. The proof (for Gauss, if 0 < s < 2, there is numerical evidence only) is based on work of Curtis, Johnson, Riess, and Rabinowitz from the early 1970s and on a refined estimate for Gauss quadrature applied to Chebyshev polynomials due to Petras (1995). The convergence rate of both quadrature rules is up to one power of n better than polynomial best approximation; hence, the classical proof strategy that bounds the error of a quadrature rule with positive weights by polynomial best approximation is doomed to fail in establishing the optimal rate.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Optimal Contours for High-Order Derivatives

As a model of more general contour integration problems we consider the numerical calculation of high-order derivatives of holomorphic functions using Cauchy's integral formula. Bornemann (2011) showed that the condition number of the Cauchy integral strongly depends on the chosen contour and solved the problem of minimizing the condition number for circular contours. In this paper we minimize the condition number within the class of grid paths of step size h using Provan's algorithm for finding a shortest enclosing walk in weighted graphs embedded in the plane. Numerical examples show that optimal rectangular paths yield small condition numbers even in those cases where circular contours are known to be of limited use, such as for functions with branch-cut singularities.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Accuracy and Stability of Computing High-Order Derivatives of Analytic Functions by Cauchy Integrals

High-order derivatives of analytic functions are expressible as Cauchy integrals over circular contours, which can very effectively be approximated, e.g., by trapezoidal sums. Whereas analytically each radius r up to the radius of convergence is equal, numerical stability strongly depends on r. We give a comprehensive study of this effect; in particular we show that there is a unique radius that minimizes the loss of accuracy caused by round-off errors. For large classes of functions, though not for all, this radius actually gives about full accuracy; a remarkable fact that we explain by the theory of Hardy spaces, by the Wiman-Valiron and Levin-Pfluger theory of entire functions, and by the saddle-point method of asymptotic analysis. Many examples and non-trivial applications are discussed in detail.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the Numerical Evaluation of Distributions in Random Matrix Theory: A Review

In this paper we review and compare the numerical evaluation of those probability distributions in random matrix theory that are analytically represented in terms of Painlevé transcendents or Fredholm determinants. Concrete examples for the Gaussian and Laguerre (Wishart) beta-ensembles and their various scaling limits are discussed. We argue that the numerical approximation of Fredholm determinants is the conceptually more simple and efficient of the two approaches, easily generalized to the computation of joint probabilities and correlations. Having the means for extensive numerical explorations at hand, we discovered new and surprising determinantal formulae for the k-th largest (or smallest) level in the edge scaling limits of the Orthogonal and Symplectic Ensembles; formulae that in turn led to improved numerical evaluations. The paper comes with a toolbox of Matlab functions that facilitates further mathematical experiments by the reader.

0 citations0 reviewsNo code linkedOpen access