def METHOD_NAME(self):
    N = 200
    system = self.system
    system.time_step = 0.01
    partcls = system.part.add(pos=system.box_l * np.random.random((N, 3)))
    kT = 2.3
    gamma = 1.5
    system.thermostat.set_dpd(kT=kT, seed=42)
    system.non_bonded_inter[0, 0].dpd.set_params(
        weight_function=0, gamma=gamma, r_cut=1.5,
        trans_weight_function=0, trans_gamma=gamma, trans_r_cut=1.5)
    system.integrator.run(10)
    system.thermostat.turn_off()
    # Reset velocities
    partcls.v = [1., 2., 3.]
    system.integrator.run(10)
    # Check that there was neither noise nor friction
    for v in partcls.v:
        for i in range(3):
            self.assertEqual(v[i], float(i + 1))
    # Turn back on
    system.thermostat.set_dpd(kT=kT, seed=42)
    # Reset velocities for faster convergence
    partcls.v = [0., 0., 0.]
    # Equilibrate
    system.integrator.run(250)
    loops = 250
    v_stored = np.zeros((loops, N, 3))
    for i in range(loops):
        system.integrator.run(10)
        v_stored[i] = partcls.v
    v_minmax = 5
    bins = 5
    error_tol = 0.012
    self.check_velocity_distribution(
        v_stored.reshape((-1, 3)), v_minmax, bins, error_tol, kT)