Проблема в броуновском движении молекул C++

[digitalis]

Нормальное распределение моделируется легко:
http://stratum.ac.ru/education/textb...lection25.html

[wqerty11]

ViktorR, Если я увеличиваю DT, то маленькие молекулы не пролетают через большую, а если увеличиваю скорость вот так:

x = -500 + rand()%(1001)

то маленькие молекулы начинают пролетать через большую, как Вы и скзаали.

[BDA]

У меня тоже RAND_MAX=32767. Изменил в коде из первого сообщения строку инициализации молекул на такую:

P p0 = { rand() % 480, rand() % 480, rand() % 10000 / 50.0 - 100, rand() % 10000 / 50.0 - 100, 2 };

Стали лететь в разные стороны. Но наблюдается проблема залета мелких молекул в крупную.

[ViktorR]

Алгоритмизация модели столкновения может показаться простой, но это далеко не так.
Почитал несколько статей и понял, что для небольшого числа тел её, от части, можно решить путём перебора.
Для большого числа тел используют ухищрения.
Например, как понял я, разбивая область на параллелограммы (или прямоугольники, или параллелепипеды - 3D: определяется используемой геометрией пространства).
Это позволяет составить 2-D или 3-D массив, в котором элементы массива - ID тел (или ссылка на список ID, если таких тел собралось в одной области несколько).
Тогда, для расчёта столкновений, достаточно просмотреть только соседние элементы массива.
Сам массив генерируется сортировкой перед просмотром столкновений.

Ещё сложнее задача будет решаться, если надо учитывать не только первый уровень (тело сталкивается с телом), но и второй - в столкновении участвует больше 2-х тел.

Так подумалось, что, для упрощения отображения, можно пересчитывать координаты тел и не отображать сами моменты столкновения, а показывать состояние системы через заданные интервалы.
Т.е., если обнаружено столкновение (тело влетело в другое), то пересчитываем их положение и затем показываем.
Правда как это будет видно визуально? Тело летело влево, а теперь летит вправо вверх, например, (излом какой-то).

PS: Ссылок достаточно много и нет смысла их тут приводить.
Для ТС: Вам надо самому принять, что должна показывать ваша модель.
Тут важнее то, что при столкновении импульс передаётся большой молекуле и эта молекула совершает случайные движения.