SolitaryThinker
db866dfb04
Sweep all 296 v2 modules: simplify verbose docstrings/comments and remove 376 dangling "(design_vN §X)" citations to the now-deleted design docs (v2/README.md is the source of truth). Comment/docstring-only — AST-verified code-identical; the CPU suite holds at 237 passed / 5 pre-existing. Also applies yapf + ruff --fix auto-fixes (import ordering, forward-ref annotation de-quoting under `from __future__ import annotations`; behavior- neutral, suite-confirmed) and adds the legitimate domain terms mot/clen/te to the codespell ignore-list. Remaining ruff (24) + mypy (68 no-untyped-call) findings are pre-existing v2 debt, untouched here.
111 lines
5.7 KiB
Python
111 lines
5.7 KiB
Python
"""Samplers as a small library; model families compose them.
|
||
|
||
Flow-match Euler is shared by the Wan and LTX-2 denoise step bodies. On GPU the real UniPC
|
||
multistep / distilled schedulers are wrapped by the torch component adapters; this numpy form
|
||
is the CPU-testable, bit-reproducible reference used by the loop tests.
|
||
|
||
Flow-match interpolation: ``x_t = (1-σ_t)·x0 + σ_t·ε``; the model predicts velocity
|
||
``v = ε - x0``. A deterministic Euler step to σ_next is ``x_next = x_t + (σ_next-σ_t)·v``
|
||
(exactly LTX-2's ``velocity=(latents-denoised)/σ; latents += velocity·dt``).
|
||
"""
|
||
from __future__ import annotations
|
||
|
||
import numpy as np
|
||
|
||
|
||
def build_flow_sigmas(num_steps: int, shift: float = 1.0, terminal: float = 0.0) -> np.ndarray:
|
||
"""σ schedule from 1→terminal over ``num_steps+1`` points, with flow-shift applied.
|
||
|
||
Flow-shift: ``σ' = shift·σ / (1 + (shift-1)·σ)`` (Wan/LTX). shift=1 is the identity.
|
||
"""
|
||
base = np.linspace(1.0, terminal, num_steps + 1, dtype=np.float64)
|
||
if shift != 1.0:
|
||
base = shift * base / (1.0 + (shift - 1.0) * base)
|
||
return base
|
||
|
||
|
||
def x0_from_velocity(x_t: np.ndarray, velocity: np.ndarray, sigma_t: float) -> np.ndarray:
|
||
"""x0 = x_t - σ_t·v (flow prediction → clean sample; Wan ``x0 = x_t - σ·model_output``)."""
|
||
return x_t - sigma_t * velocity
|
||
|
||
|
||
def build_karras_sigmas(num_steps: int,
|
||
sigma_max: float = 80.0,
|
||
sigma_min: float = 0.002,
|
||
rho: float = 7.0) -> np.ndarray:
|
||
"""Karras et al. (2022) ρ-interpolated σ schedule + the terminal σ, as Cosmos/EDM models use it.
|
||
|
||
Reproduces ``FlowMatchEulerDiscreteScheduler(use_karras_sigmas=True)`` the Cosmos pipeline configures
|
||
(``sigma_max=80, sigma_min=0.002, rho=7``): a length-``num_steps`` ramp ``σ_max→σ_min`` via
|
||
``σ_i = (max_inv_rho + i/(n-1)·(min_inv_rho-max_inv_rho))^ρ``, then ONE appended terminal value. The
|
||
scheduler appends ``0.0`` and, with ``final_sigmas_type='sigma_min'``, the Cosmos stage overwrites it
|
||
with ``σ[-2]`` to avoid a divide-by-zero in the EDM coeffs / velocity — so the terminal here is
|
||
``σ_min`` (the last ramp value), giving ``num_steps+1`` points the EDM loop integrates pairwise.
|
||
"""
|
||
ramp = np.linspace(0.0, 1.0, num_steps, dtype=np.float64)
|
||
min_inv_rho = sigma_min**(1.0 / rho)
|
||
max_inv_rho = sigma_max**(1.0 / rho)
|
||
sigmas = (max_inv_rho + ramp * (min_inv_rho - max_inv_rho))**rho # σ_max -> σ_min
|
||
return np.concatenate([sigmas, sigmas[-1:]]).astype(np.float64) # + sigma_min terminal clamp
|
||
|
||
|
||
def flow_match_euler_step(x_t: np.ndarray, velocity: np.ndarray, sigma_t: float, sigma_next: float) -> np.ndarray:
|
||
"""One deterministic Euler step along the flow ODE."""
|
||
return x_t + (sigma_next - sigma_t) * velocity
|
||
|
||
|
||
def add_noise(x0: np.ndarray, noise: np.ndarray, sigma: float) -> np.ndarray:
|
||
"""Forward flow-match interpolation: x_σ = (1-σ)·x0 + σ·noise."""
|
||
return (1.0 - sigma) * x0 + sigma * noise
|
||
|
||
|
||
def flow_sde_step_with_logprob(x_t,
|
||
velocity,
|
||
sigma_t: float,
|
||
sigma_next: float,
|
||
*,
|
||
noise=None,
|
||
prev_sample=None,
|
||
noise_scale: float = 0.7):
|
||
"""FlowGRPO-style stochastic step + Gaussian log-prob.
|
||
|
||
The RL *rollout* sampler (vs the deterministic ODE ``flow_match_euler_step`` used at serve time).
|
||
Injects noise for exploration and returns the log-prob of the realized sample under the step's
|
||
Gaussian — the quantity the FlowGRPO PPO ratio uses. ``prev_sample=None`` samples a new step;
|
||
passing ``prev_sample`` (the rollout sample) recomputes its log-prob under the *current* velocity
|
||
(the ratio's numerator at update time). Returns ``(prev_sample, log_prob, mean, eff_std)``.
|
||
"""
|
||
x_t = np.asarray(x_t, dtype=np.float64)
|
||
velocity = np.asarray(velocity, dtype=np.float64)
|
||
s = min(float(sigma_t), 0.9999)
|
||
dt = float(sigma_next) - float(sigma_t) # negative (σ decreases)
|
||
std = float(np.sqrt(s / (1.0 - s)) * noise_scale)
|
||
denom = 2.0 * max(s, 1e-6)
|
||
mean = x_t * (1.0 + std**2 / denom * dt) + velocity * (1.0 + std**2 * (1.0 - s) / denom) * dt
|
||
eff_std = max(std * np.sqrt(max(-dt, 1e-12)), 1e-6)
|
||
if prev_sample is None:
|
||
n = noise if noise is not None else np.zeros_like(x_t)
|
||
prev_sample = mean + eff_std * np.asarray(n, dtype=np.float64)
|
||
prev_sample = np.asarray(prev_sample, dtype=np.float64)
|
||
var = eff_std**2
|
||
log_prob = float(np.mean(-((prev_sample - mean)**2) / (2.0 * var) - np.log(eff_std) - 0.5 * np.log(2.0 * np.pi)))
|
||
return prev_sample.astype("float32"), log_prob, mean.astype("float32"), float(eff_std)
|
||
|
||
|
||
def flow_sde_ml_velocity(x_t, sample, sigma_t: float, sigma_next: float, *, noise_scale: float = 0.7):
|
||
"""The velocity whose deterministic SDE mean lands exactly on ``sample`` (the max-likelihood
|
||
velocity for a realized FlowGRPO sample). Moving the policy toward it, advantage-weighted, is the
|
||
policy-gradient direction — and is nonzero even at ratio==1, so it is the correct FlowGRPO update
|
||
surrogate (nudging toward the velocity the model already produced would be a no-op).
|
||
"""
|
||
x_t = np.asarray(x_t, dtype=np.float64)
|
||
sample = np.asarray(sample, dtype=np.float64)
|
||
s = min(float(sigma_t), 0.9999)
|
||
dt = float(sigma_next) - float(sigma_t)
|
||
std = float(np.sqrt(s / (1.0 - s)) * noise_scale)
|
||
denom = 2.0 * max(s, 1e-6)
|
||
a = 1.0 + std**2 / denom * dt
|
||
b = (1.0 + std**2 * (1.0 - s) / denom) * dt
|
||
b = b if abs(b) > 1e-8 else (1e-8 if b >= 0.0 else -1e-8)
|
||
return ((sample - x_t * a) / b).astype("float32")
|