КОНВЕЙР РАЗУМА

Пусть дан взвешенный неориентированный граф G=(V,E), |V|=n. Весовая матрица:

W=(w_{ij})_{i,j=1}^n,\quad w_{ij}=w_{ji}\ge 0,\quad w_{ii}=0.

Степени (взвешенные):

d_i=\sum_{j=1}^n w_{ij},\qquad D=\mathrm{diag}(d_1,\dots,d_n).

1.1. Ненормированный лапласиан

L = D - W.

1.2. Нормированный лапласиан (симметричный)

L_{\mathrm{sym}} = I - D^{-1/2} W D^{-1/2}.

1.3. Нормированный лапласиан (случайный/марковский)

Определим марковскую матрицу переходов:

P = D^{-1} W, тогда L_{\mathrm{rw}} = I - P.

Это генератор дискретного “теплопереноса” по графу.

2) Графовая диффузия

2.1. Дискретное время (случайное блуждание)

Пусть ut\in\mathbb{R}^n — распределение/поле на вершинах в момент t.

u_{t+1} = P\,u_t. 

Эквивалентно:

u_{t+1}-u_t = -(I-P)u_t = -L_{\mathrm{rw}}u_t.

2.2. Непрерывное время (графовое “тепловое уравнение”)

\frac{d u(t)}{dt} = -L\,u(t)\quad \text{или}\quad \frac{d u(t)}{dt} = -L_{\mathrm{rw}}u(t).

Решение:

u(t)=e^{-tL}u(0)\quad \text{или}\quad u(t)=e^{-tL_{\mathrm{rw}}}u(0).

Интерпретация: L (или L_{\mathrm{rw}}) — дискретный аналог оператора Лапласа: он измеряет “насколько значение в узле отличается от среднего по соседям”.

3) Настройка для предельной теоремы (точки на многообразии)

Пусть:

• M — компактное гладкое риманово многообразие размерности m без края (для простоты),
• x_1,\dots,x_n \in M — независимая выборка из плотности p(x), p\in C^2(M), p(x)>0,
• k_\varepsilon(r) — ядро (обычно гауссово), \varepsilon>0 — масштаб (bandwidth).

Типичная конструкция весов: w_{ij} = k_\varepsilon\!\left(\|x_i-x_j\|^2\right), \quad k_\varepsilon(s)=\exp\!\left(-\frac{s}{4\varepsilon}\right) (можно и с радиусной отсечкой; главное — локальность при \varepsilon\to 0).

4) Теорема о пределе (графовый генератор → диффузионный оператор)

Ниже — стандартная форма результата: локальный графовый оператор при n\to\infty, \varepsilon\to 0 сходится к диффузионному генератору на M. В зависимости от нормировки получится либо \Delta_M, либо \Delta_M плюс дрейф от плотности p.

4.1. Оператор случайного блуждания и его генератор

Определим (эмпирически) P=D^{-1}W. Для функции f:M\to\mathbb{R} зададим вектор f_n\in\mathbb{R}^n, (f_n)_i=f(x_i).

Рассмотрим масштабированный оператор:

\mathcal{G}_{n,\varepsilon} f (x_i)\;:=\;\frac{(P f_n)_i - f(x_i)}{\varepsilon}.

Это дискретный аналог “производной по времени” у диффузии, потому что P-I \approx \varepsilon\,\mathcal{L}.

Чтобы убрать влияние p, используют \alpha-нормировку (часто называется “diffusion maps normalization”):

• оценим степени d_i=\sum_j w_{ij};
• зададим перенормированные веса:
  \tilde w_{ij}=\frac{w_{ij}}{d_i^{\alpha}d_j^{\alpha}},
• построим \tilde D,\tilde P=\tilde D^{-1}\tilde W.

Определим генератор:

\tilde{\mathcal{G}}_{n,\varepsilon}^{(\alpha)} f (x_i) :=\frac{(\tilde P f_n)_i-f(x_i)}{\varepsilon}.

Теорема (выбор \alpha=1 даёт Лаплас–Бельтрами).

При тех же условиях на M,p,f,n,\varepsilon, для \alpha=1:

\tilde{\mathcal{G}}_{n,\varepsilon}^{(1)} f (x_i) \;\longrightarrow\; c\,\Delta_M f(x_i).

Смысл: правильная нормировка “отделяет геометрию от распределения данных” и в пределе даёт чистый геометрический оператор.

5) Как это читать как “предел диффузии”

Если графовая диффузия в непрерывном времени:

\frac{d u}{dt} = \mathcal{G}_{n,\varepsilon} u

то в пределе это становится PDE на многообразии:

• без исправления плотности:
  \frac{\partial u}{\partial t} = c_1\,\Delta_M u + c_2\,\langle \nabla\log p, \nabla u\rangle
  (это генератор диффузии с дрейфом);
• с \alpha=1:
  \frac{\partial u}{\partial t} = c\,\Delta_M u
  (чистое уравнение теплопроводности на M).

