         vectorSubtract(&vx2, &vy2, vx1, vy1);
         // vectorMultiply(&vx2, &vy2, 0.5);
         vectorAdd(forceX, forceY, vx2, vy2);
+        if (gid == 0)
+          DEBUG_PRINT("a:%f\t%u\n", vx2, target);
       }
     }
   }

    float vlen = vectorGetLength(vx, vy);
    vectorInverse(&vx, &vy);
    vectorMultiply(&vx, &vy, K_G * (deg[gid] + 1) / vlen);
    vectorAdd(forceX, forceY, vx, vy);
  }
}


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

