Source author record

István Faragó

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

3works
3topics
4close collaborators

Published work

3 publications

preprint2020arXiv

The linearization methods as a basis to derive the relaxation and the shooting methods

This chapter investigates numerical solution of nonlinear two-point boundary value problems. It establishes a connection between three important, seemingly unrelated, classes of iterative methods, namely: the linearization methods, the relaxation methods (finite difference methods), and the shooting methods. It has recently been demonstrated that using finite differences to discretize the sequence of linear problems obtained by quasi-linearization, Picard linearization, or constant-slope linearization, leads to the usual iteration formula of the respective relaxation method. Thus, the linearization methods can be used as a basis to derive the relaxation methods. In this work we demonstrate that the shooting methods can be derived from the linearization methods, too. We show that relaxing a shooting trajectory, i.e. an initial value problem solution, is in fact a projection transformation. The obtained function, called projection trajectory, can be used to correct the initial condition. Using the new initial condition, we can find a new shooting trajectory, and so on. The described procedure is called shooting-projection iteration (SPI). We show that using the quasi-linearization equation to relax (project) the shooting trajectory leads to the usual shooting by Newton method, the constant-slope linearization leads to the usual shooting by constant-slope method, while the Picard linearization leads to the recently proposed shooting-projection method. Therefore, the latter method can rightfully be called shooting by Picard method. A possible application of the new theoretical results is suggested and numerical computer experiments are presented. MATLAB codes are provided.

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 {tn}n=0N\{t_n\}_{n=0}^N can be modeled as the image of an equidistant grid under a smooth deformation map, i.e., tn=Φ(τn)t_n = Φ(τ_n), where τn=n/Nτ_n = n/N and the map ΦΦ is monotonically increasing with Φ(0)=0Φ(0)=0 and Φ(1)=1Φ(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 NN. Given any strongly stable multistep method, there is an N∗N^* such that the method is zero stable for N>N∗N>N^*, provided that Φ∈C2[0,1]Φ\in C^2[0,1]. Thus zero stability holds on all nonuniform grids such that adjacent step sizes satisfy hn/hn−1=1+O(N−1)h_n/h_{n-1} = 1 + \mathrm O(N^{-1}) as N→∞N\rightarrow\infty. The results are exemplified for BDF-type methods.

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.