Source author record

Imre Fekete

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

4works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

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

4 published item(s)

preprint2022arXiv

Learning the parameters of a differential equation from its trajectory via the adjoint equation

The paper contributes to strengthening the relation between machine learning and the theory of differential equations. In this context, the inverse problem of fitting the parameters, and the initial condition of a differential equation to some measurements constitutes a key issue. The paper explores an abstraction that can be used to construct a family of loss functions with the aim of fitting the solution of an initial value problem to a set of discrete or continuous measurements. It is shown, that an extension of the adjoint equation can be used to derive the gradient of the loss function as a continuous analogue of backpropagation in machine learning. Numerical evidence is presented that under reasonably controlled circumstances the gradients obtained this way can be used in a gradient descent to fit the solution of an initial value problem to a set of continuous noisy measurements, and a set of discrete noisy measurements that are recorded at uncertain times.

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.

preprint2018arXiv

Embedded error estimation and adaptive step-size control for optimal explicit strong stability preserving Runge--Kutta methods

We construct a family of embedded pairs for optimal strong stability preserving explicit Runge-Kutta methods of order $2 \leq p \leq 4$ to be used to obtain numerical solution of spatially discretized hyperbolic PDEs. In this construction, the goals include non-defective methods, large region of absolute stability, and optimal error measurement as defined in [5,19]. The new family of embedded pairs offer the ability for strong stability preserving (SSP) methods to adapt by varying the step-size based on the local error estimation while maintaining their inherent nonlinear stability properties. Through several numerical experiments, we assess the overall effectiveness in terms of precision versus work while also taking into consideration accuracy and stability.

preprint2018arXiv

On the zero-stability of multistep methods on smooth nonuniform grids

In order to be convergent, linear multistep methods must be zero stable. While constant step size theory was established in the 1950's, zero stability on nonuniform grids is less well understood. Here we investigate zero stability on compact intervals and smooth nonuniform grids. In practical computations, step size control can be implemented using smooth (small) step size changes. The resulting grid $\{t_n\}_{n=0}^N$ can be modeled as the image of an equidistant grid under a smooth deformation map, i.e., $t_n = Φ(τ_n)$, where $τ_n = n/N$ and the map $Φ$ is monotonically increasing with $Φ(0)=0$ and $Φ(1)=1$. The model is justified for any fixed order method operating in its asymptotic regime when applied to smooth problems, since the step size is then determined by the (smooth) principal error function which determines $Φ$, and a tolerance requirement which determines $N$. Given any strongly stable multistep method, there is an $N^*$ such that the method is zero stable for $N>N^*$, provided that $Φ\in C^2[0,1]$. Thus zero stability holds on all nonuniform grids such that adjacent step sizes satisfy $h_n/h_{n-1} = 1 + \mathrm O(N^{-1})$ as $N\rightarrow\infty$. The results are exemplified for BDF-type methods.