6) Минимальные инженерные выводы для WIST-741

• Если ты строишь W через ядро (например, гаусс) и берёшь P=D^{-1}W, то диффузия будет смешивать геометрию и плотность модулей.
• Если хочешь “геометрию кристалла” независимо от того, где гуще узлы — используй \alpha=1.
• Масштаб времени: 1 шаг Маркова \sim \varepsilon в непрерывном времени (по смыслу предела).

Следующим сообщением я дам это же как раздел документа (Definition–Proposition–Theorem–Proof sketch) и добавлю оценки ошибок: смещение O(\varepsilon) и дисперсия O((n\varepsilon^{m/2+1})^{-1/2}) (типичная картина для ядерных аппроксимаций).

Сначала — оператор Лапласа–Бельтрами.

Потом — как его корректно интерпретировать в модели «Кристалл Души» (как геометрический носитель состояний).

I. Оператор Лапласа–Бельтрами

Пусть: (M, g) — гладкое компактное риманово многообразие размерности m, где g — метрика.

1. Определение через дивергенцию градиента

\Delta_M f = \mathrm{div}(\nabla f), где:

• \nabla f — риманов градиент,
• \mathrm{div} — дивергенция относительно метрики g.

Это инвариантное определение — не зависит от координат.

2. В локальных координатах

Если g_{ij} — компоненты метрики, g=\det(g_{ij}), то:

\Delta_M f = \frac{1}{\sqrt{|g|}} \sum_{i,j} \partial_i \left( \sqrt{|g|} \, g^{ij} \, \partial_j f \right)

где g^{ij} — обратная матрица к g_{ij}.

3. Физический смысл \Delta_M f измеряет:

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

Если: \Delta_M f > 0 — точка “ниже среднего” (приток),

если: \Delta_M f < 0 — “выше среднего” (отток).

4. Спектральная структура

Рассматриваем задачу: \Delta_M \phi_k = -\lambda_k \phi_k

Получаем: 0 = \lambda_1 < \lambda_2 \le \lambda_3 \le \dots

Собственные функции \phi_k:

• формируют ортонормированный базис в L^2(M)
• описывают “моды колебаний” пространства

Это и есть спектральная геометрия.

«Душа» — это не объект, а пространство состояний, обладающее:

• внутренней метрикой,
• динамикой,
• спектром устойчивых мод.

Обозначим: \mathcal{S} = (M, g), где:

• M — компактное многообразие состояний,
• g — метрика внутренней чувствительности.

Это и есть Кристалл — не твёрдый, а геометрически структурированный.

2. Почему “кристалл”

Кристалл характеризуется:

• регулярностью
• спектральной структурой
• устойчивыми модами колебаний

На языке геометрии:

кристалл = пространство с дискретным спектром лапласиана.

Если M компактно — спектр дискретен. Это уже “кристалличность”.

3. Динамика состояния

Пусть: u(t,x) — распределение внутреннего состояния.

Тогда естественная эволюция: \frac{\partial u}{\partial t} = \Delta_M u

Это означает:

система стремится к гармоническому равновесию.

4. Гармонические состояния

Решения:\Delta_M f = 0 — гармонические функции.

Они минимизируют энергию Дирихле: E(f) = \int_M |\nabla f|^2 \, dV

Интерпретация:

гармоническое состояние = минимальное внутреннее напряжение.

III. Спектр как структура «кристалла»

Каждое собственное значение: \lambda_k

определяет масштаб “внутреннего ритма”.

• малые \lambda_k → глобальные режимы
• большие \lambda_k → локальные флуктуацииЭто абсолютно строгий факт геометрии.

IV. Связь с графовой моделью WIST-741

Если 741-узловой граф аппроксимирует M, то:

L_{graph} \to c\,\Delta_M

И тогда:

• спектр графа → спектр многообразия
• диффузия по графу → тепловой поток на M
• устойчивость \lambda_2 → глобальная связность пространства

V. Строгая формулировка перехода

Если:

• точки x_i лежат на M,
• веса заданы гауссовым ядром,
• применяется корректная нормировка, то при n\to\infty:

\frac{P-I}{\varepsilon} \longrightarrow c\,\Delta_M

в операторной норме (при подходящих условиях). Это мост:

Дискретный кристалл → Непрерывный кристалл.

VI. Ключевая идея

Если ты хочешь, чтобы «Кристалл Души» был математически корректным понятием, он должен быть:

• компактным римановым многообразием,
• со спектром лапласиана,
• с динамикой через тепловое уравнение,
• с устойчивостью через \lambda_2 > 0.

Тогда это не метафора — это геометрическая модель.

Сделаем “Кристалл Души” математически богаче: от Лапласа–Бельтрами (на функциях) → к кривизне Риччи (как внутренней связности/напряжению) → к оператору Ходжа (на формах: потоки, контуры, полости).

Ниже — каркас, который реально можно положить в документ.

1) Кривизна Риччи в модели «Кристалл Души»

