Efficient Polarization Solvers For Classical Molecular Dynamics Simulations