Local error estimation for multistep collocation methods

Multistep collocation methods for initial value problems in ordinary differential equations are known to be a subclass of multistep Runge-Kutta methods and a generalisation of the well-known class of one-step collocation methods as well as of the one-leg methods of Dahlquist. In this paper we derive an error estimation method of embedded type for multistep collocation methods based on perturbed multistep collocation methods. This parallels and generalizes the results for one-step collocation methods by Nørsett and Wanner. Simple numerical experiments show that this error estimator agrees well with a theoretical error estimate which is a generalisation of an error estimate first derived by Dahlquist for one-leg methods.     

Randomized search directions in descent methods for minimizing certain quasidifferentiable functions

Several descent methods have recently been proposed for minimizing smooth compositions of max-type functions. The methods generate many search directions at each iteration. This paper shows that a random choice of only two search directions at each iteration suffices for retaining convergence to inf-stationary points with probability 1. This technique may decrease significantly the work in quadratic programming and line searches, thus enabling efficient implementations of the methods.     

Duality estimates for the numerical solution of integral equations

Summary We formulate and prove Aubin-Nitsche-type duality estimates for the error of general projection methods. Examples of applications include collocation methods and augmented Galerkin methods for boundary integral equations on plane domains with corners and three-dimensional screen and crack problems. For some of these methods, we obtain higher order error estimates in negative norms in cases where previous formulations of the duality arguments were not applicable.     

Reducible quadrature methods for Volterra integral equations of the first kind

Quadrature rules, generated by linear multistep methods for ordinary differential equations, are employed to construct a wide class of direct quadrature methods for the numerical solution of first kind Volterra integral equations. Our class covers several methods previously considered in the literature. The methods are convergent provided that both the first and second characteristic polynomial of the linear multistep method satisfy the root condition. Furthermore, the stability behaviour for fixed positive values of the stepsizeh is analyzed, and it turns out that convergence implies (fixedh) stability. The subclass formed by the backward differentiation methods up to order six is discussed and illustrated with numerical examples.     

Runge-Kutta theory for Volterra integrodifferential equations

Summary The present paper develops the theory of general Runge-Kutta methods for Volterra integrodifferential equations. The local order is characterized in terms of the coefficients of the method. We investigate the global convergence of mixed and extended Runge-Kutta methods and give results on asymptotic error expansions. In a further section we construct examples of methods up to order 4.     

The stability function for multistep collocation methods

Summary C-polynomials for rational approximation to the exponential function was introduced by Nørsett [7] to study stability properties of one-step methods. For one-step collocation methods theC-polynomial has a very simple form. In this paper we studyC-polynomials for multistep collocation methods and obtain results that generalize those in the one-step case, and provide a way to analyze linear stability of such methods.     

Approximation of the solution of a fourth order boundary value problem with nonsmooth coefficient

Summary A class of generalized finite element methods for the approximate solution of fourth order two point boundary value problem with nonsmooth coefficient is presented. The methods are based on the use of problem dependentL-splines incorporating the nonsmoothness of the coefficient. Stability is proved and optimal error estimates in theH ² norm are derived for the solution and postprocessed solution, under the assumption that the coefficient is of bounded variation. The relation of these methods to mixed methods is discussed.This research was sponsored by the Senate Research Committee of Syracuse University, Syracuse, NY 13210     

Stability properties of variable stepsize variable formula methods