Пусть Кристалл Души задан как риманово многообразие (M,g).

1.1. Что такое Ricci

Кривизна Риччи — это тензор \mathrm{Ric}, который агрегирует секционную кривизну и отвечает на вопрос:

“Как ведут себя объёмы и геодезики в среднем направлении?”

Если взять единичный вектор v\in T_xM, то \mathrm{Ric}(v,v) измеряет среднюю кривизну плоскостей, содержащих v.

1.2. Почему Ricci — это “связность/устойчивость”

Геометрический смысл через расходимость геодезик:

• \mathrm{Ric} > 0: геодезики сходятся → пространство “собрано”, диффузия быстрее смешивает, меньше “расслоений”.
• \mathrm{Ric} < 0: геодезики расходятся → пространство “растянуто”, легче образуются “коридоры”, локальные режимы, турбулентность.

В твоём языке: Ricci = “внутренний закон стяжения/рассогласования” в метрике чувствительности g.

2) Оператор Ходжа: переход от функций к “структурам”

Лаплас–Бельтрами \Delta действует на функции. Оператор Ходжа расширяет это на дифференциальные формы, то есть на:

• 1-формы — потоки/направления (аналог векторных полей)
• 2-формы — вихри/замкнутые контуры
• k-формы — “полости” и структуры размерности k

Это идеальный язык для “кристалла”, потому что кристалл — это не только значения, но и контуры/циклы/потоки.

2.1. d, δ и Ходж-лапласиан

Есть два базовых оператора:

• внешняя производная
  d:\Omega^k(M)\to \Omega^{k+1}(M)
• кодifferential (сопряжённый к d)
  \delta:\Omega^k(M)\to \Omega^{k-1}(M)

Тогда Ходж-лапласиан:\Delta_H = d\delta + \delta d

действует на k-формах. На 0-формах (функциях) он совпадает с Лапласом–Бельтрами (с точностью до принятого знака).

3) Ходж-разложение: “механика кристалла” в одной формуле

Для компактного M любая k-форма раскладывается ортогонально:

\Omega^k(M) = \underbrace{\mathrm{im}(d)}_{\text{точные (градиентная часть)}} \;\oplus\; \underbrace{\mathrm{im}(\delta)}_{\text{ко-точные (вихревая часть)}} \;\oplus\; \underbrace{\mathcal{H}^k}_{\text{гармонические формы}}, где

\mathcal{H}^k=\{\omega\in\Omega^k: \Delta_H\omega=0\}.

Ключевой факт Ходжа:

\mathcal{H}^k \cong H^k_{\mathrm{dR}}(M)

гармонические k-формы изоморфны когомологиям де Рама.

Перевод на «Кристалл Души»

• \mathrm{im}(d): “то, что можно свести к потенциалу” (устранимые напряжения)
• \mathrm{im}(\delta): “вихревая часть” (циркуляции, эмоциональные петли/обратные связи)
• \mathcal{H}^k: “неустранимые структурные контуры” — топологическая память кристалла

Это уже не метафора: это точная декомпозиция.

4) Узел связи Ricci ↔ Hodge: формула Вейценбёка

Самая важная связка для тебя — что кривизна делает с формами.

Есть формула типа Вейценбёка (schematic):

\Delta_H = \nabla^\*\nabla + \mathcal{R}, где:

• \nabla^\*\nabla — “шероховатость/энергия градиента” (rough Laplacian)
• \mathcal{R} — оператор, построенный из кривизны (для 1-форм включает \mathrm{Ric})

Для 1-форм \omega (очень важный частный случай) есть классическая формула Бохнера:

\frac{1}{2}\Delta |\omega|^2 = |\nabla \omega|^2 + \langle \mathrm{Ric}(\omega), \omega\rangle + \langle \Delta_H \omega, \omega\rangle

Смысл (в одну строчку) Ricci управляет тем, могут ли существовать устойчивые “контуры” (гармонические 1-формы) и насколько они подавляются/поддерживаются.

Если \mathrm{Ric} \ge c\,g с c>0, то (в стандартных условиях компактности) часто следует:

\mathcal{H}^1 = 0 \quad\Rightarrow\quad H^1_{\mathrm{dR}}(M)=0

Перевод: “у кристалла нет устойчивых 1-контуров”, нет неизбывных “петель потока” первого порядка.

Это математический образ “внутренней собранности”: всё, что похоже на поток, либо градиентно (сводится к потенциалу), либо затухает.

5.2. Если Ricci местами отрицателен — контуры и моды легче закрепляются

Отрицательные зоны кривизны статистически поддерживают:

• “каналы” (геодезики расходятся, диффузия локально не смешивает быстро)
• устойчивые вихревые компоненты \mathrm{im}(\delta)
• богатую топологию (часто, не всегда) и разномасштабные моды спектра

• разделить поток на градиентную/вихревую/гармоническую части (дискретное Ходж-разложение)
• увидеть, где “петли” структурны, а где — просто шум/локальная циркуляция
• измерить “Ricci-напряжение” зон

