Simulating coarse-grained models of charged soft-condensed matter systems in presence of dielectric discontinuities between different media requires an efficient calculation of polarization effects. This is almost always ...