Font Size: a A A

Rigorous numerical analysis with high-order Taylor models

Posted on:2002-12-12Degree:Ph.DType:Dissertation
University:Michigan State UniversityCandidate:Hoefkens, JensFull Text:PDF
GTID:1460390011990900Subject:Mathematics
Abstract/Summary:
Interval techniques have been utilized for rigorous numerical analysis since the 1960s. However, their use has been limited by the dimensionality curse and the dependency problem. The recently developed Taylor model approach alleviates these problems and allows the use of validated numerics in a wide range of applications. To broaden the applicability of the Taylor model method, we introduce new algorithms for the inversion of functional relations and the integration of differential algebraic equations.; First we present a new method for computing verified enclosures of the inverses of given functions over large domains. The approach utilizes Taylor models and the sharpness of the enclosures scales with a high order of the domain. An integral part of the new method is the rigorous determination of invertibility of high dimensional functions over large domains, which is reduced to a verified linear algebra problem involving only first derivatives of the function of interest. Several examples highlighting various aspects of the methods are discussed.; Differential algebraic equations (DAEs) describe important problems in mechanical and chemical engineering. Existing algorithms for the integration of DAE initial value problems have traditionally been restricted to low-index systems and until recently no practical scheme for the verified integration of DAEs existed. Recognizing the antiderivation as a natural operation on Taylor models yields a method that treats DAEs within a fully differential algebraic context as implicit equations made of conventional functions and the antiderivation. The resulting integration scheme can be applied to high-index problems and allows the computation of guaranteed enclosures of final coordinates from large initial regions.; To demonstrate the general applicability of the Taylor model approach, we present results from verified asteroid orbit integrations and the theory of Hamiltonian systems. We show that the newly developed methods are practical and can indeed outperform conventional interval methods in a wide class of problems. Finally, we discuss some details of implementing interval libraries on general purpose computers and present the concept of language independent software development, which has been used for the design and implementation of the C++ and Fortran 90 interfaces to COSY Infinity.
Keywords/Search Tags:Taylor model, Rigorous
Related items