Iterative solution of linear systems in harmonic balance analysis

Harmonic balance (HB) is a steady-state simulation technique of primary interest for RF and microwave circuits. Krylov subspace methods promise efficient solution of the large linear systems that arise in HB simulators. This paper deals with an experimental investigation of GMRES and QMR, two leading Krylov subspace methods as applied to the HB problem. The problem of coordinating the linear solver´s accuracy with the error at the nonlinear level is also discussed.