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

[ViktorR]

Я так думаю, что быть экстрасенсом - неблагородное дело.
Почему вы не указываете что:

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

Произошло после того, как ...
Что вы сделали?

В предыдущем посте есть всё, что надо для понимания, а при необходимости, и изменения диапазона скоростей. Дерзайте.

[wqerty11]

ViktorR, Я сделал то, что Вы написали в прошлом посте.
А именно:
----------------------------------------
Пишем выражение:
Код:
x = -100 + rand()%(200 + 1)
и пробуем.
Т.е., выражение (целиком):
Код:
Vx = rand() % 100000 / 500.0 - 100;
заменяем на то, что написано выше.
------------------------------------------------

[ViktorR]

Ну и партизан.
А что с Vy? Тут тоже сделали замену?

x= -100 + rand()%(201)

Посмотрите как получилась эта формула и увеличьте диапазон скоростей.
Пробуйте.

[wqerty11]

ViktorR, Да, для Vy тоже сделал замену. Если пытаюсь увеличить значения скорости, опять начинаются проблемы с движением молекул.

[ViktorR]

А не сложно будет для вас показать, как вы увеличиваете скорости?

[wqerty11]

ViktorR, Я пытался увеличивать подчеркнутые значения.

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

[wqerty11]

ViktorR, Помогите, пожалуйста, решить проблему с увеличением скоростей молекул. Задание нужно сдать на проверку до конца недели. Спасибо!

[ViktorR]

Ещё раз.
У меня нет вашей среды для тестирования и нет времени всё это устанавливать (не говоря о языке ...).
Просто последовательно анализируйте данные, делайте выводы и корректируйте код.

Возможно, что кто-то ещё откликнется, хотя мне непонятно, почему другие эту тему игнорируют.
Может у них этот код просто работает без доп. усилий?

[wqerty11]

ViktorR, А что, если немного увеличить значение DT в этой строчке:

const float DT = 0.01;

Она используется в этой части программы:

void timer(int = 0)
{
 
   
    for (vector<P>::iterator i = p.begin(); i != p.end(); ++i)
        for (vector<P>::iterator j = p.begin(); j != p.end(); ++j)
            if (i != j)
            {
                float d = sqrt((i->x - j->x) * (i->x - j->x) + (i->y - j->y) * (i->y - j->y));
                if (d < i->m + j->m)
                {
                    float f = 50*(i->m + j->m - d);
                    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;
                }
            }
    for (vector<P>::iterator i = p.begin(); i != p.end(); ++i)
    {
        i->x += i->vx * DT;
        i->y += i->vy * DT;
        if (i->x < 0)
            i->x += 480;
        if (i->y < 0)
            i->y += 480;
        if (i->x > 480)
            i->x -= 480;
        if (i->y > 480)
            i->y -= 480;
    }
    display();
    glutTimerFunc(10, timer, 0);
}

Я установил значение DT = 0.03 вместо 0.01 и скорость молекул увеличилась, но они по прежнему двигаются нормально. Можно ли так поступить?

[ViktorR]

DT - это временной интервал, по истечение которого показывается следующая картинка (как я понял).
Дело в том, что в программе расчёт координат и скоростей ведётся с плавающей арифметикой, а рисовать надо с использованием целых чисел (координаты задаются в пикселах).
Увеличение DT - увеличение расстояний, на которые перемещаются молекулы.
Тут важным может оказаться тот момент, когда молекулы сталкиваются с большим телом (молекулой).
В изображении может проявиться ситуация, когда молекула влетает в большую.
Но этого не должно быть при соударении.

PS: Скорость молекул не увеличивается от DT.
От этого параметра увеличивается расстояние, на которое молекула улетает от текущей позиции.
Если текущие координаты X, Y, то следующие: X = X + Vx*DT и Y = Y + Vy*DT.

PSS: Действительно, если скорость составляет, например, 60 единиц, то за 0.01 молекула переместится на 0.6 единиц и, примерно, на один пиксел за 0.02 единицы времени.
Это значит, что формулу для получения скоростей можно поменять на

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

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

Если задавать рандомно, то будет много медленно движущихся молекул.

Если использовать нормальное распределение, то диапазон скоростей молекул будет пропорционален дисперсии, а средняя скорость составит, примерно, среднее значение такого распределения.
Алгоритм придётся поискать.
Например, есть трёхтомник Кнута, но для начинающего это будет сложно.
Мне неизвестны библиотеки в Си или С++, которые позволяют генерировать такие наборы быстро.