  A modified 3D-Ewald summation is presented for accurately simulating the
ion-dipole mixture under dielectric confinement. The method is based on the
combination of image charges and image dipoles with the conventional Ewald
summation and has a scaling O(^3/2). The accuracy and efficiency of our
algorithm are examined through numerical examples.
