Abstract In the evaluation of the correlation by means of second-order perturbation theory, the perturbative evaluation of the coefficient of each doubly excited configuration is replaced by the diagonalization of a minimally dresed 2×2 matrix. The dressing incorporates the effect of high-order exclusion principle violating diagrams in a self-consistent way. An efficient computation of those coefficients goes through the storage of intermediate summations in three-dimensional arrays. This method, the cost of which is only a few MP2 calculations, does not diverge when single bonds are broken, as shown on a few examples (Li 2, F 2). Its generalization to complete active space CI for multiple bond breakings (and even to selected CIs) is straightforward.