Harmonic averaging of smooth permittivity functions in finite-difference Poisson–Boltzmann Electrostatics