Results 1  10
of
81
Means and Averaging in the Group of Rotations
, 2002
"... In this paper we give precise definitions of different, properly invariant notions of mean or average rotation. Each mean is associated with a metric in SO(3). The metric induced from the Frobenius inner product gives rise to a mean rotation that is given by the closest special orthogonal matrix to ..."
Abstract

Cited by 63 (1 self)
 Add to MetaCart
In this paper we give precise definitions of different, properly invariant notions of mean or average rotation. Each mean is associated with a metric in SO(3). The metric induced from the Frobenius inner product gives rise to a mean rotation that is given by the closest special orthogonal matrix to the usual arithmetic mean of the given rotation matrices. The mean rotation associated with the intrinsic metric on SO(3) is the Riemannian center of mass of the given rotation matrices. We show that the Riemannian mean rotation shares many common features with the geometric mean of positive numbers and the geometric mean of positive Hermitian operators. We give some examples with closedform solutions of both notions of mean.
A Lie group variational integrator for the attitude dynamics of a rigid body with applications to the 3D pendulum
 Proceedings of the IEEE Conference on Control Applications
, 2005
"... Abstract — A numerical integrator is derived for a class of models that describe the attitude dynamics of a rigid body in the presence of an attitude dependent potential. The numerical integrator is obtained from a discrete variational principle, and exhibits excellent geometric conservation propert ..."
Abstract

Cited by 26 (19 self)
 Add to MetaCart
Abstract — A numerical integrator is derived for a class of models that describe the attitude dynamics of a rigid body in the presence of an attitude dependent potential. The numerical integrator is obtained from a discrete variational principle, and exhibits excellent geometric conservation properties. In particular, by performing computations at the level of the Lie algebra, and updating the solution using the matrix exponential, the attitude automatically evolves on the rotation group embedded in the space of matrices. The geometric conservation properties of the numerical integrator imply long time numerical stability. We apply this variational integrator to the uncontrolled 3D pendulum, that is a rigid asymmetric body supported at a frictionless pivot acting under the influence of uniform gravity. Interesting dynamics of the 3D pendulum are exposed. I.
On Magnus Integrators for TimeDependent Schrödinger Equations
 SIAM J. Numer. Anal
, 2002
"... Numerical methods based on the Magnus expansion are an ecient class of integrators for Schrodinger equations with timedependent Hamiltonian. Though their derivation assumes an unreasonably small time step size as would be required for a standard explicit integrator, the methods perform well even f ..."
Abstract

Cited by 23 (2 self)
 Add to MetaCart
Numerical methods based on the Magnus expansion are an ecient class of integrators for Schrodinger equations with timedependent Hamiltonian. Though their derivation assumes an unreasonably small time step size as would be required for a standard explicit integrator, the methods perform well even for much larger step sizes. This favorable behavior is explained, and optimalorder error bounds are derived which require no or only mild restrictions of the step size. In contrast to standard integrators, the error does not depend on higher time derivatives of the solution, which is in general highly oscillatory.
On the Global Error of Discretization Methods for HighlyOscillatory Ordinary Differential Equations
, 2000
"... Commencing from a globalerror formula, originally due to Henrici, we investigate the accumulation of global error in the numerical solution of linear highlyoscillating systems of the form y 00 + g(t)y = 0, where g(t) t!1 \Gamma! 1. Using WKB analysis we derive an explicit form of the globalerror ..."
Abstract

Cited by 20 (5 self)
 Add to MetaCart
Commencing from a globalerror formula, originally due to Henrici, we investigate the accumulation of global error in the numerical solution of linear highlyoscillating systems of the form y 00 + g(t)y = 0, where g(t) t!1 \Gamma! 1. Using WKB analysis we derive an explicit form of the globalerror envelope for RungeKutta and Magnus methods. Our results are closely matched by numerical experiments. Motivated by the superior performance of Liegroup methods, we present a modification of the Magnus expansion which displays even better longterm behaviour in the presence of oscillations.
Nonsmooth Lagrangian mechanics and variational collision integrators
 SIAM Journal on Applied Dynamical Systems
