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

Alexander Ostermann

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

ResearcherUnclaimed source record
math.NANumerical Analysisastro-ph.GAastro-ph.IMComputer VisionGraphicsInformation Theorymath.CAmath.ITphysics.comp-phphysics.plasm-ph

Catalog footprint

What is connected

14works
11topics
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

14 published item(s)

preprint2021arXiv

A $μ$-mode integrator for solving evolution equations in Kronecker form

In this paper, we propose a $μ$-mode integrator for computing the solution of stiff evolution equations. The integrator is based on a $d$-dimensional splitting approach and uses exact (usually precomputed) one-dimensional matrix exponentials. We show that the action of the exponentials, i.e. the corresponding batched matrix-vector products, can be implemented efficiently on modern computer systems. We further explain how $μ$-mode products can be used to compute spectral transforms efficiently even if no fast transform is available. We illustrate the performance of the new integrator by solving, among the others, three-dimensional linear and nonlinear Schrödinger equations, and we show that the $μ$-mode integrator can significantly outperform numerical methods well established in the field. We also discuss how to efficiently implement this integrator on both multi-core CPUs and GPUs. Finally, the numerical experiments show that using GPUs results in performance improvements between a factor of $10$ and $20$, depending on the problem.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A low-rank projector-splitting integrator for the Vlasov--Maxwell equations with divergence correction

The Vlasov--Maxwell equations are used for the kinetic description of magnetized plasmas. As they are posed in an up to 3+3 dimensional phase space, solving this problem is extremely expensive from a computational point of view. In this paper, we exploit the low-rank structure in the solution of the Vlasov equation. More specifically, we consider the Vlasov--Maxwell system and propose a dynamic low-rank integrator. The key idea is to approximate the dynamics of the system by constraining it to a low-rank manifold. This is accomplished by a projection onto the tangent space. There, the dynamics is represented by the low-rank factors, which are determined by solving lower-dimensional partial differential equations. The proposed scheme performs well in numerical experiments and succeeds in capturing the main features of the plasma dynamics. We demonstrate this good behavior for a range of test problems. The coupling of the Vlasov equation with the Maxwell system, however, introduces additional challenges. In particular, the divergence of the electric field resulting from Maxwell's equations is not consistent with the charge density computed from the Vlasov equation. We propose a correction based on Lagrange multipliers which enforces Gauss' law up to machine precision.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A splitting approach for the magnetic Schrödinger equation

The Schrödinger equation in the presence of an external electromagnetic field is an important problem in computational quantum mechanics. It also provides a nice example of a differential equation whose flow can be split with benefit into three parts. After presenting a splitting approach for three operators with two of them being unbounded, we exemplarily prove first-order convergence of Lie splitting in this framework. The result is then applied to the magnetic Schrödinger equation, which is split into its potential, kinetic and advective parts. The latter requires special treatment in order not to lose the conservation properties of the scheme. We discuss several options. Numerical examples in one, two and three space dimensions show that the method of characteristics coupled with a nonequispaced fast Fourier transform (NFFT) provides a fast and reliable technique for achieving mass conservation at the discrete level.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Splitting methods for constrained diffusion-reaction systems

We consider Lie and Strang splitting for the time integration of constrained partial differential equations with a nonlinear reaction term. Since such systems are known to be sensitive with respect to perturbations, the splitting procedure seems promising as we can treat the nonlinearity separately. This has some computational advantages, since we only have to solve a linear constrained system and a nonlinear ODE. However, Strang splitting suffers from order reduction which limits its efficiency. This is caused by the fact that the nonlinear subsystem produces inconsistent initial values for the constrained subsystem. The incorporation of an additional correction term resolves this problem without increasing the computational cost. Numerical examples including a coupled mechanical system illustrate the proven convergence results.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Leja method revisited: backward error analysis for the matrix exponential

The Leja method is a polynomial interpolation procedure that can be used to compute matrix functions. In particular, computing the action of the matrix exponential on a given vector is a typical application. This quantity is required, e.g., in exponential integrators. The Leja method essentially depends on three parameters: the scaling parameter, the location of the interpolation points, and the degree of interpolation. We present here a backward error analysis that allows us to determine these three parameters as a function of the prescribed accuracy. Additional aspects that are required for an efficient and reliable implementation are discussed. Numerical examples that illustrate the performance of our Matlab code are included.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fast algorithms for morphological operations using run-length encoded binary images

This paper presents innovative algorithms to efficiently compute erosions and dilations of run-length encoded (RLE) binary images with arbitrary shaped structuring elements. An RLE image is given by a set of runs, where a run is a horizontal concatenation of foreground pixels. The proposed algorithms extract the skeleton of the structuring element and build distance tables of the input image, which are storing the distance to the next background pixel on the left and right hand sides. This information is then used to speed up the calculations of the erosion and dilation operator by enabling the use of techniques which allow to skip the analysis of certain pixels whenever a hit or miss occurs. Additionally the input image gets trimmed during the preprocessing steps on the base of two primitive criteria. Experimental results show the advantages over other algorithms. The source code of our algorithms is available in C++.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Splitting methods for time integration of trajectories in combined electric and magnetic fields

The equations of motion of a single particle subject to an arbitrary electric and a static magnetic field form a Poisson system. We present a second-order time integration method which preserves well the Poisson structure and compare it to commonly used algorithms, such as the Boris scheme. All the methods are represented in a general framework of splitting methods. We use the so-called $ϕ$ functions, which give efficient ways for both analyzing and implementing the algorithms. Numerical experiments show an excellent long term stability for the new method considered.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An almost symmetric Strang splitting scheme for the construction of high order composition methods

