A-stable and L-stable high-order integration methods for solving stiff differential equations

  • Authors:
  • Emad Gad;Michel Nakhla;Ramachandra Achar;Yinghong Zhou

  • Affiliations:
  • School of Information Technology and Engineering, University of Ottawa, Ottawa, ON, Canada;Department of Electronics, Carleton University, Ottawa, ON, Canada;Department of Electronics, Carleton University, Ottawa, ON, Canada;Department of Electronics, Carleton University, Ottawa, ON, Canada

  • Venue:
  • IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems
  • Year:
  • 2009

Quantified Score

Hi-index 0.03

Visualization

Abstract

This paper describes a new A- and L-stable integration method for simulating the time-domain transient response of nonlinear circuits. The proposed method, which is based on the Obreshkov formula, can be made of arbitrary high order while maintaining the A-stability property. The new method allows for the adoption of higher order integration methods for the transient analysis of electronic circuits while enabling them to take larger step sizes without violating the stability, leading to faster simulations. The method can be run in an L-stable mode to handle circuits with extremely stiff equations. Necessary theoretical foundations, implementation details, error-control mechanisms, and computational results are presented.