def METHOD_NAME(self, s):
    corrector, corrector2, kernel, maxerror = s
    sim = rebound.Simulation()
    rebound.data.add_outer_solar_system(sim)
    for p in sim.particles:
        p.vx += 1.
    com = sim.calculate_com()
    sim.integrator = "whfast"
    sim.ri_whfast.corrector = corrector 
    sim.ri_whfast.corrector2 = corrector2
    sim.ri_whfast.kernel = kernel
    sim.dt = 0.0123235235*sim.particles[1].P  
    e0 = sim.energy()
    sim.integrate(1000.*2.*3.1415)
    e1 = sim.energy()
    self.assertLess(math.fabs((e0-e1)/e1),maxerror)
    com1 = sim.calculate_com()
    self.assertLess(math.fabs((com.x+com.vx*sim.t-com1.x)/(com1.x+com1.y)),1e-12)
    self.assertLess(math.fabs((com.y+com.vy*sim.t-com1.y)/(com1.x+com1.y)),1e-12)