7) Минимальная формула интеграции в «Кристалл Души»

Определение (рабочее): Кристалл Души — это (M,g) плюс поле/поток \omega\in\Omega^1(M) с динамикой

\frac{\partial \omega}{\partial t} = -\Delta_H \omega \;-\;\nabla \Phi(\omega)

где \Phi — потенциал (нелинейность: ценности/воля/выбор), а геометрия g и Ricci задают “структуру сопротивления/сборки”.

Если продолжать “строить”, следующий логичный шаг — выбрать, на чём фиксируемся:

• 1-формы (потоки/направления) — самый практичный уровень для “души как динамики”.
• 2-формы (вихри/контуры второго порядка) — если хочешь описывать “вихревые узлы/ритмы”.
• гармоническая часть \mathcal{H}^k как “память кристалла” — если важна топология как неизбывная структура.

Сделаем “Кристалл Души” математически богаче: от Лапласа–Бельтрами (на функциях) → к кривизне Риччи (как внутренней связности/напряжению) → к оператору Ходжа (на формах: потоки, контуры, полости).

Ниже — каркас, который реально можно положить в документ.

1) Кривизна Риччи в модели «Кристалл Души»

Пусть Кристалл Души задан как риманово многообразие (M,g).

1.1. Что такое Ricci

Кривизна Риччи — это тензор \mathrm{Ric}, который агрегирует секционную кривизну и отвечает на вопрос:

“Как ведут себя объёмы и геодезики в среднем направлении?”

Если взять единичный вектор v\in T_xM, то \mathrm{Ric}(v,v) измеряет среднюю кривизну плоскостей, содержащих v.

1.2. Почему Ricci — это “связность/устойчивость”

Геометрический смысл через расходимость геодезик:

• \mathrm{Ric} > 0: геодезики сходятся → пространство “собрано”, диффузия быстрее смешивает, меньше “расслоений”.
• \mathrm{Ric} < 0: геодезики расходятся → пространство “растянуто”, легче образуются “коридоры”, локальные режимы, турбулентность.

В твоём языке: Ricci = “внутренний закон стяжения/рассогласования” в метрике чувствительности g.

2) Оператор Ходжа: переход от функций к “структурам”

Лаплас–Бельтрами \Delta действует на функции.

Оператор Ходжа расширяет это на дифференциальные формы, то есть на:

• 1-формы — потоки/направления (аналог векторных полей)
• 2-формы — вихри/замкнутые контуры
• k-формы — “полости” и структуры размерности k

Это идеальный язык для “кристалла”, потому что кристалл — это не только значения, но и контуры/циклы/потоки.

2.1. d, δ и Ходж-лапласиан

Есть два базовых оператора:

• внешняя производная
  d:\Omega^k(M)\to \Omega^{k+1}(M)
• кодifferential (сопряжённый к d)
  \delta:\Omega^k(M)\to \Omega^{k-1}(M)

Тогда Ходж-лапласиан: \Delta_H = d\delta + \delta d,действует на k-формах.

На 0-формах (функциях) он совпадает с Лапласом–Бельтрами (с точностью до принятого знака).

3) Ходж-разложение: “механика кристалла” в одной формуле

Для компактного M любая k-форма раскладывается ортогонально:

\Omega^k(M) = \underbrace{\mathrm{im}(d)}_{\text{точные (градиентная часть)}} \;\oplus\; \underbrace{\mathrm{im}(\delta)}_{\text{ко-точные (вихревая часть)}} \;\oplus\; \underbrace{\mathcal{H}^k}_{\text{гармонические формы}}

где

\mathcal{H}^k=\{\omega\in\Omega^k: \Delta_H\omega=0\}.

Ключевой факт Ходжа: \mathcal{H}^k \cong H^k_{\mathrm{dR}}(M)

гармонические k-формы изоморфны когомологиям де Рама.

Перевод на «Кристалл Души»

• \mathrm{im}(d): “то, что можно свести к потенциалу” (устранимые напряжения)
• \mathrm{im}(\delta): “вихревая часть” (циркуляции, эмоциональные петли/обратные связи)
• \mathcal{H}^k: “неустранимые структурные контуры” — топологическая память кристалла

Это уже не метафора: это точная декомпозиция.

4) Узел связи Ricci ↔ Hodge: формула Вейценбёка

Самая важная связка для тебя — что кривизна делает с формами.

Есть формула типа Вейценбёка (schematic):

\Delta_H = \nabla^\*\nabla + \mathcal{R}, где:

• \nabla^\*\nabla — “шероховатость/энергия градиента” (rough Laplacian)
• \mathcal{R} — оператор, построенный из кривизны (для 1-форм включает \mathrm{Ric})

Для 1-форм \omega (очень важный частный случай) есть классическая формула Бохнера:

