Deflated Restarting of Exponential Integrator Method with an Implicit Regularization for Efficient Transient Circuit Simulation

Exponential integrator (EI) method based on Krylov subspace approximation is a promising method for large-scale transient circuit simulation. However, it suffers from the singularity problem and consumes large subspace dimensions for stiff circuits when using the ordinary Krylov subspace. Restarting schemes are commonly applied to reduce the subspace dimension, but they also slow down the convergence and degrade the overall computational efficiency. In this paper, we first devise an implicit and sparsity-preserving regularization technique to tackle the singularity problem facing EI in the ordinary Krylov subspace framework. Next, we analyze the root cause of the slow convergence of the ordinary Krylov subspace methods when applied to stiff circuits. Based on the analysis, we propose a deflated restarting scheme, compatible with the above regularization technique, to accelerate the convergence of restarted Krylov subspace approximation for EI methods. Numerical experiments demonstrate the effectiveness of the proposed regularization technique, and up to 50% convergence improvements for Krylov subspace approximation compared to the non-deflated version.

[1]  Quan Chen,et al.  A new tightly-coupled transient electro-thermal simulation method for power electronics , 2016, 2016 IEEE/ACM International Conference on Computer-Aided Design (ICCAD).

[2]  C. Tischendorf,et al.  Structural analysis of electric circuits and consequences for MNA , 2000 .

[3]  L. Bergamaschi,et al.  Interpolating discrete advection-diffusion propagators at Leja sequences , 2004 .

[4]  Xiang Zhang,et al.  Simulation Algorithms With Exponential Integration for Time-Domain Analysis of Large-Scale Power Delivery Networks , 2015, IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems.

[5]  Quan Chen,et al.  Circuit simulation via matrix exponential method for stiffness handling and parallel processing , 2012, 2012 IEEE/ACM International Conference on Computer-Aided Design (ICCAD).

[6]  Quan Chen,et al.  A fast time-domain EM-TCAD coupled simulation framework via matrix exponential , 2012, 2012 IEEE/ACM International Conference on Computer-Aided Design (ICCAD).

[7]  Awad H. Al-Mohy,et al.  Computing the Action of the Matrix Exponential, with an Application to Exponential Integrators , 2011, SIAM J. Sci. Comput..

[8]  V. Il'in Projection Methods in Krylov Subspaces , 2019, Journal of Mathematical Sciences.

[9]  André Gaul,et al.  Recycling Krylov subspace methods for sequences of linear systems , 2014 .

[10]  Y. Saad Analysis of some Krylov subspace approximations to the matrix exponential operator , 1992 .

[11]  Timo Reis,et al.  Surveys in Differential-Algebraic Equations III , 2015 .

[12]  Scott H. Cameron PIECE-WISE LINEAR APPROXIMATIONS , 1966 .

[13]  Oliver G. Ernst,et al.  A Restarted Krylov Subspace Method for the Evaluation of Matrix Functions , 2006, SIAM J. Numer. Anal..

[14]  Quan Chen,et al.  Time-Domain Analysis of Large-Scale Circuits by Matrix Exponential Method With Adaptive Control , 2012, IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems.

[15]  Leon O. Chua,et al.  Computer-Aided Analysis Of Electronic Circuits , 1975 .

[16]  Chung-Kuan Cheng,et al.  MATEX: A distributed framework for transient simulation of power distribution networks , 2014, 2014 51st ACM/EDAC/IEEE Design Automation Conference (DAC).