Revolution EDA allows use of Harmonic balance simulation in Xyce to calculate all the relevant parameters such as IIP3, IIP2 and IP1dB even when the simulations do not converge when the circuit behaves strongly nonlinear at high RF power levels. Have a look at Https://
reveda.eu website for more information.