Add unit test for gravity calculation (rotation curve)
Unify gravity code between single-thread and OpenMP versions