Preconditioned Iterative Solves in Model Reduction of Second Order Linear Dynamical Systems
From MaRDI portal
Publication:6274271
arXiv1606.01216MaRDI QIDQ6274271FDOQ6274271
Authors: Navneet Pratap Singh, Kapil Ahuja, Heike Faßbender
Publication date: 3 June 2016
Abstract: Recently a new algorithm for model reduction of second order linear dynamical systems with proportional damping, the Adaptive Iterative Rational Global Arnoldi (AIRGA) algorithm, has been proposed. The main computational cost of the AIRGA algorithm is in solving a sequence of linear systems. These linear systems do change only slightly from one iteration step to the next. Here we focus on efficiently solving these systems by iterative methods and the choice of an appropriate preconditioner. We propose the use of relevant iterative algorithm and the Sparse Approximate Inverse (SPAI) preconditioner. A technique to cheaply update the SPAI preconditioner in each iteration step of the model order reduction process is given. Moreover, it is shown that under certain conditions the AIRGA algorithm is stable with respect to the error introduced by iterative methods. Our theory is illustrated by experiments. It is demonstrated that SPAI preconditioned Conjugate Gradient (CG) works well for model reduction of a one dimensional beam model with AIRGA algorithm. Moreover, the computation time of preconditioner with update is on an average 2/3 rd of the computation time of preconditioner without update. With average timings running into hours for very large systems, such savings are substantial.
Iterative numerical methods for linear systems (65F10) Transformation and reduction of ordinary differential equations and systems, normal forms (34C20) Stability and convergence of numerical methods for ordinary differential equations (65L20)
This page was built for publication: Preconditioned Iterative Solves in Model Reduction of Second Order Linear Dynamical Systems
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6274271)