А если тупо ребра лить в файл? Для передачи энергии от одной точки к другой - вроде устраивает, просто проход по файлу. Правда извлечь/получить ребра "по запросу" будет большая забота, ну пока неясно нужно ли это.

Ну и вроде "фокусируется" на том как грамотно создавать/объединять точки. Здесь наверно лучше сразу планировать "как задействовать все ядра"

Ну хорошо, вот маленькая подзадачка, чисто математика. Граф задуман для "перекачки энергии" между точками, псевдокод

C++ (Qt)
void Transfer( Point * p0, Point * p1 )
{
  auto e01 = p0->Energy() * p0->GetIn() * p1->GetOut();
  auto e10 = p1->Energy() * p1->GetIn() * p0->GetOut();
  p0->AddEnergy(e10);
  p1->AddEnergy(e01);
}
 

Для простоты геттеры просто возвращают константы. Да, но это одна итерация. Вопрос: могу ли я учесть все итерации (напр 5) за одно вычисление/проход ? Или нужно 5 проходов (что не так уж дешево)

Спасибо

Я тут на один канал подписан, там, в частности, есть проблема графа.. Чувак просто чудеса творит https://www.youtube.com/user/foo52ru
Может кому интересно будет

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

Ну вот опять пример непоставленной проблемы.. Что Вы ожидаете, от тех, кто не в этой свферы находится, какое решение?

Ну вот, мучить детей "тензорным анализом" - то да, а как нужна простая но необходимая вещь - сразу смылся   Ось малюнок (аттач). Лучи показаны точками, типа "каждый сидит в своей ячейке", размеры ячеек неодинаковы, на краях они заметно больше, это и соответствует тому что по нормали точка получает больше света. Вес всех лучей одинаков, рез-т получается простым осреднением.

Да, но как посчитать размеры ячеек чтобы перейти к конкретному коду? На картинке есть подсказка

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

1) Из точки выбросить заданное кол-во случайных лучей M так чтобы они группировались вокруг заданной нормали N как показано в моем предыдущем посте - по нормали больше лучей, по касательной меньше. Т.е. требуется заменить "больше/меньше" на конкретные формулы, а еще лучше на алгоритм, код напишем - не проблема. Вот тогда все увидят какая это красивая и мощная дисциплина - тензорный анализ 

2) Как всегда, направления лучей нужны случайные, но .. не совсем. Нужно уметь выбрасывать в центр каждой ячейки, а потом случайно изменить напр-е но так чтобы луч остался в той же ячейке. Так мы сможем парить адаптивно, а главное - "квалифицировать" луч выброшенный из др точки (в какой он ячейке)

Точка на поверхности единичной сферы с центром (0, 0, 0) и есть компоненты луча, т.е. это не решение, а повтор задачи др словами (плюс понты "в чем проблема")

Ну вот опять - запутать, запугать, задрочить простейшую вещь так чтобы простой человек почесал репу и сказал "не, ну нафиг, я программист, а не математик" Народная примета: если автор статьи резво сыпет интегралами - скорее всего пустышка.
Еще мерзкая черта - не давать определений. Кто такие тета и фи? Откуда человек "от сохи" (читай "сборки") может их знать?

В общем, вежливо (здесь научным работникам надо отдать должное), но настойчиво внушается мысль что все это, мол, "не программистского ума дело". Увы, это стандартный подход, Вы не хуже остальных.

Однако это совсем не так. Все не то что "просто", а "очень просто" если уловить идею решения. Подсказка: внизу нарисована плоская сетка, проекция на касательную плоскость. Площади проекций всех ячеек одинаковы (увы, никто не допер ). Оказывается проекции ячеек/лучей распределены равномерно внутри круга (сами ячейки конечно нет). Дальше "элементарщина", для нормали по Y

C++ (Qt)
QVector3D RandomDiffuseRay( void )
{
// случайная длина в плоскости XZ
  float R = RandomFloat();
 
// случайное напр-е - угол в плоскости XZ
  float angle = RandomFloat() * M_PI * 2;
 
  return QVector3D(cos(angle) * R, sqrt(1 - R * R), sin(angle) * R);
}

Найти в какую ячейку попал луч - тоже просто, надо только запастись табличкой "широт" для заданного разрешения, и потом по ней lower_bound(y). А "долгота" = atan2(z, x). Правда это для нормали (0, 1, 0), но до заданной уж как-то довернем, там "достаточно пройти дальше по ссылке" (як каже молодь  )

Как видим, никакие интегралы не нужны. Зачем же их упорно рисуют?