We propose a generalization of nonlinear stability of numerical one-step integrators to Riemannian manifolds in the spirit of Butcher's notion of B-stability. Taking inspiration from Simpson-Porco and Bullo, we introduce non-expansive systems on such manifolds and define B-stability of integrators. In this first exposition, we provide concrete results for a geodesic version of the Implicit Euler (GIE) scheme. We prove that the GIE method is B-stable on Riemannian manifolds with non-positive sectional curvature. We show through numerical examples that the GIE method is expansive when applied to a certain non-expansive vector field on the 2-sphere, and that the GIE method does not necessarily possess a unique solution for large enough step sizes. Finally, we derive a new improved global error estimate for general Lie group integrators.
翻译:我们提出将数值单步积分器的非线性稳定性推广到黎曼流形上,沿袭Butcher关于B稳定性的概念。受Simpson-Porco和Bullo的启发,我们在这类流形上引入非扩张系统,并定义了积分器的B稳定性。在首次阐述中,我们针对隐式欧拉(GIE)方案的测地线版本给出了具体结果。我们证明了在具有非正截面曲率的黎曼流形上,GIE方法是B稳定的。通过数值算例表明,当应用于2-球面上的特定非扩张向量场时,GIE方法具有扩张性;且对于足够大的步长,GIE方法不一定存在唯一解。最后,我们推导出一般李群积分器的一个新的改进全局误差估计。