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

[wqerty11]

ViktorR, Задача нужна для университета. И я знаю как выполняется вывод на экран, просто не понятен этот момент в тех местах, где они меняются

[ViktorR]

void timer(int = 0)
{
 ...
                    i->vx += f * (i->x - j->x) / d / i->m * DT;
                    i->vy += f * (i->y - j->y) / d / i->m * DT;
                    j->vx -= f * (i->x - j->x) / d / j->m * DT;
                    j->vy -= f * (i->y - j->y) / d / j->m * DT;
...
}
 
int main(int argc, char** argv)
{
    P p0 = { 480 / 2, 480 / 2, 0, 0, 40 };
    p.push_back(p0);
    for (int i = 0; i < 100; ++i)
    {
        P p0 = { rand() % 480, rand() % 480, rand() % 100000 / 500.0 - 100, rand() % 100000 / 500.0 - 100, 1 };
        p.push_back(p0);
    }

Вот что делает код в выделенных строках?
В одном случае параметры вектора инициализируются, а потом изменяются в void timer(int = 0).
Вот надо и посмотреть, для небольшого числа шариков, как выполнилась инициализация (какие числа получились) и как они были изменены.
Если все начальные параметры совпадают при инициализации, то как работает rand()?
Если тут всё нормально, то как работает процедура изменения параметров?
Почему ваши частицы, как вы пишете:

Молекулы постоянно летят из правого нижнего угла окна в левый верхний.

Смотрите код пошагово. Разберитесь с каждой операцией.

[wqerty11]

ViktorR, Здравствуйте. Попытался разобраться в коде, и выяснил, что проблема в выделенных значениях. Чем они больше, тем хуже летят молекулы. Например, если значение 10, то всё вроде в порядке. В чем может быть проблема?

P p0 = { rand() % 480, rand() % 480, rand() % 100000 / 500.0 - 100, rand() % 100000 / 500.0 - 100, 1 };

[wqerty11]

wqerty11, Короче, проблема 100% в этой строке, осталось понять как её исправить.

[ViktorR]

Вы написали, что у вас

Значение RAND_MAX=32767

rand() % 100000 - это, как я понял из описания, генерация с масштабирующим множителем.
В этом случае формируется случайное число в интервале от 0 до 100000.
У меня нет возможности это проверять какие числа генерируются, а вы не желаете сделать вывод характеристик нескольких векторов.
В чём затруднение? Почему нельзя сделать вывод характеристик?
Ведь код можно поправить и написать так:

for (int i = 0; i < 5; ++i)   // только для пяти точек
    {
      Xc = rand() % 480;
      Yc = rand() % 480;
      Vx =  rand() % 100000 / 500.0 - 100;
      Vy =  rand() % 100000 / 500.0 - 100;
      <тут сделать вывод на экран этих параметров>

        P p0 = { Xc, Yc, Vx, Vy, 1 };
        p.push_back(p0);
    }

Сделайте вывод, посмотрите на него, и поймите, почему так работает.

[wqerty11]

ViktorR, Вот что получилось:

[ViktorR]

Т.е., все проекции скоростей отрицательные. Это значит, что частицы двигаются влево и вверх.
rand() % 100000 / 500.0 - всегда меньше 100, которую отнимаем.
Разберитесь с генератором случайных чисел.
rand() % 100000 - это какой диапазон генерируемых чисел?

Какой диапазон скоростей мог бы быть в вашей модели?
Возможно, что можно делить не на 500.0, а 50.0.
Вот нашёл, например:
генерация в диапазоне [a, b):

x = a + rand()%(b-a)  // Правая граница открытая: меньше b

генерация в диапазоне [a, b]:

x = a + rand()%(b-a+1)  // Правая граница закрытая: меньше или равно (не больше) b

Прими для себя некоторые значения для a и b (например, -20, 20) и проверь.

PS: Не стесняйся, заменяй предложенные в видео варианты на свои.
Экспериментируй. Не факт, что автор видео крутой и кудрявый.
Возможно, что он что-то своё двигал (ну так на тот момент сложилось).

[wqerty11]

ViktorR, У автора видео в компиляторе RAND_MAX = 2147483647, а у меня 32767. Думаю проблема в этом.

[ViktorR]

Не заморачивайся этим видео. Пробуй свои варианты. Осваивайся.
Там не всё так тривиально, поскольку в выражении rand() % 100000
100000 - это масштабирующий коэффициент.
Т.е., если я правильно понял, то генерируются числа в диапазоне от 0 до 100000 в независимости от
значения RAND_MAX.
После деления на 500.0 получаем диапазон от 0 до 200.
Вычитая 100 должны получить распределение от -100 до 100.
Пишем выражение:

x = -100 + rand()%(200 + 1)

и пробуем.
Т.е., выражение (целиком):

Vx =  rand() % 100000 / 500.0 - 100;

заменяем на то, что написано выше.

[wqerty11]

ViktorR, Движение частичек вроде наладилось, но их скорость теперь достаточно низкая.