“…The calculations are fast and accurate and in obvious contrast to the Ewald summation methods where one needs to evaluate at least 100 -200 terms, involving error func- tions etc., the calculations are also simpler. In particular, if one is doing Monte Carlo simulations for a magnetic array the positions of the particles remain fixed for the entire simulation and only the directions of the magnetic moment vector changes in each step of the simulation [21]. In that case the space summations A, B, C, D, E, F and G involved in the dipole interaction energy between each pair of particles i and j in the array can be stored and recalled in each MC step.…”