Polarisation energy via iterative method is implemented but not tested. The forces are not implemented yet.