A posteriori estimates of inverse operators for initial value problems in linear ordinary differential equations

We present constructive a posteriori estimates of inverse operators for initial value problems in linear ordinary differential equations (ODEs) on a bounded interval. Here, “constructive” indicates that we can obtain bounds of the operator norm in which all constants are explicitly given or are represented in a numerically computable form. In general, it is difficult to estimate these inverse operators a priori. We, therefore, propose a technique for obtaining a posteriori estimates by using Galerkin approximation of inverse operators. This type of estimation will play an important role in the numerical verification of solutions for initial value problems in nonlinear ODEs as well as for parabolic initial boundary value problems.     

Boundary value methods for the solution of differential-algebraic equations

Summary Boundary value techniques for the solution of initial value problems of ODEs, despite their apparent higher cost, present some important advantages over initial value methods. Among them, there is the possibility to have greater accuracy, to control the global error, and to have an efficient parallel implementation.In this paper, the same techniques are applied to the solution of linear initial value problems of DAEs. We have considered three term numerical methods (Midpoint, Simpson, and an Adams type method) in order to obtain a block tridiagonal linear system as a discrete problem.Convergence results are stated in the case of constant coefficients, and numerical examples are given on linear time-varying problems.Work supported by the Ministero della Ricerca Scientifica, 40% project, and by the C.N.R. (contract of research # 92.00535.CT01)     

On dynamic iteration for delay differential equations

In this paper we study dynamic iteration techniques for systems of nonlinear delay differential equations. After pointing out a close connection to the truncated infinite embedding, as proposed by Feldstein, Iserles, and Levin, we give a proof of the superlinear convergence of the simple dynamic iteration scheme. Then we propose a more general scheme that in addition allows for a decoupling of the equations into disjoint subsystems, just like what we are used to from dynamic iteration schemes for ODEs. This scheme is also shown to converge superlinearly.     

Defect-based local error estimators for splitting methods,with application to Schrödinger equations,Part I: The linear case

We introduce a defect correction principle for exponential operator splitting methods applied to time-dependent linear Schrödinger equations and construct a posteriori local error estimators for the Lie–Trotter and Strang splitting methods. Under natural commutator bounds on the involved operators we prove asymptotical correctness of the local error estimators, and along the way recover the known a priori convergence bounds. Numerical examples illustrate the theoretical local and global error estimates.     

Optimal order results for a class of regularization methods using unbounded operators

A class of regularization methods using unbounded regularizing operators is considered for obtaining stable approximate solutions for ill-posed operator equations. With an a posteriori as well as an a priori parameter choice strategy, it is shown that the method yields the optimal order. Error estimates have also been obtained under stronger assumptions on the generalized solution. The results of the paper unify and simplify many of the results available in the literature. For example, the optimal results of the paper include, as particular cases for Tikhonov regularization, the main result of Mair (1994) with an a priori parameter choice, and a result of Nair (1999) with an a posteriori parameter choice. Thus the observations of Mair (1994) on Tikhonov regularization of ill-posed problems involving finitely and infinitely smoothing operators is applicable to various other regularization procedures as well. Subsequent results on error estimates include, as special cases, an optimal result of Vainikko (1987) and also some recent results of Tautenhahn (1996) in the setting of Hilbert scales.     

Convergence results for the MATLAB ode23 routine

We prove convergence results on finite time intervals, as the user-defined tolerance τ→0, for a class of adaptive timestepping ODE solvers that includes the ode23 routine supplied in MATLAB Version 4.2. In contrast to existing theories, these convergence results hold with error constants that are uniform in the neighbourhood of equilibria; such uniformity is crucial for the derivation of results concerning the numerical approximation of dynamical systems. For linear problems the error estimates are uniform on compact sets of initial data. The analysis relies upon the identification of explicit embedded Runge-Kutta pairs for which all but the leading order terms of the expansion of the local error estimate areO(∥f(u∥)²). This work was partially supported by NSF Grant DMS-95-04879.     

A Test Set for stiff Initial Value Problem Solvers in the open source software R: Package deTestSet

In this paper we present the R package deTestSet that includes challenging test problems written as ordinary differential equations (ODEs), differential algebraic equations (DAEs) of index up to 3 and implicit differential equations (IDEs). In addition it includes 6 new codes to solve initial value problems (IVPs). The R package is derived from the Test Set for Initial Value Problem Solvers available at http://www.dm.uniba.it/~testset which includes documentation of the test problems, experimental results from a number of proven solvers, and Fortran subroutines providing a common interface to the defining problem functions. Many of these facilities are now available in the R package deTestSet, which comprises an R interface to the test problems and to most of the Fortran solvers. The package deTestSet is free software which is distributed under the GNU General Public License, as part of the R open source software project.     

Stability of some boundary value methods for the solution of initial value problems

The stability properties of three particular boundary value methods (BVMs) for the solution of initial value problems are considered. Our attention is focused on the BVMs based on the midpoint rule, on the Simpson method and on an Adams method of order 3. We investigate their BV-stability regions by considering the scalar test problem and constant stepsize. The study of the conditioning of the coefficient matrix of the discrete problem is extended to the case of variable stepsize and block ODE problems. We also analyse an appropriate choice for the stepsize for stiff problems. Numerical tests are reported to evidentiate the effectiveness of the BVMs and the differences among the BVMs considered.Work supported by the Ministero della Ricerca Scientifica, 40% project, and C.N.R. (contract of research # 92.00535.01).     

Convergence of linear multistep and one-leg methods for stiff nonlinear initial value problems

To prove convergence of numerical methods for stiff initial value problems, stability is needed but also estimates for the local errors which are not affected by stiffness. In this paper global error bounds are derived for one-leg and linear multistep methods applied to classes of arbitrarily stiff, nonlinear initial value problems. It will be shown that under suitable stability assumptions the multistep methods are convergent for stiff problems with the same order of convergence as for nonstiff problems, provided that the stepsize variation is sufficiently regular.     

Adaptive boundary element methods for strongly elliptic integral equations

Summary Integral operators are nonlocal operators. The operators defined in boundary integral equations to elliptic boundary value problems, however, are pseudo-differential operators on the boundary and, therefore, provide additional pseudolocal properties. These allow the successful application of adaptive procedures to some boundary element methods. In this paper we analyze these methods for general strongly elliptic integral equations and obtain a-posteriori error estimates for boundary element solutions. We also apply these methods to nodal collocation with odd degree splines. Some numerical examples show that these adaptive procedures are reliable and effective.This work was carried out while Dr. De-hao Yu was an Alexander-von-Humboldt-Stiftung research fellow at the University of Stuttgart in 1987, 1988     

Harmonic analysis for star graphs and the spherical coordinate trapezoidal rule

Novel ideas in harmonic analysis are used to analyze the trapezoidal rule integration for two spheres. Sampling in spherical coordinates links three levels of harmonic analysis. Eigenfunctions of a nonstandard manifold Laplacian descend by restriction, first to a differential graph Laplacian, and then to difference operators. Trapezoidal rule integration with appropriate sampling is exact on eigenspaces of the manifold Laplacian, a fact which leads to trapezoidal rule error estimates on Sobolev-style spaces of functions. Singular functions with accurate trapezoidal rule integrals are identified, and a simplified analysis of smooth function numerical integration is provided.     

Design of software for ODEs

We discuss the design of software that is easy to use for simple problems, but still capable of solving complicated problems. Nowadays a design should accommodate related, but fundamentally different, tasks such as solving DAEs and DDEs in addition to IVPs and BVPs for ODEs. Major issues are illustrated with the MATLAB ODE Suite and selected solvers written in Fortran 90 and Maple.     

Sturm-Liouville systems with rational Weyl functions: Explicit formulas and applications

This paper solves explicitly the direct and inverse spectral problems for Sturm-Liouville systems with rational Weyl functions. As in the previous papers of the authors on canonical and pseudocanonical systems, the method employed in the present paper is based on state space techniques and uses the idea of realization from mathematical system theory. For a subclass of rational potentials a property of modified bispectrality is proven. A wide class of solutions of the matrix Korteweg-de Vries equation is derived.To the memory of M.G. Kreîn, one of the founding fathers of modern spectral theory of differential operators, with admiration and gratitude.     

On the a posteriori estimates for inverse operators of linear parabolic equations with applications to the numerical enclosure of solutions for nonlinear problems

We consider the guaranteed a posteriori estimates for the inverse parabolic operators with homogeneous initial-boundary conditions. Our estimation technique uses a full-discrete numerical scheme, which is based on the Galerkin method with an interpolation in time by using the fundamental solution for semidiscretization in space. In our technique, the constructive a priori error estimates for a full discretization of solutions for the heat equation play an essential role. Combining these estimates with an argument for the discretized inverse operator and a contraction property of the Newton-type formulation, we derive an a posteriori estimate of the norm for the infinite-dimensional operator. In numerical examples, we show that the proposed method should be more efficient than the existing method. Moreover, as an application, we give some prototype results for numerical verification of solutions of nonlinear parabolic problems, which confirm the actual usefulness of our technique.     

Backward Euler method for abstract time-dependent parabolic equations with variable domains

We consider semidiscretizations in time, based on the backward Euler method, of an abstract, non-autonomous parabolic initial value problem where , , is a family of sectorial operators in a Banach space X. The domains are allowed to depend on t. Our hypotheses are fulfilled for classical parabolic problems in the , , norms. We prove that the semidiscretization is stable in a suitable sense. We get optimal estimates for the error even when non-homogeneous boundary values are considered. In particular, the results are applicable to the analysis of the semidiscretizations of time-dependent parabolic problems under non-homogeneous Neumann boundary conditions. Received October 17, 1997 / Revised version received April 17, 1998     

Locally extra-optimal regularizing algorithms and a posteriori estimates of the accuracy for ill-posed problems with discontinuous solutions

Local a posteriori estimates of the accuracy of approximate solutions to ill-posed inverse problems with discontinuous solutions from the classes of functions of several variables with bounded variations of the Hardy or Giusti type are studied. Unlike global estimates (in the norm), local estimates of accuracy are carried out using certain linear estimation functionals (e.g., using the mean value of the solution on a given fragment of its support). The concept of a locally extra-optimal regularizing algorithm for solving ill-posed inverse problems, which has an optimal in order local a posteriori estimate, was introduced. A method for calculating local a posteriori estimates of accuracy with the use of some distinguished classes of linear functionals for the problems with discontinuous solutions is proposed. For linear inverse problems, the method is bases on solving specialized convex optimization problems. Examples of locally extra-optimal regularizing algorithms and results of numerical experiments on a posteriori estimation of the accuracy of solutions for different linear estimation functionals are presented.     

Asymptotically correct error estimation for collocation methods applied to singular boundary value problems

We discuss an a posteriori error estimate for collocation methods applied to boundary value problems in ordinary differential equations with a singularity of the first kind. As an extension of previous results we show the asymptotical correctness of our error estimate for the most general class of singular problems where the coefficient matrix is allowed to have eigenvalues with positive real parts. This requires a new representation of the global error for the numerical solution obtained by piecewise polynomial collocation when applied to our problem class.     

Einschließungsaussagen bei Differentialoperatoren zweiter Ordnung durch punktweise Ungleichungen

Summary We consider non-inverse positive differential operators of the second order and show that error estimates can be obtained through methods of inverse positivity alone. This means a significant simplification if the operator has only one negative eigenvalue. In addition we show how to obtain Range-Domain implications for operators with mixed boundary conditions.

The author would like to thank the European Research Office of the United States Army for their kind interest     

A posteriori error estimates for mixed finite element solutions of convex optimal control problems

In this paper, we present an a posteriori error analysis for mixed finite element approximation of convex optimal control problems. We derive a posteriori error estimates for the coupled state and control approximations under some assumptions which hold in many applications. Such estimates can be used to construct reliable adaptive mixed finite elements for the control problems.     

Basic theory of fractional differential equations

In this paper, the basic theory for the initial value problem of fractional differential equations involving Riemann–Liouville differential operators is discussed employing the classical approach. The theory of inequalities, local existence, extremal solutions, comparison result and global existence of solutions are considered.