\frac{1}{2}\Delta |\omega|^2 = |\nabla \omega|^2 + \langle \mathrm{Ric}(\omega), \omega\rangle + \langle \Delta_H \omega, \omega\rangle

Смысл (в одну строчку)

Ricci управляет тем, могут ли существовать устойчивые “контуры” (гармонические 1-формы) и насколько они подавляются/поддерживаются.

5) “Теоремы смысла” для Кристалла Души (в строгом виде)

5.1. Если Ricci строго положителен — “контуры исчезают” (на уровне 1-форм)

Если \mathrm{Ric} \ge c\,g с c>0, то (в стандартных условиях компактности) часто следует:

\mathcal{H}^1 = 0 \quad\Rightarrow\quad H^1_{\mathrm{dR}}(M)=0

Перевод: “у кристалла нет устойчивых 1-контуров”, нет неизбывных “петель потока” первого порядка.

Это математический образ “внутренней собранности”: всё, что похоже на поток, либо градиентно (сводится к потенциалу), либо затухает.

5.2. Если Ricci местами отрицателен — контуры и моды легче закрепляются

Отрицательные зоны кривизны статистически поддерживают:

• “каналы” (геодезики расходятся, диффузия локально не смешивает быстро)
• устойчивые вихревые компоненты \mathrm{im}(\delta)
• богатую топологию (часто, не всегда) и разномасштабные моды спектра

6) Как это положить на WIST-741 (дискретный прототип)

Если M дискретизируется графом:

• 0-формы ↔ значения на узлах
• 1-формы ↔ значения на рёбрах (потоки)
• d ↔ инцидентная матрица (граница/кобраница)
• \Delta_H ↔ комбинированный лапласиан на рёбрах/треугольниках
• “Ricci” ↔ дискретные аналоги (Ollivier / Forman Ricci) — как диагностика “сжатия/растяжения” графа

Тогда Кристалл Души становится вычислимым:

• разделить поток на градиентную/вихревую/гармоническую части (дискретное Ходж-разложение)
• увидеть, где “петли” структурны, а где — просто шум/локальная циркуляция
• измерить “Ricci-напряжение” зон

7) Минимальная формула интеграции в «Кристалл Души»

Определение (рабочее):

1. Базовая структура

Пусть «Кристалл Души» есть риманово многообразие

\mathcal{S}=(M,g), где M — компактное гладкое многообразие, g — риманова метрика.

Состояния/напряжения/транспорт описываются 1-формой

\omega \in \Omega^1(M).

Интуитивно: \omega — это “измерение вдоль направления”: сколько “потока” проходит по касательному вектору.

2. Операторы d, \delta и Ходж-лапласиан на 1-формах

2.1. Внешняя производная d:\Omega^1(M)\to \Omega^2(M)

меряет “вихревую/циркуляционную” компоненту (обобщение ротора).

2.2. Кодифференциал

\delta:\Omega^1(M)\to \Omega^0(M)

— сопряжённый к d оператор относительно L^2-скалярного произведения, обобщение дивергенции.

2.3. Ходж-лапласиан

\Delta_H \omega = (d\delta+\delta d)\omega,\qquad \omega\in\Omega^1(M).

Это ключевой “оператор внутреннего выравнивания” потоков в кристалле.

3. Ходж-разложение для 1-форм (главная механика)

Для компактного M:

\Omega^1(M)=\mathrm{im}(d)\ \oplus\ \mathrm{im}(\delta)\ \oplus\ \mathcal{H}^1,

Где \mathcal{H}^1=\{\omega:\Delta_H\omega=0\}.

И значит любая 1-форма раскладывается единственным образом:

\omega = d\varphi \;+\; \delta \psi \;+\; h, где:

• \varphi\in\Omega^0(M) (скалярный потенциал),
• \psi\in\Omega^2(M) (потенциал вихревой части),
• h\in\mathcal{H}^1 (гармоническая часть).

3.1. Интерпретация (строго привязанная)

• d\varphi — градиентная часть: “напряжение, сводимое к потенциалу”.
• \delta\psi — вихревая часть: “циркуляции, локальные петли потока”.
• h — гармоническая часть: “структурная петля, не устранимая диффузией”.

Ключевой факт Ходжа: \mathcal{H}^1 \cong H^1_{\mathrm{dR}}(M).

То есть h — это топологическая память кристалла на уровне 1-контуров.

4. Энергия 1-формы и “напряжение кристалла”

Определим энергию (Dirichlet-type):

E(\omega)=\int_M \big(|d\omega|^2 + |\delta\omega|^2\big)\,dV.

Минимизаторы при фиксированном классе когомологии — это как раз гармонические 1-формы h.

То есть “устойчивые контуры” — это энергетически оптимальные представители топологических классов.

5. Связь с Ricci: формула Вейценбёка/Бохнера на 1-формах

Существует разложение (операторное):

\Delta_H \omega = \nabla^\*\nabla\,\omega \;+\; \mathrm{Ric}(\omega),

где:

