Быстрые nbody алгоритмы / решения (с opengl/c ++/??)

Я продолжаю работать (C++, opengl) проект, где у меня должно быть много частиц, которые влияют друг на друга, если я корректен, это называют nbody проблемой. Делает кто-то знает то, что решения там для алгоритмов как это.

Я знаю barnes алгоритм хижины, и возможно я могу посмотреть вокруг openCL, хотя я только задаюсь вопросом, использовали ли Вы, возможно, другие решения.

Код, который я создам, будет иметь много из:

for(int i = 0; i < num_particles; ++i) {
  for(int j = i+1, j < num_particles; ++j)
     dist = distance(particles[i],particles[j]);
     if(dist > limit) {....}
  }
}

С наилучшими пожеланиями, Поллукс