, 2003
"... Abstract. Variational techniques are used to analyze the problem of rigidbody dynamics with impacts. The theory of smooth Lagrangian mechanics is extended to a nonsmooth context appropriate for collisions, and it is shown in what sense the system is symplectic and satisfies a Noetherstyle momentum ..."
Abstract

Cited by 19 (7 self)
 Add to MetaCart
Abstract. Variational techniques are used to analyze the problem of rigidbody dynamics with impacts. The theory of smooth Lagrangian mechanics is extended to a nonsmooth context appropriate for collisions, and it is shown in what sense the system is symplectic and satisfies a Noetherstyle momentum conservation theorem. Discretizations of this nonsmooth mechanics are developed by using the methodology of variational discrete mechanics. This leads to variational integrators which are symplecticmomentum preserving and are consistent with the jump conditions given in the continuous theory. Specific examples of these methods are tested numerically, and the longtime stable energy behavior typical of variational methods is demonstrated.
Newton’s method on Riemannian manifolds: convariant alpha theory
 IMA J. Numer. Anal
, 2003
"... In this paper, Smale’s α theory is generalized to the context of intrinsic Newton iteration on geodesically complete analytic Riemannian and Hermitian manifolds. Results are valid for analytic mappings from a manifold to a linear space of the same dimension, or for analytic vector fields on the mani ..."
Abstract

Cited by 18 (2 self)
 Add to MetaCart
In this paper, Smale’s α theory is generalized to the context of intrinsic Newton iteration on geodesically complete analytic Riemannian and Hermitian manifolds. Results are valid for analytic mappings from a manifold to a linear space of the same dimension, or for analytic vector fields on the manifold. The invariant γ is defined by means of high order covariant derivatives. Bounds on the size of the basin of quadratic convergence are given. If the ambient manifold has negative sectional curvature, those bounds depend on the curvature. A criterion of quadratic convergence for Newton iteration from the information available at a point is also given. 1 Introduction and main results. Numerical problems posed in manifolds arise in many natural contexts. Classical examples are given by the eigenvalue problem, the symmetric eigenvalue problem, invariant subspace computations, minimization problems with orthogonality constraints, optimization problems with equality constraints... etc. In the first
Hamilton–Pontryagin integrators on Lie groups
, 2007
"... Abstract In this paper, structurepreserving timeintegrators for rigid bodytype mechanical ..."
Abstract

Cited by 17 (6 self)
 Add to MetaCart
Abstract In this paper, structurepreserving timeintegrators for rigid bodytype mechanical
Generalized Polar Decompositions for the Approximation of the Matrix Exponential
, 2000
"... In this paper we describe the use of the theory of generalized polar decompositions (MuntheKaas, Quispel & Zanna 2000b) to approximate a matrix exponential. The algorithms presented in this paper have the property that, if Z 2 g, a Lie algebra of matrices, then the approximation for exp(Z) resides i ..."
Abstract

Cited by 14 (5 self)
 Add to MetaCart
In this paper we describe the use of the theory of generalized polar decompositions (MuntheKaas, Quispel & Zanna 2000b) to approximate a matrix exponential. The algorithms presented in this paper have the property that, if Z 2 g, a Lie algebra of matrices, then the approximation for exp(Z) resides in G, the matrix Lie group of g. This property is very relevant when solving Liegroup ODEs and is not usually fulfilled by standard approximations to the matrix exponential. We propose algorithms based on a splitting of Z into matrices having a very simple structure, usually one row and one column (or a few rows and a few columns), whose exponential is computed very cheaply to machine accuracy. The proposed methods have a complexity of O \Gamma n 3 \Delta , the constant is small, depending on the order and the Lie algebra g. The algorithms are recommended in cases where it is of fundamental importance that the approximation for the exponential resides in G, and when the order of approx...
Geometric Integration and Its Applications
 in Handbook of numerical analysis
, 2000
"... This paper aims to give an introduction to the relatively new eld of geometric integration. ..."
Abstract

Cited by 12 (1 self)
 Add to MetaCart
This paper aims to give an introduction to the relatively new eld of geometric integration.