The small shear rate response of electrorheological suspensions. I. Simulation in the point–dipole limit