Low-rank one-shot horizon predictor
Implementation & benchmark of arXiv:2607.13289 — Model predictive control for laser thermal processing: operator learning, closed-loop validation, and out-of-distribution analysis
Source paper: Model predictive control for laser thermal processing: operator learning, closed-loop validation, and out-of-distribution analysis arXiv:2607.13289 ⓘ · analyzed Aug 30, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Replace an autoregressive rollout of a learned dynamical model with a branch-trunk factorization that predicts all future steps simultaneously. The branch network encodes the future action sequence, while the trunk network encodes the current state and query coordinates; their inner products produce the complete horizon. This removes repeated state updates during inference and gives a compact differentiable model for planning.
Formulas
Mathematical statement
The paper uses the finite-rank separable representation \(\widehat{T}_{\max,k+i}=\sum_{j=1}^{p}B_{ij}(\mathbf{u}_{k})T_j(\mathbf{y}_{k})+b_i\). Here \(N_h\) is the prediction horizon, \(\mathbf{u}_k\) is the encoded future input sequence, \(\mathbf{y}_k\) is the encoded current state and query information, \(B(\mathbf{u}_k)\in\mathbb{R}^{N_h\times p}\) is the branch output, \(T(\mathbf{y}_k)\in\mathbb{R}^{p}\) is the trunk output, \(p\) is the latent basis rank, and \(\mathbf b\) is a trainable horizon bias. This is a low-rank separable approximation of a nonlinear map from an input sequence and state/query pair to an output trajectory. The transferable property is that one branch evaluation and one trunk evaluation produce every horizon output instead of recursively evaluating a one-step model \(N_h\) times. For a vector-valued forecast with output dimension \(d\), use \(B(\mathbf u)\in\mathbb R^{N_h\times p\times d}\) and \(\widehat{Y}_{i,:}=\sum_{j=1}^{p}B_{i,j,:}(\mathbf u)T_j(\mathbf y)+b_{i,:}\).
Implementation notes
Integrate this at the world-model or differentiable-planning interface, replacing an autoregressive model call that currently predicts the next state, feeds it back, and repeats for \(N_h\) steps. Construct training examples from contiguous trajectories: the branch input is the full future action sequence \(\mathbf u_t=(a_t,\ldots,a_{t+N_h-1})\) plus optional exogenous inputs; the trunk input is the current state \(\mathbf y_t=s_t\) and normalized query coordinates \(q_i=i/N_h\). Use two MLPs with latent width \(p\). The branch output should have shape \([B,N_h,p,d]\), and the trunk output shape \([B,p]\). Compute predictions with Yhat = einsum('bjpd,bp->bjd', Bcoef, T) + bias. Pseudocode is Bcoef=Branch(u); T=Trunk(concat(s,q)); Yhat=einsum('bjpd,bp->bjd',Bcoef,T)+bias; loss=mse(Yhat,Y)+lambda*mean(Bcoef**2). The mathematical component is the finite-rank separable factorization; rank, network widths, and regularization are empirical choices. Train on teacher-generated trajectory windows with normalized state and action channels. The first experiment should use Pendulum or CartPole with horizons 8, 16, and 32, comparing against an equally sized one-step MLP rolled out autoregressively. Measure wall-clock latency, GPU memory, horizon-wise MSE, gradient-based action-planning reward, and error under shifted action distributions. Success means similar or better long-horizon error with lower latency and fewer sequential kernel launches; report loss and planning performance at equal FLOPs.
Verification
Stage 1 · Toy mechanism gate: Passed ✓
Stage 2 · Mechanism transferred to benchmark: Confirmed ✓
Stage 2 · Practical benchmark result: Significant win ✓
Methodology: 8 paired seeds, permutation test (p<0.05) against a learning-rate-tuned baseline. How verification works
Stage 1 — Mechanism check agent confidence 9/10
Built a reproducible low-rank branch/trunk one-shot horizon predictor and compared it with an autoregressive one-step MLP. The finite-rank check matched the omitted-singular-value prediction to numerical precision, reaching approximately zero error at rank 7 for a rank-7 operator. Measured inference latency showed the predicted horizon scaling, while the H=16 one-shot model also achieved lower test MSE and lower latency in the synthetic benchmark.
- Agent confidence
- 9/10
- Baseline
- Rank sweep matched the singular-value-tail prediction; H=16 rollout mean horizon MSE 0.002229, final-step MSE 0.002600, and latency 9.376 ms per batch of 128.
- Idea
- Rank-7 observed relative MSE 2.22e-30 versus 2.38e-30 predicted; H=16 one-shot mean horizon MSE 0.000528, final-step MSE 0.000601, and latency 0.744 ms per batch of 128, giving 12.61x lower measured latency. Baseline/idea timing ratios for H=4,8,16,32,64 were 2.60, 5.56, 10.27, 14.50, and 40.62.
Limitations: The mechanism sweep uses a synthetic finite matrix and exact SVD rather than a trained nonlinear operator. The learned benchmark uses a synthetic scalar dynamical system, one horizon, and non-equal parameter/FLOP budgets; it does not test Pendulum or CartPole, action-gradient planning reward, shifted action distributions, memory, or a real world-model interface.
How to run: python3 experiment.py > results.json
Stage 2 — Benchmark vs. tuned baseline (arXiv:2607.13289)
Beats the tuned baseline by 86.3% (p=0.0081, wins 8 / 8 paired seeds; metric: lower is better).
- Benchmark
- Controlled_pendulum_horizon
- Model
- local_shared_horizon_system
- Paired seeds
- 8
- Baseline mean
- 0.0484
- Idea mean
- 0.0067
- Effect (Δ)
- -0.0418 (−86.3%; negative = idea better)
- Wins
- 8 / 8 paired seeds
- p-value
- 0.0081 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- n/a
- Mechanism
- Confirmed ✓
- Practical effect
- Helps
- Baseline tuning
- swept over 3 configs
Only the custom 8-step controlled pendulum horizon was tested. Parameter/FLOP parity, horizons beyond 8, action-gradient planning reward, shifted action distributions, and real-world or standard Pendulum environments were not tested.
How to run: python3 bench_run.py
Verdict computed by deterministic test code from paired-seed statistics — not by the language model.
Stage 2 — Benchmark latest Worked ✓
Agent confidence: 8/10
Implemented a structurally matched custom controlled-pendulum horizon track and compared an autoregressive one-step MLP rollout with an end-to-end low-rank branch/trunk predictor. Across 8 paired seeds, the idea significantly improved test MSE and reduced measured inference latency; the mechanism signature confirmed one-shot versus sequential behavior. The idea used more parameters, so the result is not an equal-parameter win.
- Baseline
- Best lr=0.006; mean test MSE 0.0484111, std 0.0098726; 2642 parameters; 3.3496 ms/batch32.
- Delta mean
- -0.0417583
- Idea
- Best lr=0.006; mean test MSE 0.00665281, std 0.0011661; 14876 parameters; 0.7940 ms/batch32; 4.218x speedup.
- Idea wins
- 8
- N pairs
- 8
- Permutation p value
- 0.0081
Benchmark result
Beats the tuned baseline by 86.3% (p=0.0081, wins 8 / 8 paired seeds; metric: lower is better).
- Benchmark
- Controlled_pendulum_horizon
- Model
- local_shared_horizon_system
- Paired seeds
- 8
- Baseline mean
- 0.0484
- Idea mean
- 0.0067
- Effect (Δ)
- -0.0418 (−86.3%; negative = idea better)
- Wins
- 8 / 8 paired seeds
- p-value
- 0.0081 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- n/a
- Mechanism
- Confirmed ✓
- Practical effect
- Helps
- Baseline tuning
- swept over 3 configs
Only the custom 8-step controlled pendulum horizon was tested. Parameter/FLOP parity, horizons beyond 8, action-gradient planning reward, shifted action distributions, and real-world or standard Pendulum environments were not tested.
How to run: python3 bench_run.py
Artifacts
- bench_report.json 3.7 KB View Raw
- bench_run.py 4.7 KB View Raw
- custom_horizon_dynamics.py 1.2 KB View Raw
- experiment.py 6.1 KB View Raw
- patch.py 1.2 KB View Raw
- report.md 1.6 KB View
- report_bench_2026-08-30T175229.md 3.6 KB View
- results.json 3.0 KB View Raw
- results2.json 3.8 KB View Raw