Source author record

Lajos Lóczi

Lajos Lóczi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA math.DS Numerical Analysis

Catalog footprint

What is connected

6works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Linear multistep methods and global Richardson extrapolation

In this work, we study the application the classical Richardson extrapolation (RE) technique to accelerate the convergence of sequences resulting from linear multistep methods (LMMs) for solving initial-value problems of systems of ordinary differential equations numerically. The advantage of the LMM-RE approach is that the combined method possesses higher order and favorable linear stability properties in terms of $A$- or $A(α)$-stability, and existing LMM codes can be used without any modification.

preprint2016arXiv

Dense output for strong stability preserving Runge-Kutta methods

We investigate dense output formulae (also known as continuous extensions) for strong stability preserving (SSP) Runge-Kutta methods. We require that the dense output formula also possess the SSP property, ideally under the same step-size restriction as the method itself. A general recipe for first-order SSP dense output formulae for SSP methods is given, and second-order dense output formulae for several optimal SSP methods are developed. It is shown that SSP dense output formulae of order 3 and higher do not exist, and that in any method possessing a second-order SSP dense output, the coefficient matrix A has a zero row.

preprint2016arXiv

Strong stability preserving explicit linear multistep methods with variable step size

Strong stability preserving (SSP) methods are designed primarily for time integration of nonlinear hyperbolic PDEs, for which the permissible SSP step size varies from one step to the next. We develop the first SSP linear multistep methods (of order two and three) with variable step size, and prove their optimality, stability, and convergence. The choice of step size for multistep SSP methods is an interesting problem because the allowable step size depends on the SSP coefficient, which in turn depends on the chosen step sizes. The description of the methods includes an optimal step-size strategy. We prove sharp upper bounds on the allowable step size for explicit SSP linear multistep methods and show the existence of methods with arbitrarily high order of accuracy. The effectiveness of the methods is demonstrated through numerical examples.

preprint2014arXiv

Discretizing the transcritical and pitchfork bifurcations -- conjugacy results

We present two case studies in one-dimensional dynamics concerning the discretization of transcritical (TC) and pitchfork (PF) bifurcations. In the vicinity of a TC or PF bifurcation point and under some natural assumptions on the one-step discretization method of order $p\ge 1$, we show that the time-$h$ exact and the step-size-$h$ discretized dynamics are topologically equivalent by constructing a two-parameter family of conjugacies in each case. As a main result, we prove that the constructed conjugacy maps are ${\cal{O}}(h^p)$-close to the identity and these estimates are optimal.

preprint2014arXiv

Propagation of internal errors in explicit Runge--Kutta methods and internal stability of SSP and extrapolation methods

In practical computation with Runge--Kutta methods, the stage equations are not satisfied exactly, due to roundoff errors, algebraic solver errors, and so forth. We show by example that propagation of such errors within a single step can have catastrophic effects for otherwise practical and well-known methods. We perform a general analysis of internal error propagation, emphasizing that it depends significantly on how the method is implemented. We show that for a fixed method, essentially any set of internal stability polynomials can be obtained by modifying the implementation details. We provide bounds on the internal error amplification constants for some classes of methods with many stages, including strong stability preserving methods and extrapolation methods. These results are used to prove error bounds in the presence of roundoff or other internal errors.

preprint2013arXiv

On the absolute stability regions corresponding to partial sums of the exponential function

Certain numerical methods for initial value problems have as stability function the nth partial sum of the exponential function. We study the stability region, i.e., the set in the complex plane over which the nth partial sum has at most unit modulus. It is known that the asymptotic shape of the part of the stability region in the left half-plane is a semi-disk. We quantify this by providing disks that enclose or are enclosed by the stability region or its left half-plane part. The radius of the smallest disk centered at the origin that contains the stability region (or its portion in the left half-plane) is determined for $1\le n\le 20$. Bounds on such radii are proved for $n\ge 2$; these bounds are shown to be optimal in the limit $n\to +\infty$. We prove that the stability region and its complement, restricted to the imaginary axis, consist of alternating intervals of length tending to $π$, as $n\to\infty$. Finally, we prove that a semi-disk in the left half-plane with vertical boundary being the imaginary axis and centered at the origin is included in the stability region if and only if $n\equiv 0 \mod 4$ or $n\equiv 3\mod 4$. The maximal radii of such semi-disks are exactly determined for $1\le n\le 20$.