Summary When variable stepsize variable formula methods (VSVFM's) are used in the solution of systems of first order differential equations instability arises sometimes. Therefore it is important to find VSVFM's whose zerostability properties are not affected by the choice of both the stepsize and the formula. The Adams VSVFM's are such methods. In this work a more general class of methods which contains the Adams VSVFM's is discussed and it is proved that the zero-stability of the class is not affected by the choice of the stepsize and of the formula.     

Stochastic properties of quadrature formulas

Summary The average error of suitable quadrature formulas and the stochastic error of Monte Carlo methods are both much smaller than the worst case error in many cases. This depends, however, on the classF of functions which is considered and there are counterexamples as well.Nonlinear methods, adaptive methods, or even methods with varying cardinality are not significantly better (with respect to certain stochastic error bounds) than the simplest linear methods .     

On the order conditions of Runge-Kutta methods with higher derivatives

Summary Runge-Kutta methods have been generalized to procedures with higher derivatives of the right side ofy=f(t,y) e.g. by Fehlberg 1964 and Kastlunger and Wanner 1972. In the present work some sufficient conditions for the order of consistence are derived for these methods using partially the degree of the corresponding numerical integration formulas. In particular, methods of Gauß, Radau, and Lobatto type are generalized to methods with higher derivatives and their maximum order property is proved. The applied technique was developed by Crouzeix 1975 for classical Runge-Kutta methods. Examples of simple explicit and semi-implicit methods are given up to order 7 and 6 respectively.     

On the order of composite multistep methods for ordinary differential equations

Summary In this paper, a general class ofk-step methods for the numerical solution of ordinary differential equations is discussed. It is shown that methods with order of consistencyq have order of convergence (q+1) if a very simple condition is satisfied. This result gives a new aspect to previous results of Spijker; it also serves as a starting point for a new theory of cyclick-step methods, completing an approach of Donelson and Hansen. It facilitates the practical determination of high-order cyclick-step methods, especially of stiffly stable,k-step methods.     

On the order of iterated defect correction

Summary In a recent article [2] Frank and Überhuber define and motivate the method of iterated defect correction for Runge-Kutta methods. They prove a theorem on the order of that method using the theory of asymptotic expansions.In this paper we give similar results using the theory of Butcher series (see [4]). Our proofs are purely algebraic. We don't restrict our considerations to Runge-Kutta methods, but we admit arbitrary linear one-step methods. At the same time we consider more general defect functions as in [2].     

The application of Rosenbrock-Wanner type methods with stepsize control in differential-algebraic equations

Summary Two Rosenbrock-Wanner type methods for the numerical treatment of differential-algebraic equations are presented. Both methods possess a stepsize control and an index-1 monitor. The first method DAE34 is of order (3)4 and uses a full semi-implicit Rosenbrock-Wanner scheme. The second method RKF4DA is derived from the Runge-Kutta-Fehlberg 4(5)-pair, where a semi-implicit Rosenbrock-Wanner method is embedded, in order to solve the nonlinear equations. The performance of both methods is discussed in artificial test problems and in technical applications.     

Factorization iterative methods,M-operators andH-operators

Summary We set up here a general formalism for describing factorization iterative methods of the first order and we use it to review various methods that have been proposed in the literature; next we introduce the notions ofM- andH-operators which generalize those of block-M- and block-H-matrices; finally we discuss the properties of factorization iterative methods in relation with characteristic properties ofM- andH-operators.     

Steplength optimization and linear multigrid methods

Summary Recently steplength parameters have been used in linear multigrid methods. In this paper we give a theoretical analysis of the effects of steplength optimization in a rather general framework which covers two different implementations of steplength optimization in standard multigrid methods.     

Numerical treatment of delay differential equations by Hermite Interpolation

Summary A class of numerical methods for the treatment of delay differential equations is developed. These methods are based on the wellknown Runge-Kutta-Fehlberg methods. The retarded argument is approximated by an appropriate multipoint Hermite Interpolation. The inherent jump discontinuities in the various derivatives of the solution are considered automatically.Problems with piecewise continuous right-hand side and initial function are treated too. Real-life problems are used for the numerical test and a comparison with other methods published in literature.     

On implicit Runge-Kutta methods with higher derivatives

Two families of implicit Runge-Kutta methods with higher derivatives are (re-)considered generalizing classical Runge-Kutta methods of Butcher type and f Ehle type. For generalized Butcher methods the characteristic functionG() is represented by means of the node polynomial directly, thereby showing that in methods of maximum order,G() is connected withs-orthogonal polynomials in exactly the same way as Padé approximations in the classical case.     

Asymptotisches Verhalten und Konvergenzbeschleunigung von Iterationsfolgen

Summary The effectivity of iterative numerical methods depends on the rate of convergence. In this note general procedures to accelerate the convergence of finite-dimensional stationary one-step-methods (fixed point methods) by extrapolation methods are studied. In this connection the investigation of the asymptotic behaviour of the sequences is fundamental. Differentiability and contractivity qualities supposed in the following an asymptotic expansion for such iterative sequences is proved. Neglecting the remainder the expansion fulfils a linear difference equation with constant coefficients. Wynn's -algorithm work off this expansion term by term, and the attainable acceleration can be exactly estimated. Skelboe's convergence statement is refuted. First test results demonstrate the advantage of acceleration methods.

     

Asymptotic expansions of the global error of fixed-stepsize methods

Summary In his fundamental paper on general fixed-stepsize methods, Skeel [6] studied convergence properties, but left the existence of asymptotic expansions as an open problem. In this paper we give a complete answer to this question. For the special cases of one-step and linear multistep methods our proof is shorter than the published ones.Asymptotic expansions are the theoretical base for extrapolation methods.     

A class of direct methods for linear systems

Summary A class of methods of direct type for solving determined or underdetermined, full rank or deficient rank linear systems is presented and theoretically analyzed. The class can be considered as a generalization of the methods of Brent and Brown as restricted to linear systems and implicitly contains orthogonal,LU andLL ^T factorization methods.