         ux = momentum * ux - eta * gx * dx;
         uy = momentum * uy - eta * gy * dy;
 
-        pts[i] += ux;
-        pts[i + nnodes + 1] += uy;
+        pts[i] += -eta * dx;
+        pts[i + nnodes + 1] += -eta * dy;
 
         old_forces[i] = ux;
         old_forces[N + i] = uy;

                        volatile float * __restrict old_forces) // (N x 2)
{
  register int i, inc;
  register float dx, dy, ux, uy, gx, xy;

  // iterate over all bodies assigned to thread
  // TODO: fix momentum at step 0

}


__global__
__launch_bounds__(1024, 1)
void ClearKernel2(int nnodesd, volatile int * __restrict startd, volatile float * __restrict massd)
{
