Variational approaches provide a general way to derive many of the partial differential equations of physics [10]. They also are a starting point for devising efficient numerical solution techniques. We summarize the variational formulation of electrostatics next and show how this leads to flexible numerical techniques for solving the Poisson Eq. Similar ideas will be presented for the Poisson-Boltzmann (PB) Eq. in later sections.