In this paper we consider splitting methods for nonlinear ordinary differential equations in which one of the (partial) flows that results from the splitting procedure can not be computed exactly. Instead, we insert a well-chosen state $y_{\star}$ into the corresponding nonlinearity $b(y)y$, which results in a linear term $b(y_{\star})y$ whose exact flow can be determined efficiently. Therefore, in the spirit of splitting methods, it is still possible for the numerical simulation to satisfy certain properties of the exact flow. However, Strang splitting is no longer symmetric (even though it is still a second order method) and thus high order composition methods are not easily attainable. We will show that an iterated Strang splitting scheme can be constructed which yields a method that is symmetric up to a given order. This method can then be used to attain high order composition schemes. We will illustrate our theoretical results, up to order six, by conducting numerical experiments for a charged particle in an inhomogeneous electric field, a post-Newtonian computation in celestial mechanics, and a nonlinear population model and show that the methods constructed yield superior efficiency as compared to Strang splitting. For the first example we also perform a comparison with the standard fourth order Runge--Kutta methods and find significant gains in efficiency as well better conservation properties.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Convergence analysis of Strang splitting for Vlasov-type equations

A rigorous convergence analysis of the Strang splitting algorithm for Vlasov-type equations in the setting of abstract evolution equations is provided. It is shown that under suitable assumptions the convergence is of second order in the time step τ. As an example, it is verified that the Vlasov-Poisson equation in 1+1 dimensions fits into the framework of this analysis. Also, numerical experiments for the latter case are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Explicit exponential Runge-Kutta methods of high order for parabolic problems

Exponential Runge-Kutta methods constitute efficient integrators for semilinear stiff problems. So far, however, explicit exponential Runge-Kutta methods are available in the literature up to order 4 only. The aim of this paper is to construct a fifth-order method. For this purpose, we make use of a novel approach to derive the stiff order conditions for high-order exponential methods. This allows us to obtain the conditions for a method of order 5 in an elegant way. After stating the conditions, we first show that there does not exist an explicit exponential Runge-Kutta method of order 5 with less than or equal to 6 stages. Then, we construct a fifth-order method with 8 stages and prove its convergence for semilinear parabolic problems. Finally, a numerical example is given that illustrates our convergence bound.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Exponential Integrators on Graphic Processing Units

In this paper we revisit stencil methods on GPUs in the context of exponential integrators. We further discuss boundary conditions, in the same context, and show that simple boundary conditions (for example, homogeneous Dirichlet or homogeneous Neumann boundary conditions) do not affect the performance if implemented directly into the CUDA kernel. In addition, we show that stencil methods with position-dependent coefficients can be implemented efficiently as well. As an application, we discuss the implementation of exponential integrators for different classes of problems in a single and multi GPU setup (up to 4 GPUs). We further show that for stencil based methods such parallelization can be done very efficiently, while for some unstructured matrices the parallelization to multiple GPUs is severely limited by the throughput of the PCIe bus.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Analysis of exponential splitting methods for inhomogeneous parabolic equations

We analyze the convergence of the exponential Lie and exponential Strang splitting applied to inhomogeneous second-order parabolic equations with Dirichlet boundary conditions. A recent result on the smoothing properties of these methods allows us to prove sharp convergence results in the case of homogeneous Dirichlet boundary conditions. When no source term is present and natural regularity assumptions are imposed on the initial value, we show full-order convergence of both methods. For inhomogeneous equations, we prove full-order convergence for the exponential Lie splitting, whereas order reduction to 1.25 for the exponential Strang splitting. Furthermore, we give sufficient conditions on the inhomogeneity for full-order convergence of both methods. Moreover our theoretical convergence results explain the severe order reduction to 0.25 of splitting methods applied to problems involving inhomogeneous Dirichlet boundary conditions. We include numerical experiments to underline the sharpness of our theoretical convergence results.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Convergence analysis of a discontinuous Galerkin/Strang splitting approximation for the Vlasov--Poisson equations

A rigorous convergence analysis of the Strang splitting algorithm with a discontinuous Galerkin approximation in space for the Vlasov--Poisson equations is provided. It is shown that under suitable assumptions the error is of order $\mathcal{O}(τ^2+h^q +h^q / τ)$, where $τ$ is the size of a time step, $h$ is the cell size, and $q$ the order of the discontinuous Galerkin approximation. In order to investigate the recurrence phenomena for approximations of higher order as well as to compare the algorithm with numerical results already available in the literature a number of numerical simulations are performed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Investigation of the recombination of the retarded shell of "born-again" CSPNe by time-dependent radiative transfer models

A standard planetary nebula stays more than 10 000 years in the state of a photoionized nebula. As long as the timescales of the most important ionizing processes are much smaller, the ionization state can be characterized by a static photoionization model and simulated with codes like CLOUDY (Ferland et al. 1998). When the star exhibits a late Helium flash, however, its ionizing flux stops within a very short period. The star then re-appears from itsopaque shell after a few years (or centuries) as a cold giant star without any hard ionizing photons. Describing the physics of such behavior requires a fully time-dependent radiative transfer model. Pollacco (1999), Kerber et al. (1999) and Lechner & Kimeswenger (2004) used data of the old nebulae around V605 Aql and V4334 Sgr to derive a model of the pre-outburst state of the CSPN in a static model. Their argument was the long recombination time scale for such thin media. With regard to these models Schoenberner (2008) critically raised the question whether a significant change in the ionization state (and thus the spectrum) has to be expected after a time of up to 80 years, and whether static models are applicable at all.

0 citations0 reviewsNo code linkedOpen access