The Jordan Structure of Residual Dynamics Used to Solve Linear Inverse Problems
Chein-Shan Liu1, Su-Ying Zhang2, Satya N. Atluri3
Department of Civil Engineering, National Taiwan University, Taipei, Taiwan. E-mail: liucs@ntu.edu.tw
College of Physics and Electronic Engineering, Shanxi University, Taiyuan, China.
Center for Aerospace Research & Education, University of California, Irvine
With a detailed investigation of n linear algebraic equations Bx=b, we find that the scaled residual dynamics for y∈Sn−1 is equipped with four structures: the Jordan dynamics, the rotation group SO(n), a generalized Hamiltonian formulation, as well as a metric bracket system. Therefore, it is the first time that we can compute the steplength used in the iterative method by a novel algorithm based on the Jordan structure. The algorithms preserving the length of y are developed as the structure preserving algorithms (SPAs), which can significantly accelerate the convergence speed and are robust enough against the noise in the numerical solutions of ill-posed linear inverse problems.
Liu, C., Zhang, S., Atluri, S. N. (2012). The Jordan Structure of Residual Dynamics Used to Solve Linear Inverse Problems. CMES-Computer Modeling in Engineering & Sciences, 88(1), 29–48.
This work is licensed under a Creative Commons Attribution 4.0 International License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.