Solving the Discrete Euler-Arnold Equations for the Generalized Rigid Body Motion

by   João R. Cardoso, et al.

We propose three iterative methods for solving the Moser-Veselov equation, which arises in the discretization of the Euler-Arnold differential equations governing the motion of a generalized rigid body. We start by formulating the problem as an optimization problem with orthogonal constraints and proving that the objective function is convex. Then, using techniques from optimization on Riemannian manifolds, the three feasible algorithms are designed. The first one splits the orthogonal constraints using the Bregman method, whereas the other two methods are of the steepest-descent type. The second method uses the Cayley-transform to preserve the constraints and a Barzilai-Borwein step size, while the third one involves geodesics, with the step size computed by Armijo's rule. Finally, a set of numerical experiments are carried out to compare the performance of the proposed algorithms, suggesting that the first algorithm has the best performance in terms of accuracy and number of iterations. An essential advantage of these iterative methods is that they work even when the conditions for applicability of the direct methods available in the literature are not satisfied.


page 1

page 2

page 3

page 4


An Euler-type method for Volterra integro-differential equations

We describe an algorithm, based on Euler's method, for solving Volterra ...

Deep Euler method: solving ODEs by approximating the local truncation error of the Euler method

In this paper, we propose a deep learning-based method, deep Euler metho...

Optimal-order error estimates for Euler discretization of high-index saddle dynamics

High-index saddle dynamics provides an effective means to compute the an...

A boundary-penalized isogeometric analysis for second-order hyperbolic equations

Explicit time-marching schemes are popular for solving time-dependent pa...

A Bregman-Kaczmarz method for nonlinear systems of equations

We propose a new randomized method for solving systems of nonlinear equa...

A Direct Method to Learn States and Parameters of Ordinary Differential Equations

Though ordinary differential equations (ODE) are used extensively in sci...

A Data-Driven Line Search Rule for Support Recovery in High-dimensional Data Analysis

In this work, we consider the algorithm to the (nonlinear) regression pr...

Please sign up or login with your details

Forgot password? Click here to reset