精读笔记
Problem Setting
[Riemannian Multilevel Optimization with Application to Constrained Energy Minimization Problems]（arXiv preprint / 2026）
论文处理的是一类离散化后的约束能量最小化：变量落在流形上，例如 Stiefel、ellipsoid、Bernoulli manifold。真正的问题不是写出 Riemannian gradient descent，而是如何让 coarse level 参与当前 fine-level optimization，而不只是作为初始化。
关键矛盾是：multilevel optimization 需要在 coarse 层解一个便宜的 surrogate 来修正当前 fine iterate，但流形上 coarse correction 不能直接相减、不能直接 prolongate 成可行方向，也不能保证对 fine objective 下降。欧氏 MGOPT 的核心对象是 linear correction term 和 restricted gradient；在流形约束下，这两个对象都需要重新定义。
以前路线各有缺口。Cascadic 方法只做 coarse-to-fine continuation，不对当前 fine objective 做一阶一致修正；投影式 Riemannian multigrid 可以在嵌入流形上工作，但通常依赖特定结构；已有 Bernoulli / low-rank Riemannian multilevel 方法更像局部实例化，缺少统一 transfer 语言和收敛理论。本文实际要补的是“Riemannian MGOPT 的一般一阶框架”。
Motivation
作者的核心观察是：MGOPT 成立的关键不是欧氏空间，也不是网格本身，而是 coarse surrogate 在当前点与 fine objective 的一阶信息一致。只要能把 fine gradient 以几何一致的方式限制到 coarse tangent space，并把 coarse displacement 送回 fine tangent space，就可以保留 multilevel correction 的本质。
已有路线缺的是一个抽象层：点如何跨层，tangent vector 如何跨层，gradient correction 如何写成 intrinsic object，以及什么条件下 coarse solve 的下降能推出 fine-level 下降。本文把这些问题统一成 point transfer、vector transfer、lifting map、metric-adjoint Galerkin condition 的组合。
这也是为什么作者选择 Nash coarse model 而不是单纯 nonlinear equation / FAS 表述：optimization formulation 更容易允许 inexact coarse solve、line search 和 descent proof。这个选择很务实，理论上也更适合 Riemannian first-order convergence analysis。
Core Idea
论文真正核心是把欧氏 Nash coarse model 中的 z-y 替换为 coarse manifold 上的 lifting L_y^H(z)，把线性 correction 写成 Riemannian pairing，并让 correction vector w=grad f_H(y)-R_x^y grad f_h(x) 强制 coarse model 在 y 处满足 grad q(y)=R grad f_h(x)。也就是说，coarse objective 可以是独立离散化的能量，但它在当前 coarse point 的一阶信息被 fine objective 校准。
这改变了 coarse model 的建模方式：coarse 层不再只是 lower-resolution objective，也不是 fine objective 的代数 pullback，而是一个被当前 fine first-order residual 修正过的几何 surrogate。它引入的 inductive bias 是“粗层负责低频 / 大尺度 nonlinear structure，fine 层 gradient 负责当前一阶一致性”。这正是 multigrid 在优化里的有效机制，只是本文把它提升到 manifold tangent bundle 层面。
与 prior 的本质区别在于：本文不是为某个流形手写 correction，而是把 transfer operator 的设计空间抽象出来，并指出 Galerkin adjoint 条件是 coherence 和 descent proof 的关键接口。这个抽象比具体实验更重要。
Method
1. Riemannian coarse model：解决的是欧氏 coarse correction 在流形上没有 displacement 和 linear term 的问题。lifting map L_y(z) 把 coarse 点 z 表示为 y 处 tangent displacement；修正项 <w,L_y(z)> 把 fine gradient 信息注入 coarse objective。核心变化是 coarse solve 的输出变成一个 tangent correction，而不是一个需要硬投影的点差。
2. First-order coherence：解决的是 coarse surrogate 可能引入伪 stationary point 的问题。通过 w=grad f_H(y)-R grad f_h(x)，得到 grad q(y)=R grad f_h(x)。这保证如果 fine gradient 的 coarse component 非零，coarse model 在 y 处也看见这个非零一阶信号。
3. Vector transfer operators：解决的是 fine/coarse tangent spaces 不可直接比较的问题。论文给出 restriction-based geometric、restriction-based algebraic、prolongation-based、projection-based consistent/inconsistent 等几类构造。机制上重要的不是这些公式本身，而是 P 与 R 是否满足 metric-adjoint Galerkin condition；满足时 descent proof 和 metric-independence 才比较干净。
4. Descent and line search：解决的是 coarse correction 是否真能降低 fine objective。若 coarse objective 在局部 R-convex，且 coarse solve 产生 q(z)<q(y)，则 prolongated direction 对 fine objective 是下降方向。实际算法再用 fine-level line search 接管 step length，避免 coarse model 过度承诺。
5. Coarse correction condition：解决的是何时 coarse level 仍有用。只有 restricted fine gradient 足够大时才触发 coarse correction；并且不允许连续 coarse correction，强制 fine smoothing 介入。这是理论证明需要的结构，也是 multigrid 直觉中的 smoothing/coarse-correction 交替。
Key Insight / Why It Works
最关键的 insight 是：multilevel acceleration 在这里不是靠 coarse objective 近似 fine objective 的函数值，而是靠 coarse model 对当前 fine objective 的一阶对齐。只要 low-dimensional coarse space 捕获了当前误差中的大尺度成分，coarse solve 就能用更低成本给出一个方向；fine line search 和 smoothing 负责把这个方向变成安全下降。
最可能的核心贡献是 coarse model + metric-compatible transfer 的抽象接口。它把“跨层优化”从欧氏 residual correction 推广到 tangent-level first-order information correction。Proposition 3.2 和 3.3 是方法成立的骨架：前者保证 coherence，后者保证 descent。没有这两点，后面的 multilevel algorithm 只是工程启发式。
metric-independence 的结论也有价值，但需要谨慎理解。它说明固定 P 且 R=P* 时，coarse model 可写成 differential 形式，因而不依赖具体 metric 表达；但实际性能仍强烈依赖 P 是否与 metric / problem geometry 匹配。Bernoulli 实验正好说明：形式上的一致性不够，非均匀 Fisher-Rao metric 下 transfer scaling 会直接决定 line search 是否崩掉。
加速来源主要是 multilevel / scaling，而不是更强的 optimization oracle。它通过 coarse solve 增加 test-time compute，但把 compute 移到低维层；本质是 memory reuse / latent multiscale structure exploitation，而不是改变一阶方法的渐近性质。文中也承认连续 coarse correction 被禁止，因此渐近 rate 很可能继承 single-level smoother；实际收益更多是 transient regime 的误差快速压缩。
辅助部分包括多种 transfer operator catalogue 和三个应用推导。它们证明框架可落地，但不都是同等核心。尤其某些实验里不同 transfer 的结果几乎相同，说明在 Stiefel / ellipsoid 上收益可能主要来自 coarse hierarchy 和已有 preconditioner，而非某个新 transfer 公式。
Relation To Prior Work
这篇论文属于 MGOPT / FAS optimization formulation 与 Riemannian optimization 的交叉谱系。它最接近 Nash MGOPT、Wen-Goldfarb line-search multigrid、Sutti-Vandereycken 的 Riemannian multigrid line search，以及 Bernoulli manifold 上的 multilevel geometric optimization。
与 cascadic multigrid 的本质差异是：cascadic 只是把 coarse solution prolongate 作为 finer initialization，不校准当前 fine gradient；本文的 coarse model 是围绕当前 fine iterate 构造的 correction model，因此它是在线的 current-iterate correction，而不是 continuation。
与 FAS / MGOPT 的关系是：思想并不新，真正新增的是把 linear correction term 几何化。lifting map 替代 z-y，vector restriction 替代 R，metric pairing 替代 Euclidean inner product。这是已有思想的 Riemannian 重写，但不是机械替换，因为 descent proof 需要 R-convexity 和 Galerkin adjoint 条件。
与 [44] 和 [33] 相比，实质创新在统一性和理论闭合。[44] 更偏 embedded / projection-based low-rank setting；[33] 更偏 Bernoulli 与 prolongation-based transfer。本文把 intrinsic/extrinsic、restriction/prolongation/projection 几类 transfer 放在同一个一阶 coherence 框架里，并给出 global stationarity proof。这是工程上可迁移的抽象，而不只是多一个应用。
Dataset / Evaluation
evaluation 覆盖面是本文较强的部分：Kohn-Sham、Gross-Pitaevskii、binary continuous cuts 三个问题的几何、离散化和能量结构差异很大。它确实支持“框架不是只适用于单一流形或单一网格”的 claim。尤其 Bernoulli / Fisher-Rao case 暴露出 metric-compatible transfer 的必要性，比单纯报 speedup 更有信息量。
但实验没有完全回答增益归因。Kohn-Sham 中 baseline 包含 H1RGD / H1RCG，multilevel 方法也使用 H1 preconditioning 和 nested coarse solvers；Gross-Pitaevskii 使用 energy-adaptive metric；continuous cuts 的结果强依赖 transfer scaling 和 Armijo 行为。因此 speedup 不能简单归因于 Riemannian coarse model 本身，而是 coarse hierarchy、preconditioner、line search、coarse solve schedule 的组合收益。
benchmark 主要验证 computational efficiency 和 broad applicability，没有验证更强的泛化意义。这里的“generalizable”是数学框架可实例化到不同 manifold/discretization，而不是算法无需问题知识即可自动有效。每个应用仍需要手工设计 point maps、transfer operators、metric choices 和 coarse objectives。
没有真实部署层面的实验问题，因为这是数值优化论文；但从 numerical evidence 看，claim 的边界应表述为：在存在自然多尺度结构且 transfer 设计合理时，能显著降低 time-to-stationarity。
Limitation
第一，descent guarantee 依赖 coarse objective 的局部 R-convexity 或 fallback 检查。对 Kohn-Sham、Gross-Pitaevskii 这类非凸能量，文中更多是通过局部/实践机制维持下降；全局非凸结构下 coarse correction 何时可靠，文中未充分说明。
第二，方法把难点从“如何优化 fine problem”部分转移到“如何设计跨层几何 transfer”。point restriction/prolongation 并非自动存在，也不一定数值稳定。Bernoulli 例子说明，错误 scaling 的 transfer 会让 line search 退化，甚至使 coarse correction 频繁但无效。
第三，global convergence 只到一阶 stationarity，且证明依赖 Wolfe、compact sublevel set、Lipschitz pullback gradient、bounded gradient step size，以及禁止连续 coarse correction。这个理论并不解释为什么更快，也不给 rate。加速机制仍主要是经验性的 multilevel scaling。
第四，level 数、coarse solve 精度、coarse condition 参数 η/μ 都是关键 engineering knobs。文中展示 adaptive cycle 比 fixed cycle 好，但没有系统理论说明如何选。不同应用的 sweet spot 不同，说明方法上限受 hierarchy cost model 强烈约束。
第五，metric-independence 容易被过度解读。coarse model 表达可以 metric-independent，但实际 direction quality 不独立于 metric；特别是在强非均匀 metric 上，transfer 是否 metric-conjugated 决定性能。这一点是优点，也是限制。
Takeaway
- 1. 最值得记住的是：Riemannian multilevel optimization 的核心接口应放在 tangent-level first-order coherence，而不是 point-level interpolation。
- 跨层点映射只是生成 vector transfer 的手段。
- 2. Nash coarse model 的 Riemannian 版本是一个可迁移模板：用 lifting 写 displacement，用 metric-adjoint transfer 写 gradient correction，用 line search 保证 fine-level safety。
- 这个模板可用于其他 constrained variational problems。
一句话总结
这篇论文把 MGOPT 的一阶一致 coarse correction 几何化为 Riemannian coarse model 与 metric-compatible tangent transfer，是 Riemannian optimization 中 multilevel scaling 机制的一次统一化和理论化，而不是一个全新的单层优化算法。