MEMD with varying epsilon

This is the same pull request as before, but like Olaf suggested I added a unit test for setting the local permittivity in MEMD.

It now includes a (not perfect) testcase, documentation, and the code changes.