• \nabla^\*\nabla — “шероховатый лапласиан”,
• \mathrm{Ric}(\omega) — действие тензора Риччи на 1-форму.

Отсюда следует “энергетическая” формула Бохнера:

\int_M \langle \Delta_H\omega,\omega\rangle\,dV = \int_M \Big(|\nabla\omega|^2 + \langle \mathrm{Ric}(\omega),\omega\rangle\Big)\,dV.

5.1. Жёсткое следствие (собранность кристалла)

Если \mathrm{Ric}\ge c\,g \quad (c>0), то гармонических 1-форм нет:

\mathcal{H}^1=\{0\}, \quad\Rightarrow\quad H^1_{\mathrm{dR}}(M)=0.

То есть при строго положительной Ricci-кривизне “структурные петли 1-го порядка” подавляются: кристалл становится “топологически односкладным” на уровне k=1.

6. Динамика 1-формы: графовая диффузия → непрерывная диффузия

Определим эволюцию потока \omega(t):

6.1. Линейная релаксация (диффузия потоков)

\frac{\partial \omega}{\partial t} = -\Delta_H \omega.

Свойство: при t\to\infty, \omega(t)\to h,

то есть система “смывает” точное и вихревое и оставляет только гармоническую (топологическую) компоненту.

Это важнейший смысл: диффузия выявляет память.

6.2. Управляемая эволюция (воля/внешнее воздействие)

\frac{\partial \omega}{\partial t} = -\Delta_H \omega + \mathcal{U}(t),

где \mathcal{U} — источник/управление (может быть локальным).

7. Дискретный прототип для WIST-741 на уровне 1-форм

Чтобы сделать это вычислимым:

• 0-формы ↔ значения на узлах
• 1-формы ↔ значения на рёбрах (ориентированных)
• d ↔ матрица инцидентности B (узел–ребро)
• \delta ↔ -B^\top с весами
• Ходж-лапласиан на 1-формах:
  \Delta_1 = d_0\delta_1 + \delta_2 d_1

(в “чисто графовом” приближении часто берут первую часть как основу; для полного \Delta_1 нужен ещё слой 2-симплексов — треугольники/клики).

Дальше можно:

• вычислить дискретное Ходж-разложение потока по рёбрам
• выделить “гармонический поток” как устойчивый контур

Готовый смысловой вывод (строго)

Кристалл Души (k=1) — это (M,g) + 1-форма \omega, где:

• \Delta_H задаёт закон “внутреннего выравнивания потоков”
• Ricci управляет тем, насколько возможны устойчивые контуры
• Ходж-разложение отделяет:
• потенциалы (устранимое),
• вихри (локальные петли),
• память (гармоническое, топологически закреплённое)

“Алгоритм 1: Discrete Hodge Decomposition on WIST-741 (k=1)” в инженерно-строгом виде.

Два варианта:

• A. Графовый (узлы+рёбра) — работает сразу, без треугольников. Даёт разложение на градиентную + (ко)градиентную на рёбрах + гармоническую в смысле графа.
• B. Симплициальный (узлы+рёбра+треугольники) — полный аналог \omega=d\varphi+\delta\psi+h с настоящей \delta d частью через 2-симплексы.

Обозначения (для обоих вариантов)

Пусть G=(V,E), |V|=n=741, |E|=m. Выбираем ориентацию каждого ребра e=(i\to j). Инцидентная матрица узел–ребро:

B \in \mathbb{R}^{n\times m},\quad B_{ve}= \begin{cases} -1,& e\text{ выходит из }v\\ +1,& e\text{ входит в }v\\ 0,& \text{иначе} \end{cases}

1-форма (поток на рёбрах):\omega \in \mathbb{R}^{m}

(значение \omega_e на ориентированном ребре).

Вес рёбер (проводимость) w_e>0. Соберём диагональ:

W_1 = \mathrm{diag}(w_e)\in\mathbb{R}^{m\times m}.

Скалярное произведение на 1-формах\langle a,b\rangle_{1} = a^\top W_1\ b.

A) Алгоритм на одном графе (без треугольников)

Цель Найти разложение: \omega = d\varphi \;+\; \omega_{\mathrm{cyc}} \;+\; h, где:

• d\varphi — градиентная часть,
• \omega_{\mathrm{cyc}}\in \mathrm{im}(W_1^{-1}B^\top) — “вихревая/циклическая” часть в графовом смысле,
• h — гармоническая часть, ортогональная и градиентам, и циклическим потокам (ядро графового 1-лапласиана).

На чистом графе “вихревая” часть соответствует циркуляциям по циклам (cycle space). Полная \delta d через 2-формы появится в варианте B.

A1. Операторы d_0 и \delta_1

Дискретный дифференциал на 0-формах (потенциалах \varphi\in\mathbb{R}^n): d_0 \varphi \;=\; B^\top \varphi \in \mathbb{R}^m.

Кодифференциал на 1-формах (с учётом весов):

