"Intelligence is the optimal balance between memorizing the past and compressing the future."
Status: Currently scaling a sub-quadratic continuous state-space model across an arbitrary number of interconnected GPUs. Expect severe delays in correspondence due to global gradient synchronization.
The current discourse in Artificial Intelligence is highly fractured: autoregressive models scale memorization, while pure reinforcement learning struggles with sample efficiency. My research bridges these domains by framing intelligence through the lens of Active Inference, Geometric Deep Learning, and Continuous State-Space Models (SSMs).
Let
- Geometric Deep Learning: Exploiting symmetry and invariance (group theory) in neural architectures to learn efficiently from non-Euclidean manifolds.
- Continuous State-Space Models: Pushing beyond standard attention mechanisms using structured sequence models (e.g., Mamba, S4) for unbounded context horizons.
- Active Inference: Designing agents that don't just passively predict the world, but take actions to minimize the divergence between their predictions and sensory reality.
-
Applied Coffee Optimization: Theorem 1 states there exists an optimal learning rate
$\eta \in (0, 1]$ such that a model converges precisely when my coffee mug empties. The proof is trivial and left as an exercise to the GPU cluster.
Intelligence fundamentally requires measuring the compatibility between an observation
Let
By treating an action
Note: If $\mathcal{F}$ diverges to infinity, it is scientifically customary to blame the batch size.
[Expand] Proof of Asymptotic Global Convergence for Active Inference Policies
Theorem: Let an autonomous agent be governed by a continuous-time stochastic differential equation (SDE) over a smooth Riemannian manifold
Proof:
Let the internal state dynamics be described by the Itô SDE:
ds_t = f(s_t, a_t) dt + \Sigma(s_t) dW_t
where
\frac{\partial p}{\partial t} = -\nabla \cdot (f(s,a)p) + \frac{1}{2} \sum_{i,j} \frac{\partial^2}{\partial s_i \partial s_j} (\Sigma \Sigma^T)_{ij} p
The agent's objective is to minimize the expected free energy path integral over an infinite horizon
\mathcal{F}(\pi) = \mathbb{E}_{Q(o, s | \pi)} \left[ \int_0^\infty e^{-\gamma t} \left( \ln Q(s_t) - \ln P(o_t, s_t) \right) dt \right]
To find the optimal control policy
\gamma V(s) = \min_{a} \left\{ \mathcal{F}(s, a) + \nabla_s V(s)^T f(s,a) + \frac{1}{2} \text{Tr}\left( \Sigma(s)\Sigma(s)^T \nabla_{ss}^2 V(s) \right) \right\}
By taking the functional derivative with respect to the action
\frac{\partial \mathcal{F}}{\partial a} + \left( \frac{\partial f}{\partial a} \right)^T \nabla_s V(s) = 0
To prove global asymptotic stability, we propose
Taking the orbital derivative along the system trajectories:
\dot{V}(s_t) = \nabla_s V^T \dot{s}_t = \nabla_s V^T f(s_t, a^*_t) \leq -\mathcal{F}(s_t, a^*_t) < 0
Because
Thus, the agent perfectly predicts and dictates its environment, securing indefinite survival.
[2608.09112]Bounding Variational Free Energy in Continuous State-Space Models (Under Review)[2511.03450]Autoregressive Collapse: Why Next-Token Prediction Cannot Yield General Intelligence[2502.11899]Non-Euclidean Manifold Traversals in Hierarchical JEPAs
To support my open-source work and fund my compute clusters, I package my internal production tools and deep-dive technical notes into high-value bundles. If you find my research helpful, check them out:
- 🧠 Advanced AI Agent Architecture: Learn true LLM orchestration with ReAct loops and state management.
- 🌐 Mastering Distributed Systems in Go: Deep technical dives into Raft, Gossip, and consensus algorithms.
- ☁️ The Production AWS Infrastructure Boilerplate: Production-grade Terraform VPC + EKS boilerplate for rapid deployments.
- ⚛️ The Enterprise Next.js Boilerplate: Optimized React/Next.js 14 template for quick startup bootstrapping.
- 🗄️ The Data Engineering Airflow Toolkit: Local Airflow + Postgres Docker setup with production ETL DAGs.
- 🛡️ The Cybersecurity Script Bundle: Async Python port scanners and web log brute-force analyzers.
- 🤖 The Machine Learning Engineer's Starter Kit: PyTorch AMP training and FastAPI inference boilerplate.
- ⚡ The Ultimate Developer Productivity Pack: High-efficiency Bash aliases, Git hooks, and CI/CD pipelines.
For those seeking to escape the local minima of standard Deep Learning tutorials, I recommend the following foundational texts:
- Active Inference: The Free-Energy Principle: A Unified Brain Theory? — Karl Friston (2010)
- Geometric DL: Geometric Deep Learning: Grids, Groups, Graphs, Geodesics, and Gauges — Bronstein et al. (2021)
- State-Space Models: Mamba: Linear-Time Sequence Modeling with Selective State Spaces — Gu & Dao (2023)
Scaling these dynamic systems requires a robust engineering foundation to handle high-dimensional manifolds without bottlenecking.
Lemma 1: The Attention Bottleneck Standard self-attention scales at
$\mathcal{O}(N^2)$ with respect to sequence length$N$ . Without sub-quadratic architectures or continuous state-spaces, attempting infinite-context reasoning is isomorphic to heating the earth.
On Neuro-Symbolic Integration: Deep learning is unparalleled at statistical pattern matching (System 1), but struggles with rigid logical deduction (System 2). The next generation of models must embed symbolic constraints directly into the differentiable loss landscape.
- Professional Network: LinkedIn
- Code & Implementations: See public repositories below.
Hyper-Optimized Neural SDE Implementation (Variational Free Energy Predictor):
import torch
import torch.nn as nn
from torch.nn import functional as F
from typing import Tuple
@torch.compile(mode="reduce-overhead")
class NeuralSDEPredictor(nn.Module):
"""
Continuous-time latent state predictor modeling the Itô SDE:
d(s_t) = f_θ(s_t, a_t)dt + g_φ(s_t)dW_t
Optimized for heavily batched, non-Euclidean manifold traversals via
Lie group regularized integrators.
"""
def __init__(self, d_model: int = 4096, d_action: int = 1024):
super().__init__()
# Drift network (f_θ) with SwiGLU activations and RMSNorm
self.drift_proj = nn.Linear(d_model + d_action, d_model * 2, bias=False)
self.drift_out = nn.Linear(d_model, d_model, bias=False)
self.norm = nn.RMSNorm(d_model * 2)
# Diffusion network (g_φ) modeling irreducible aleatoric uncertainty
self.log_diffusion = nn.Parameter(torch.zeros(d_model))
# Spectral normalization to enforce Lipschitz continuity for HJB stability
nn.utils.parametrizations.spectral_norm(self.drift_out)
def drift(self, s_t: torch.Tensor, a_t: torch.Tensor) -> torch.Tensor:
"""Computes the deterministic drift vector field."""
x = torch.cat([s_t, a_t], dim=-1)
x = self.norm(self.drift_proj(x))
x, gate = x.chunk(2, dim=-1) # SwiGLU gating
return self.drift_out(x * F.silu(gate))
def forward(self, s_t: torch.Tensor, a_t: torch.Tensor, dt: float = 1e-3) -> Tuple[torch.Tensor, torch.Tensor]:
"""
Euler-Maruyama integration step with auxiliary penalty
for variational free energy minimization.
"""
f_t = self.drift(s_t, a_t)
g_t = torch.exp(self.log_diffusion)
dW = torch.randn_like(s_t) * (dt ** 0.5) # Wiener process increments
s_next = s_t + (f_t * dt) + (g_t * dW)
# Regularization penalty derived from the Fokker-Planck density evolution
penalty = 0.5 * torch.sum(g_t ** 2, dim=-1).mean()
return s_next, penalty


