The implementation of an improved small-signal AC simulation capability in the general-purpose device simulator PISCES-IIB is described. The preconditioned generalized conjugate residual (GCR) algorithm has been implemented, which allows AC simulations to be performed up to any frequency without convergence problems, although at great computational expense. The current implementation of PISCES-IIB uses the block successive overrelaxation algorithm for AC simulations. This algorithm fails at high frequencies, thus making it impossible to determine accurately the cutoff frequency and switching speed of high-frequency devices. The preconditioned GCR algorithm has been implemented using two different preconditioners. A comparison of the three methods shows that they are most efficient at different frequency ranges, allowing programming of an automatic switching algorithm that chooses the most efficient simulation method.>