     float vlen = vectorGetLength(vx, vy);
     vectorInverse(&vx, &vy);
     vectorMultiply(&vx, &vy, K_G * (deg[gid] + 1) / vlen);
+    if (gid == 0)
+      DEBUG_PRINT("g:%f\n", vx);
     vectorAdd(forceX, forceY, vx, vy);
   }
 }

        vectorSubtract(&vx2, &vy2, vx1, vy1);
        // vectorMultiply(&vx2, &vy2, 0.5);
        vectorAdd(forceX, forceY, vx2, vy2);
      }
    }
  }


    fa2CleanEdgeMemory(&edgeData);
  }
  speedVectorFinish(*averageSpeedX, *averageSpeedY, numgraphs,
      vertices->numvertices);
  fa2CleanGeneralMemory(&data);
}

