Numerical Algorithms as Dynamcal Systems - Structure Preservation, Convergence Theory, and Rediscretization

该项目的重点是通过将数值分析中的迭代算法与特定的微分方程系统联系起来来改进迭代算法。三十年前，人们意识到计算特征值的QR算法可以看作是Toda格的时间- t映射。这开辟了将动力系统的定性方法引入算法的可能性，并有可能通过更有效的离散化来加速算法。Toda格的关键方面是它是一个连续共轭，它保留了上海森伯格形式，因此保留了特征值。在最近的研究中，PI通过开发保留双对角结构的Lotka-Volterra系统，在奇异值分解的计算中取得了类似的结果，这导致了SVD计算的进步。本项目研究了其他保持辛和哈密顿结构的动力系统，这在许多应用领域都是关键的。最终目标是研究其几何结构与现有数值算法之间的联系，建立严格的动力学行为数学理论，并开发可能的结构保持再离散化，以提高数值分析迭代的鲁棒性、速度和准确性。保存结构的动力系统是自然而普遍存在的。物理世界中的守恒定律和工程系统中的约束力学只是两个例子。结构保存在计算中也是必要的，因为它使更有效的算法成为可能，提高物理可行性和可解释性，并且在长期行为中更健壮。提出的研究将数值算法重塑为模拟相应迭代方案的结构保持特性的微分系统。理解连续统系统的整体动力学可以揭示相关离散系统的收敛性质，也有助于将连续统系统重新离散化为具有更好数值性质的新算法。广泛的应用将受益于这些系统的性质和几何结构的研究，这些研究基本上包括所有需要结构保存的学科，包括经典力学和量子力学、模型简化、可逆系统和分子动力学。

The focus of the project is to improve iterative algorithms in numerical analysis by linking them to particular systems of differential equations. Thirty years ago, it was realized that the QR algorithm for calculating eigenvalues can be considered as a time-T map of the Toda lattice. This opens up the possibility of bringing the qualitative methods of dynamical systems to bear on the algorithm, and to potentially speed up the algorithm by more efficient discretizations. The crucial aspect of the Toda lattice is that it is a continuous conjugation that preserves upper Hessenberg form and therefore the eigenvalues. In recent research, the PI has achieved similar results with the calculation of the singular value decomposition by developing a Lotka-Volterra system that preserves bidiagonal structures, which has led to advances in computing the SVD. This project investigates other dynamical systems preserving symplectic and Hamiltonian structure, which are pivotal in many areas of applications. The ultimate goal is to investigate the connection between their geometric structures and existing numerical algorithms, to establish a rigorous mathematical theory on dynamical behaviors, and to develop possible structure-preserving re-discretizations to improve robustness, speed, and accuracy of iterations in numerical analysis. Structure-preserving dynamical systems are natural and ubiquitous. Conservation laws in the physical world and constrained mechanics in engineered systems are just two examples. Structure preservation is also imperative in computation because it makes possible more efficient algorithms, improves physical feasibility and interpretability, and is more robust in long-term behavior. The proposed research recasts numerical algorithms as differential systems that mimic the structure-preserving properties of the corresponding iterative schemes. Understanding the overall dynamics of the continuum system can shed light on the convergence properties of the related discrete counterparts, and can also contribute to the re-discretization of the continuum system into a new algorithm with better numerical properties. A wide range of applications stands to benefit from the study of properties and geometric structure of these systems, which essentially include all disciplines that entail structure preservation, including classical and quantum mechanics, model reduction, reversible systems, and molecular dynamics.