In this paper we consider the numerical solution of the large system of nonlinear differential equations. We assume that the system simulates the semiconductor circuit. We apply the well known event driven techniques to get some approximation of solution fast. We extend those techniques by considering pairs of nodes, instead of single nodes, as usual. The “pairwise” solution is more efficient in tightly coupled circuits. The improvement of efficiency of solution is important optimizing the parameters of circuits. In many cases we need the global optimization of equation parameters. Then, the modeling speed could be the main factor of successful application.