This procedure was repeated 200 times. Factors were initialized by the final solution values, and all auxiliary equations were maintained active. The re-sampling was performed first by generating a matrix of weights with Poissonian distribution ri. Only the rows corresponding to a positive ri are sampled.