Division Operation Simulation

Práca sa zaoberá numerickou integráciou a operáciou delenia. Najskôr je čitateľ oboznámený s numerickým riešením diferenciálnych rovníc s operáciou delenia pomocou Taylorovej rady. Ďalej je vysvetlený princíp delenia v hardvéri algoritmom SRT a je predstavený návrh sériovo-paralelného a paralelného deliaceho integrátora v pevnej rádovej čiarke. Praktickým cieľom práce je implementácia paralelného deliaceho integrátora a vytvorenie programového simulátora tohoto integrátora.This work deals with numerical integration and division operation. The reader is acquainted with the numerical solution of differential equations using division by the Taylor series. Next is explained the principle of SRT division in hardware and introduction of draft of design series-parallel and parallel division integrator in fixed point arithmetic. The practical aim of this work is implementation parallel division integrator and development of a software simulation of this integrator.

Multiple Operation Simulation

Tato práce se zabývá numerickou integrací. Nejprve je popsáno numerické řešení diferenciálních rovnic použitím metody Taylorovy řady. Poté jsou popsány jednotlivé varianty integrátorů. V praktické části je popsán návrh dvouvstupého integrátoru násobení a dále jeho realizace v prostředí FPGA. Pro tento integrátor je také vytvořen simulátor znázorňující jeho funkci.This work deals with the numeric integration. The reader is acquainted with the numerical solution of differential equations using Taylor series method. Then describes the different variants integrators. The practical part describes the design double-input integrator with multiplication and its implementation in an FPGA. For this integrator is also provided a simulator demonstrate its function.

Shuttle program. Onorbit navigation integrator results for typical shuttle orbits

Three types of navigation onorbit numerical integrators were evaluated: (1) power integrators with no delta-V incorporation, just coasting (using Taylor series expansion integrators); (2) coasting integrators using the Cowell method of special perturbations; and (3) coasting integrator using the Pines variation of parameter perturbation method. Results show that the super G integrator is a very simple and effective for 2 and 4 second time steps. Since IMU delta-V data can be easily incorporated in the integration scheme, its use as the standard onorbit navigation propagator for the maintenance of the current state was implemented in the onboard navigation software. The Pines formulation method with a Runge-Kutta-Gill fourth-order integrator method produces excellent results up to 300 second time steps. On orbit prediction with this method was implemented in the onboard onorbit navigation scheme. The Runge-Kutta third order, using Cowell's method, is an excellent general purpose determination integrator for time steps up to a 60 second duration

Algebraic Structures and Stochastic Differential Equations driven by Levy processes

We construct an efficient integrator for stochastic differential systems driven by Levy processes. An efficient integrator is a strong approximation that is more accurate than the corresponding stochastic Taylor approximation, to all orders and independent of the governing vector fields. This holds provided the driving processes possess moments of all orders and the vector fields are sufficiently smooth. Moreover the efficient integrator in question is optimal within a broad class of perturbations for half-integer global root mean-square orders of convergence. We obtain these results using the quasi-shuffle algebra of multiple iterated integrals of independent Levy processes.Comment: 41 pages, 11 figure

Some self starting integrators for x Prime equals f (x, t)

The integration is discussed of the vector differential equation X = F(x, t) from time t sub i to t sub (i = 1) where only the values of x sub i are available for the the integration. No previous values of x or x prime are used. Using an orbit integration problem, comparisons are made between Taylor series integrators and various types and orders of Runge-Kutta integrators. A fourth order Runge-Kutta type integrator for orbital work is presented, and approximate (there may be no exact) fifth order Runge-Kutta integrators are discussed. Also discussed and compared is a self starting integrator ising delta f/delta x. A numerical method for controlling the accuracy of integration is given, and the special equations for accurately integrating accelerometer data are shown

Algebraic structure of stochastic expansions and efficient simulation

We investigate the algebraic structure underlying the stochastic Taylor solution expansion for stochastic differential systems.Our motivation is to construct efficient integrators. These are approximations that generate strong numerical integration schemes that are more accurate than the corresponding stochastic Taylor approximation, independent of the governing vector fields and to all orders. The sinhlog integrator introduced by Malham & Wiese (2009) is one example. Herein we: show that the natural context to study stochastic integrators and their properties is the convolution shuffle algebra of endomorphisms; establish a new whole class of efficient integrators; and then prove that, within this class, the sinhlog integrator generates the optimal efficient stochastic integrator at all orders.Comment: 19 page

Efficient integration of the variational equations of multi-dimensional Hamiltonian systems: Application to the Fermi-Pasta-Ulam lattice

We study the problem of efficient integration of variational equations in multi-dimensional Hamiltonian systems. For this purpose, we consider a Runge-Kutta-type integrator, a Taylor series expansion method and the so-called `Tangent Map' (TM) technique based on symplectic integration schemes, and apply them to the Fermi-Pasta-Ulam $\beta$ (FPU-$\beta$) lattice of $N$ nonlinearly coupled oscillators, with $N$ ranging from 4 to 20. The fast and accurate reproduction of well-known behaviors of the Generalized Alignment Index (GALI) chaos detection technique is used as an indicator for the efficiency of the tested integration schemes. Implementing the TM technique--which shows the best performance among the tested algorithms--and exploiting the advantages of the GALI method, we successfully trace the location of low-dimensional tori.Comment: 14 pages, 6 figure

B-series methods are exactly the affine equivariant methods

Butcher series, also called B-series, are a type of expansion, fundamental in the analysis of numerical integration. Numerical methods that can be expanded in B-series are defined in all dimensions, so they correspond to \emph{sequences of maps}---one map for each dimension. A long-standing problem has been to characterise those sequences of maps that arise from B-series. This problem is solved here: we prove that a sequence of smooth maps between vector fields on affine spaces has a B-series expansion if and only if it is \emph{affine equivariant}, meaning it respects all affine maps between affine spaces