def METHOD_NAME(self):
    self.system.part.add(pos=[0, 0, 0], q=1)
    self.system.part.add(pos=[0, 0, 1], q=-1)
    mmm1d = self.MMM1D(prefactor=1.0, maxPWerror=1e-20)
    self.system.electrostatics.solver = mmm1d
    self.assertTrue(mmm1d.is_tuned)
    self.system.integrator.run(steps=0, recalc_forces=True)
    self.check_with_analytical_result(prefactor=1.0, accuracy=0.0004)