Some of the early calculations of the polarization vector were based on 
direct inversion and the implementation of the conjugate gradient method by Petravic and Kuo-Petravic. Subsequently, many other conjugate gradient methods have been tested. Advances in the preconditioning of linear systems of equations arising in the DDA setup have also been reported.