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

Winfried Auzinger

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

ResearcherUnclaimed source record
math.NANumerical Analysisphysics.comp-ph

Catalog footprint

What is connected

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

11 published item(s)

preprint2021arXiv

Efficient Magnus-type integrators for solar energy conversion in Hubbard models

Strongly interacting electrons in solids are generically described by Hubbardtype models, and the impact of solar light can be modeled by an additional time-dependence. This yields a finite dimensional system of ordinary differential equations (ODE)s of Schrödinger type, which can be solved numerically by exponential time integrators of Magnus type. The efficiency may be enhanced by combining these with operator splittings. We will discuss several different approaches of employing exponential-based methods in conjunction with an adaptive Lanczos method for the evaluation of matrix exponentials and compare their accuracy and efficiency. For each integrator, we use defect-based local error estimators to enable adaptive time-stepping. This serves to reliably control the approximation error and reduce the computational effort

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Time adaptive Zassenhaus splittings for the Schrödinger equation in the semiclassical regime

Time dependent Schrödinger equations with conservative force field U commonly constitute a major challenge in the numerical approximation, especially when they are analysed in the semiclassical regime. Extremely high oscillations originate from the semiclassical parameter, and call for appropriate methods. We propose to employ a combination of asymptotic Zassenhaus splitting with time adaptivity. While the former turns the disadvantage of the semiclassical parameter into an advantage, leading to highly efficient methods with low error constants, the latter enables to choose an optimal time step and to speed up the calculations when the oscillations subside. We support the results with numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Computable upper error bounds for Krylov approximations to matrix exponentials and associated $φ$-functions

An a posteriori estimate for the error of a standard Krylov approximation to the matrix exponential is derived. The estimate is based on the defect (residual) of the Krylov approximation and is proven to constitute a rigorous upper bound on the error, in contrast to existing asymptotical approximations. It can be computed economically in the underlying Krylov space. In view of time-stepping applications, assuming that the given matrix is scaled by a time step, it is shown that the bound is asymptotically correct (with an order related to the dimension of the Krylov space) for the time step tending to zero. This means that the deviation of the error estimate from the true error tends to zero faster than the error itself. Furthermore, this result is extended to Krylov approximations of $φ$-functions and to improved versions of such approximations. The accuracy of the derived bounds is demonstrated by examples and compared with different variants known from the literature, which are also investigated more closely. Alternative error bounds are tested on examples, in particular a version based on the concept of effective order. For the case where the matrix exponential is used in time integration algorithms, a step size selection strategy is proposed and illustrated by experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Adaptive high-order splitting methods for systems of nonlinear evolution equations with periodic boundary conditions

We assess the applicability and efficiency of time-adaptive high-order splitting methods applied for the numerical solution of (systems of) nonlinear parabolic problems under periodic boundary conditions. We discuss in particular several applications generating intricate patterns and displaying nonsmooth solution dynamics. First we give a general error analysis for splitting methods for parabolic problems under periodic boundary conditions and derive the necessary smoothness requirements on the exact solution in particular for the Gray-Scott equation and the Van der Pol equation. Numerical examples demonstrate the convergence of the methods and serve to compare the efficiency of different time-adaptive splitting schemes and of splitting into either two or three operators, based on appropriately constructed a posteriori local error estimators.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Adaptive splitting methods for nonlinear Schrödinger equations in the semiclassical regime

The error behavior of exponential operator splitting methods for nonlinear Schr{ö}dinger equations in the semiclassical regime is studied. For the Lie and Strang splitting methods, the exact form of the local error is determined and the dependence on the semiclassical parameter is identified. This is enabled within a defect-based framework which also suggests asymptotically correct a~posteriori local error estimators as the basis for adaptive time stepsize selection. Numerical examples substantiate and complement the theoretical investigations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Convergence of a Strang splitting finite element discretization for the Schrödinger-Poisson equation

Operator splitting methods combined with finite element spatial discretizations are studied for time-dependent nonlinear Schrödinger equations. In particular, the Schrödinger-Poisson equation under homogeneous Dirichlet boundary conditions on a finite domain is considered. A rigorous stability and error analysis is carried out for the second-order Strang splitting method and conforming polynomial finite element discretizations. For sufficiently regular solutions the classical orders of convergence are retained, that is, second-order convergence in time and polynomial convergence in space is proven. The established convergence result is confirmed and complemented by numerical illustrations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Practical splitting methods for the adaptive integration of nonlinear evolution equations. Part I: Construction of optimized schemes and pairs of schemes

We present a number of new contributions to the topic of constructing efficient higher-order splitting methods for the numerical integration of evolution equations. Particular schemes are constructed via setup and solution of polynomial systems for the splitting coefficients. To this end we use and modify a recent approach for generating these systems for a large class of splittings. In particular, various types of pairs of schemes intended for use in adaptive integrators are constructed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Setup of Order Conditions for Splitting Methods

This article is based on earlier papers where an approach based on Taylor expansion and the structure of its leading term as an element of a free Lie algebra was described for the setup of a system of order conditions for operator splitting methods. Along with a brief review of these materials and some theoretical background, we discuss the implementation of the ideas from these papers in computer algebra, in particular using Maple 18. A parallel version of such a code is described.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Symbolic Manipulation of Flows of Nonlinear Evolution Equations, with Application in the Analysis of Split-Step Time Integrators

We describe a package realized in the Julia programming language which performs symbolic manipulations applied to nonlinear evolution equations, their flows, and commutators of such objects. This tool was employed to perform contrived computations arising in the analysis of the local error of operator splitting methods. It enabled the proof of the convergence of the basic method and of the asymptotical correctness of a defect-based error estimator. The performance of our package is illustrated on several examples.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Fast stray field computation on tensor grids

A direct integration algorithm is described to compute the magnetostatic field and energy for given magnetization distributions on not necessarily uniform tensor grids. We use an analytically-based tensor approximation approach for function-related tensors, which reduces calculations to multilinear algebra operations. The algorithm scales with N^(4/3) for N computational cells used and with N^(2/3) (sublinear) when magnetization is given in canonical tensor format. In the final section we confirm our theoretical results concerning computing times and accuracy by means of numerical examples.

0 citations0 reviewsNo code linkedOpen access