def METHOD_NAME(self):
    import espressomd.electrostatics
    self.system.time_step = 0.01
    self.add_charged_particles()
    for node_grid, p3m_params in FFT_PLANS[self.n_nodes]:
        self.system.cell_system.node_grid = node_grid
        solver = espressomd.electrostatics.P3M(
            prefactor=2, accuracy=1e-6, tune=False, **p3m_params)
        self.system.electrostatics.solver = solver
        ref_energy = -75.871906
        p3m_energy = self.system.analysis.energy()['coulomb']
        self.system.electrostatics.clear()
        np.testing.assert_allclose(p3m_energy, ref_energy, rtol=1e-4)