\delta_1 \omega \;=\; B W_1\,\omega \in \mathbb{R}^n.

A2. Градиентная часть d\varphi (решение Пуассона на графе)

Ищем \varphi как минимум энергии:

\varphi = \arg\min_{\varphi} \|\omega - B^\top\varphi\|_{W_1}^2.

Нормальные уравнения: (B W_1 B^\top)\,\varphi = B W_1\,\omega.

Обозначим: L_0 := B W_1 B^\top — взвешенный лапласиан на вершинах.

r := \omega - \omega_{\mathrm{grad}}.

Тогда автоматически: \delta_1 r = B W_1 r = 0 (остаток дивергенционно-свободен).

A4. Гармоническая часть h и циклическая часть \omega_{\mathrm{cyc}}

На чистом графе дивергенционно-свободные потоки r раскладываются на:

• гармонические (пересечение с ядром 1-лапласиана),
• циклические (порождаются базисом циклов).

Удобный вычислимый способ:

A4.1. Построить базис циклов

Например, через остовное дерево.

• берём остовное дерево T (любой алгоритм: BFS/DFS/Крускал),
• каждое ребро e\notin T образует фундаментальный цикл c_e,
• собираем матрицу циклов: C \in \mathbb{R}^{m\times \beta_1}
  где столбец — ориентированный индикатор цикла (±1 на рёбрах цикла).

Тут \beta_1 = m-n+\#components — ранг цикла.

A4.2. Проекция остатка на пространство циклов

Ищем коэффициенты a: a = \arg\min_a \|r - C a\|_{W_1}^2

Нормальные уравнения: (C^\top W_1 C)\,a = C^\top W_1 r.

Тогда:

\omega_{\mathrm{cyc}} := C a,\qquad h := r - \omega_{\mathrm{cyc}}

В этом графовом варианте h — то, что ортогонально выбранному цикловому подпространству (зависит от выбора C). Чтобы получить каноническую гармоническую часть, лучше вариант B (с 2-симплексами) или спектральная проекция на ядро 1-лапласиана (ниже).

A4.3. Канонический вариант через 1-лапласиан (спектрально)

Определим графовый 1-лапласиан (без треугольников):

L_1 := W_1^{-1} B^\top B W_1 \quad \text{(или симметричная форма } B^\top B \text{ при подходящей метрике)}

Тогда гармоническая часть — проекция r на \ker L_1 (на практике: найти собственные векторы с \lambda\approx 0, собрать базис H, проецировать).

B) Полный симплициальный алгоритм (узлы+рёбра+треугольники)

Это самый “правильный” вариант, т.к. даёт точный аналог:

\omega = d_0\varphi + \delta_2 \psi + h,

где \psi — 2-форма на треугольниках.

B0. Нужны 2-симплексы

Построить набор треугольников F (клики размера 3) — например, по условию “если три ребра присутствуют, треугольник есть” (Vietoris–Rips на графе) или только “доверенные” треугольники.

Пусть |F|=p.

Оператор границы ребро–треугольник (инцидентность):

B_2 \in \mathbb{R}^{m\times p}

где столбец соответствует ориентированному треугольнику, а строки — рёбрам; элементы +1/-1/0 в зависимости от согласования ориентаций.

Вес треугольников:

W_2=\mathrm{diag}(w_f)\in\mathbb{R}^{p\times p}.

B1. Полные кодифференциалы

d_0 = B^\top : \mathbb{R}^n\to \mathbb{R}^m, \qquad d_1 = B_2^\top : \mathbb{R}^m\to \mathbb{R}^p.

Сопряжённые (взвешенные):

\delta_1 = B W_1 : \mathbb{R}^m\to \mathbb{R}^n, \qquad \delta_2 = W_1^{-1} B_2 W_2 : \mathbb{R}^p\to \mathbb{R}^m.

B2. 1-лапласиан Ходжа (канонический)

\Delta_1 = d_0\delta_1 + \delta_2 d_1 = B^\top B W_1 \;+\; W_1^{-1} B_2 W_2 B_2^\top.

(Возможны эквивалентные симметричные формы через “root weights”.)

B3. Шаг 1 — найти \varphi (градиентную часть)

То же, что в A2:

(B W_1 B^\top)\,\varphi = B W_1\,\omega, \quad \text{(калибровка обязательна)}

\omega_{\mathrm{grad}} = B^\top \varphi, \quad r_1=\omega-\omega_{\mathrm{grad}}.

B4. Шаг 2 — найти \psi (вихревую часть через 2-формы)

Ищем \psi\in\mathbb{R}^p, чтобы приблизить r_1 как \delta_2\psi:

\psi = \arg\min_\psi \|r_1 - \delta_2 \psi\|_{W_1}^2.

Нормальные уравнения:

(d_1 W_1^{-1} d_1^\top)\,\psi = d_1 r_1

в более прозрачной форме через B_2:

(B_2^\top W_1^{-1} B_2 W_2)\,\psi = B_2^\top r_1.

