Reliable Computer Simulation Methods for Electrostatic Biomolecular Models Based on the Poisson–Boltzmann Equation