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

Steffen Börm

Steffen Börm contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical AnalysisDiscrete MathematicsMathematical Software

Trust snapshot

Quick read

Trust 21 - Emerging
7works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

7 published item(s)

preprint2022arXiv

Distributed $\mathcal{H}^2$-matrices for boundary element methods

Standard discretization techniques for boundary integral equations, e.g., the Galerkin boundary element method, lead to large densely populated matrices that require fast and efficient compression techniques like the fast multipole method or hierarchical matrices. If the underlying mesh is very large, running the corresponding algorithms on a distributed computer is attractive, e.g., since distributed computers frequently are cost-effective and offer a high accumulated memory bandwidth. Compared to the closely related particle methods, for which distributed algorithms are well-established, the Galerkin discretization poses a challenge, since the supports of the basis functions influence the block structure of the matrix and therefore the flow of data in the corresponding algorithms. This article introduces distributed $\mathcal{H}^2$-matrices, a class of hierarchical matrices that is closely related to fast multipole methods and particularly well-suited for distributed computing. While earlier efforts required the global tree structure of the $\mathcal{H}^2$-matrix to be stored in every node of the distributed system, the new approach needs only local multilevel information that can be obtained via a simple distributed algorithm, allowing us to scale to significantly larger systems. Experiments show that this approach can handle very large meshes with more than $130$ million triangles efficiently.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Algorithmic Counting of Zero-Dimensional Finite Topological Spaces With Respect to the Covering Dimension

Taking the covering dimension dim as notion for the dimension of a topological space, we first specify thenumber zdim_{T_0}(n) of zero-dimensional T_0-spaces on {1,...,n}$ and the number zdim(n) of zero-dimensional arbitrary topological spaces on {1,\ldots,n} by means oftwo mappings po and P that yieldthe number po(n) of partial orders on {1,...,n} and the set P(n) of partitions of {1,...,n}, respectively. Algorithms for both mappings exist. Assuming one for po to be at hand, we use our specification of zdim_{T_0}(n) and modify one for P in such a way that it computes zdim_{T_0}(n) instead of P(n). The specification of zdim(n) then allows to compute this number from zdim_{T_0}(1) to zdim_{T_0}(n) and the Stirling numbers of the second kind S(n,1) to S(n,n). The resulting algorithms have been implemented in C and we also present results of practical experiments with them. To considerably reduce the running times for computing zdim_{T_0}(n), we also describe a backtracking approach and its parallel implementation in C using the OpenMP library.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fast large-scale boundary element algorithms

Boundary element methods (BEM) reduce a partial differential equation in a domain to an integral equation on the domain's boundary. They are particularly attractive for solving problems on unbounded domains, but handling the dense matrices corresponding to the integral operators requires efficient algorithms. This article describes two approaches that allow us to solve boundary element equations on surface meshes consisting of several millions of triangles while preserving the optimal convergence rates of the Galerkin discretization.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Approximation of integral operators by Green quadrature and nested cross approximation

We present a fast algorithm that constructs a data-sparse approximation of matrices arising in the context of integral equation methods for elliptic partial differential equations. The new algorithm uses Green's representation formula in combination with quadrature to obtain a first approximation of the kernel function and then applies nested cross approximation to obtain a more efficient representation. The resulting $\mathcal{H}^2$-matrix representation requires $\mathcal{O}(n k)$ units of storage for an $n\times n$ matrix, where $k$ depends on the prescribed accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An exact solver for simple ${\mathcal H}$-matrix systems

Hierarchical matrices (usually abbreviated ${\mathcal H}$-matrices) are frequently used to construct preconditioners for systems of linear equations. Since it is possible to compute approximate inverses or $LU$ factorizations in ${\mathcal H}$-matrix representation using only ${\mathcal O}(n \log^2 n)$ operations, these preconditioners can be very efficient. Here we consider an algorithm that allows us to solve a linear system of equations given in a simple ${\mathcal H}$-matrix format \emph{exactly} using ${\mathcal O}(n \log^2 n)$ operations. The central idea of our approach is to avoid computing the inverse and instead use an efficient representation of the $LU$ factorization based on low-rank updates performed with the well-known Sherman-Morrison-Woodbury equation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Computing the eigenvalues of symmetric H2-matrices by slicing the spectrum

The computation of eigenvalues of large-scale matrices arising from finite element discretizations has gained significant interest in the last decade. Here we present a new algorithm based on slicing the spectrum that takes advantage of the rank structure of resolvent matrices in order to compute m eigenvalues of the generalized symmetric eigenvalue problem in $\mathcal{O}(n m \log^αn)$ operations, where $α>0$ is a small constant.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Efficient arithmetic operations for rank-structured matrices based on hierarchical low-rank updates

Many matrices appearing in numerical methods for partial differential equations and integral equations are rank-structured, i.e., they contain submatrices that can be approximated by matrices of low rank. A relatively general class of rank-structured matrices are $\mathcal{H}^2$-matrices: they can reach the optimal order of complexity, but are still general enough for a large number of practical applications. We consider algorithms for performing algebraic operations with $\mathcal{H}^2$-matrices, i.e., for approximating the matrix product, inverse or factorizations in almost linear complexity. The new approach is based on local low-rank updates that can be performed in linear complexity. These updates can be combined with a recursive procedure to approximate the product of two $\mathcal{H}^2$-matrices, and these products can be used to approximate the matrix inverse and the LR or Cholesky factorization. Numerical experiments indicate that the new method leads to preconditioners that require $\mathcal{O}(n)$ units of storage, can be evaluated in $\mathcal{O}(n)$ operations, and take $\mathcal{O}(n \log n)$ operations to set up.

0 citations0 reviewsNo code linkedOpen access