ПривеТ! Решил поделиться с читателями своими небольшими экспериментами с системами частиц в трехмерном пространстве. За основу взял публикацию на Хабре об экспериментах с частицами в 2D пространстве.
Начнем со ссылки на статью которая меня и подтолкнула к действию. Но была еще одна причина, с недавнего времени все больше перехожу в своих экспериментах на Ubuntu, многое радует начиная со свободной установки ОС и далее по списку плюсов. Минусы есть, у меня это борьба с драйверами при нестандартной установке ОС вроде двух видеокарт и нескольких мониторов.
За основу взял с++, добавил поддержку CUDA как вычислительной платформы, частиц много и центральный процессор явно не справится и графический движок Ogre3D к ним в компанию. Кашу заварили, буду приправлять повествование Gif анимацией и цитатами из статьи донора о 2D варианте.
"Сначала я пошёл по стопам игры «жизнь»: у каждой частицы есть счётчик “перенаселения”, который равен сумме обратных квадратов расстояний до других частиц. Если этот счётчик меньше определённого предела, то есть соседей мало, то частица притягивается к другим частицам, а если соседей много — отталкивается. Если частицы пересекаются, то они отталкиваются в любом случае, чтобы не проходить сквозь друг друга.
Случайно раскидываем частицы по полю и смотрим, что выйдет."
Расширенный вариант видео
Теперь немного о логике происходящего в программе. Создается массив частиц с определёнными параметрами. Часть параметров отвечают за физические свойства: радиус, масса. Часть за создание связей между частицами, такими как тип частиц, кол-во связей частицы с другими частицами и тд. В среднем использовал в симуляции от 1000 до 3000 частиц, большее значение приводило к торможению картинки, так как хотелось считать в реальном времени. Поэтому все это добро передаётся в память видеокарты, где уже GPU в режиме сильного распараллеливания обрабатывает три основных подпрограммы: симуляция движения частиц, обработка столкновений частиц(соударений) и образование-разрушение связей между частицами.
"Меняем правила игры. Больше не будем считать соседей. Пусть частицы будут просто притягиваться или отталкиваться в зависимости от их типов. Если все частицы одного типа, то тут всего 2 варианта: они либо все отталкиваются, либо все притягиваются."
"Не будем сильно менять правила. Вместо этого добавим новую фичу. Теперь частицы на небольшом расстоянии будут образовывать связи. Если частицы связаны, то они постоянно притягиваются друг к другу. Это притяжение не ослабевает с расстоянием. Но, если расстояние выше определённого порога, то связь рвётся."
Далее начался процесс проб и ошибок. Во первых программирование на GPU требует дополнительной внимательности в части так называемого "races on read-modify-write of shared data", что означает что могут быть проблемы когда несколько потоков одновременно пытаются изменить переменную. Возникают нестабильности вроде таких:
Затем надо было ограничить область пространства в котором происходит эксперимент, первое что пришло на ум это куб. Но спасибо багам первых версий программы, частицы смело расползались из него, образуя что-то подобное космическим станциям из фантастики 70-х.
Как говорится если не куб, то пусть будет шар:
Здесь логика в том, что после отлёта от центра на частицу начинает действовать сила обратная ее движению, что даже похоже на гравитацию.
__device__ static int Links[3][3] = {{1,0,1},{1,0,0},{1,1,0}};
__device__ static int LinksField[3][3] = {{0,0,0},{0,0,1},{0,0,0}};
__device__ static int LinkTypeSize[3] = {3,8,2};
__device__ static int LinkTypePP[3][3] = {{0,1,1},{8,1,6},{0,1,1}};Привел кусочек кода как есть, это матрицы взаимодействия трех типов частиц друг с другом.
-первая строчка это принцип образования связей: 1 есть возможность установить связь с другой частицей, 0 соответственно нет.
-вторая строчка- это закон притяжения отталкивания частиц. По сути в этом варианте все отталкиваются друг от друга, кроме притяжения второго типа частиц к третьему. Можно конечно зеркалить матрицу, но в данном случае так.
-третья строчка количество общих связей частиц по типам.
-четвёртая это количество связей частицы с каждым типом частиц.Первый тип например сам с собой не может образовать связь.
Получает что-то такое в сильном движении:
расширенная версия видео
Пришлось вводить потерю энергии при столкновениях.
Стало слишком все статично, уменьшаем силы взаимодействия частиц.
Если посмотреть видео видно как происходит постройка каркаса. Потом он остаётся все равно статичным.
Меняем правила в части коэфф матриц.
расширенная версия видео
Видим образование подобия органических молекул. Даже есть бензольные кольца.
Но надо попробовать повторить вариант 2D автора, чьи цитаты в теле статье. Есть вариант такой "жизни"
или такой
Что можно добавить. Исходники.
И еще ролик, на самом деле завораживающая картина, то кажутся молекулы, то нейросети, но жизни пока там нет. Как вариант развития генетику нужно добавлять, но непонятно, что может являться фитнес функцией для этих "существ".
Надеюсь кому-то понравится и он шагнёт из 3D дальше, как автор шагнул из 2D. Главное для себя повторил уроки по CUDA и Ogre3D под Ubuntu.
И еще будут идеи пишите, может что интересное придумаем )
Автор: Geo Evclid
