BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Folkmar Bornemann

Folkmar Bornemann contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NAmath-phmath.MPmath.PRmath.CVNumerical Analysiscond-mat.stat-mechmath.CAphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

10 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
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&#39;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&#39;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