A Multi-Stage Adaptive Sampling Scheme for Passivity Characterization of Large-Scale Macromodels