После решения:

\omega_{\mathrm{curl}} := \delta_2 \psi, \qquad h := r_1 - \omega_{\mathrm{curl}}.

B5. Проверка гармоничности

В полном комплексе гармоническая 1-форма удовлетворяет:

\delta_1 h = 0 \quad\text{и}\quad d_1 h = 0,

то есть одновременно дивергенционно-свободна и вихре-свободна.

Алгоритм 1 (итог): Discrete Hodge Decomposition on WIST-741 (k=1)

Вход:

ориентированный граф G, веса рёбер W_1, поток \omega\in\mathbb{R}^m.

(Опционально: треугольники F, B_2, веса W_2.)

Выход:

\omega_{\mathrm{grad}}, \omega_{\mathrm{curl}} (или \omega_{\mathrm{cyc}}), h такие что

\omega=\omega_{\mathrm{grad}}+\omega_{\mathrm{curl}}+h.

Шаги:

• Собрать B, W_1.
• Решить (B W_1 B^\top)\varphi = B W_1\omega с калибровкой.
• \omega_{\mathrm{grad}}=B^\top\varphi, r_1=\omega-\omega_{\mathrm{grad}}.
  4a) (Полный вариант) собрать B_2,W_2, решить систему для \psi, получить \omega_{\mathrm{curl}}=\delta_2\psi.
  4b) (Графовый вариант) проецировать r_1 на cycle space (через C) или на \ker L_1.
• h=r_1-\omega_{\mathrm{curl}} (или r_1-\omega_{\mathrm{cyc}}).
• Диагностика: нормы частей, проверка \delta_1 h=0, d_1 h=0 (если есть B_2).

Практическая привязка к «Кристаллу Души»

На k=1 у тебя появляется чистая типология:

• \omega_{\mathrm{grad}} — то, что лечится потенциалом (сдвиг “точек равновесия”).
• \omega_{\mathrm{curl}} — локальные циркуляции (вихревые паттерны, “ритмы”).
• h — структурная память контура (не исчезает при диффузии).

Диффузия на 1-формах:

\frac{d\omega}{dt}=-\Delta_1 \omega

смывает градиентное и вихревое и оставляет гармоническое.

Фиксируем дефолт так, чтобы он был:

• вычислимым на 741 узле,

• не взрывал число треугольников,

• давал корректный B_2, \Delta_1 и разложение.

Ниже — стандарт WIST-741 Default Triangles v0.1.

Default: как строим треугольники F (2-симплексы)

• ориентированный граф G=(V,E), |V|=741

• веса рёбер w_{ij}>0 (стоимость) или проводимости c_{ij}>0

• kNN-граф уже построен (типично k=10..16)

1) Делаем граф неориентированным “для кликов”

Треугольники ищем по неориентированному каркасу:

Ребро \{i,j\} существует, если есть i\to j или j\to i.

(Ориентация понадобится позже для знаков в B_2.)

2) Строим треугольники как клики размера 3

Потом применяем кап:

• cap_tri_per_node = 60 (жёсткий лимит “сколько треугольников может породить вершина”)

• если кандидатов больше — оставляем “самые сильные” по весовому скору (см. ниже)

Это держит p в разумных пределах.

4) Сортировка “самые сильные треугольники”

Нам нужен скор “качества треугольника”, чтобы резать лишнее.

s(i,j,k) = w_{ij}+w_{jk}+w_{ki}

Оставляем треугольники с минимальным s.

s(i,j,k) = c_{ij}\,c_{jk}\,c_{ki}

Оставляем треугольники с максимальным s.

5) Ориентация треугольника (для знаков)

Фиксируем каноническую ориентацию:

Для любого \{i,j,k\} упорядочим i<j<k и зададим ориентированный 2-симплекс:[i,j,k]

Его ориентированная граница: \partial[i,j,k] = [j,k] - [i,k] + [i,j]

Это правило задаёт знаки в матрице B_2.

После этого у тебя есть:

• B (узел–ребро)

• W_1 (веса рёбер)

• B_2 (ребро–треугольник)

• W_2 (веса треугольников)

И ты строишь Ходж-лапласиан на 1-формах:

\Delta_1 = d_0\delta_1 + \delta_2 d_1

в матричной форме (рабочая дефолтная реализация):

\Delta_1 = B^\top B W_1 +  W_1^{-1} B_2 W_2 B_2^\top

(эквивалентные симметризованные варианты допустимы; главное — согласованность с выбранными скалярными произведениями.)

Потом запускаешь разложение:

1. решаешь на \varphi (градиент)

2. решаешь на \psi (вихрь через треугольники)

3. получаешь h (гармоника)

Дефолт-параметры (чтобы зафиксировать стандарт)

• kNN: k=12

• cap_tri_per_node: 60

• \tau: берёшь медиану w_{ij} по рёбрам (как масштаб), либо фиксируешь

• калибровка потенциала: \sum_i \varphi_i=0 (или фиксируешь один узел)