精读笔记
Problem Setting
《On the sharp linear convergence rate of the circumcentered–reflection method on subspaces》（arXiv preprint / 2026）研究两个子空间 U,V 上 CRM 的 sharp 收敛率。表面任务是求 P_{U∩V}(x_0)，但论文真正处理的是：CRM 这个非线性 circumcenter 映射在子空间几何下到底以什么 rate 收缩，以及已有 c_F 分析为什么不是 tight。
困难点有两个。第一，CRM 不是线性迭代；不能像 MAP/DRM 那样直接读谱半径。第二，主角度中最坏方向不是单纯 Friedrichs angle 对应方向，而是最小正角 θ_F 与最大角 θ_p 的混合方向。已有分析只用 c_F=cos θ_F，忽略了 θ_p 参与的 conditioning 信息，因此上界偏粗。
关键矛盾是：CRM 是完全 parameter-free 的几何方法，但最优 relaxed MAP rate ρ_V 通常需要知道 θ_F、θ_p 才能设 μ⋆。本文要证明的是，CRM 在初始化到 V 后无需知道这些角，却自动执行了达到同一 worst-case rate 的 adaptive step。
Motivation
已有路线的不足不是“不知道 CRM 收敛”，而是“不知道 CRM 到底快多少”。原始 CRM rate c_F 与 DRM rate 相同，解释力很弱；MAP 的 sharp rate c_F^2 已知，而 relaxed MAP 的最优 rate ρ_V 也已由 Bauschke et al. 给出，但这些都是线性/linesearch 视角。
作者的核心观察是：CRM 虽然定义为三点 circumcenter，但当当前点在 V 内时，circumcenter 被几何约束强制落在 MAP 残差线 v+span(Tv−v) 上。也就是说，CRM 在 V 上隐藏地做了 exact linesearch。这个观察填补了关键缺口：把“几何 circumcenter”与“最优 adaptive relaxation”放到同一个对象里。
动机上最重要的是证明这个 equivalence 不是经验现象，而是由子空间反射的等距性和 principal-angle 结构强制产生的；一旦识别出这一点，sharp rate 就不再来自 CRM 特殊代数，而来自一个标准但此前未用于该问题的 Kantorovich 二阶矩不等式。
Core Idea
论文的核心思想可以压缩为一句：在 V 上，CRM 不是二维 circumcenter 方法，而是沿自伴算子 M=(Id−P_VP_U)|_V 的负残差方向做无参数 exact linesearch。这个重建改变了分析对象：从非线性三点几何映射，变为一个由 M 的两个 moment 控制的一步最优投影。
这种建模方式引入的 inductive bias 是“当前误差只通过 M 的谱分布被处理”。CRM 的 μ_v=<w,Mw>/<w,M^2w> 不需要知道谱端点，但在 worst-case ray 上恰好等于全局最优 relaxation μ⋆=2/(a+b)。因此它既是局部 exact line minimizer，又在全局 worst-case 上匹配最优常数。
与 prior 的本质区别：prior 中 ρ_V 是调参后的 linear relaxation 或 explicit linesearch 的 rate；本文证明纯几何定义的 CRM 在合适初始化后 pointwise 等价于这些 linesearch，并且通过 Kantorovich extremizer 证明 sharpness 和最优性。这不是提出新算法，而是重新解释一个旧算法的真实机制和上限。
Method
1. CRM-on-V 的 linesearch identity：证明 C_T(v)=v+μ_v(Tv−v)。它解决的是 CRM 非线性不可直接谱分析的问题；必要性在于把 circumcenter 映射压缩到一维可控方向；核心变化是 CRM 从 geometric operator 变成 adaptive relaxed MAP。
2. 自伴算子 M 的谱化：定义 M=(Id−P_VP_U)|_V，利用 principal-angle frame 得到 M f_k=sin^2 θ_k f_k。它解决的是 rate 常数从何而来的问题；必要性在于把所有几何量统一为 M 的 moment；核心变化是 θ_F 和 θ_p 成为谱区间端点。
3. Kantorovich inequality：对 <w,Mw>^2/(||w||^2<w,M^2w>) 给下界。它解决 sharp worst-case bound；必要性在于 exact linesearch 的残差正好是这一二阶矩比；核心变化是 rate 不再靠粗 operator norm，而由谱区间 condition number 给出。
4. Sharpness / optimality / non-improvability：构造 worst-case ray sqrt(b)f_{s+1}+sqrt(a)f_p，并证明 over-reflection、single-step parameter-free variants 都无法突破 ρ_V。它解决的是“这个 rate 是否只是分析 artifact”的问题；核心变化是把 ρ_V 从 upper bound 提升为 tight barrier。
5. Memory 与 parametric alternatives：Chebyshev 半迭代可达 ρ_Cheb<ρ_V，但需要 memory 和角度知识；linear two-parameter family 在全空间有 1/3 floor，pre-project 后退化到 relaxed MAP。它们主要用于定位 CRM 的方法谱系，而不是主贡献。
Key Insight / Why It Works
最关键的有效原因是 CRM 在 V 上自动完成了 steepest-line minimization：它选择的 μ_v 正是使 ||v+μ(Tv−v)−v̄|| 最小的 μ。这个 μ_v 由当前点的几何量 dist(v,U) 和 ||Tv−v|| 决定，不需要主角度，但数学上等价于 moment ratio <w,Mw>/<w,M^2w>。因此 CRM 的“无参数性”不是魔法，而是 circumcenter equidistance 条件自动解出了局部最优 relaxation。
真正核心贡献是 Kantorovich 视角。它解释了为什么最坏方向必须混合谱端点，而不是只看 Friedrichs eigenspace。Kantorovich extremizer sqrt(b)u_a+sqrt(a)u_b 同时给出 sharp rate、equality case、over-reflection obstruction 和 single-step lower bound。这个统一性很强，说明论文的增益不是 numerical trick，而是对隐藏谱结构的正确识别。
哪些部分可能只是辅助：大量关于 CDR family、over-reflection、warm-start、实验图的分析更多是在圈定边界。它们有价值，但不是产生 ρ_V 的原因。真正产生 ρ_V 的只有两件事：CRM-on-V 等价 linesearch，以及 M 的 Kantorovich moment bound。
从机制归因看，这不是 scaling、不是 data coverage、不是 retrieval，也不是 benchmark engineering；它是一个 better inductive bias / representation alignment 结果：把迭代限制到 V 后，误差表示与自伴谱结构对齐。若不 warm-start 到 V，这种 alignment 消失，直接 CRM 的 rate 回到 c_F；这也说明 ρ_V 的提升不是 CRM 全局几何的普遍性质，而是特定状态空间选择带来的结构化收益。
Relation To Prior Work
最接近的是三条线：MAP/relaxed MAP，Gubin–Polyak–Raik adaptive projection/linesearch，和 CRM/DRM 反射几何方法。本文不是提出一个全新迭代，而是证明这些线在 V 上发生了 pointwise collapse：A_T(v)=B_T(v)=C_T(v)。
与 Bauschke et al. [8] 的关系尤其关键。ρ_V 不是本文新发现的常数；[8] 已经把它识别为 relaxed MAP family 和 B_T 的最优 rate。本文新增的是：几何 CRM 也无参数达到同一常数，并且给出 sharp equality ray 和 single-step parameter-free optimality。这是实质创新，因为它把 CRM 从“经验上更快的 circumcenter 方法”重新定位为“隐式 optimal adaptive relaxation”。
与原始 CRM 分析 [14,15] 的差异是 rate 机制不同。原分析在全空间给 c_F；本文证明只要初始化到 V，CRM 进入一个完全不同的谱 regime，rate 严格优于 c_F^2（除正交边界）。因此它不是简单改进常数，而是指出旧分析分析的是错误/过大的状态空间。
与 Chebyshev、AAMR/GAP 的关系：这些方法能在某些 regime 更快，但需要 memory 或角度调参。CRM 的位置更像 parameter-free memoryless baseline 的 sharp frontier。它不支配所有 tuned 方法，但明确占据“无参数单步方法的最优点”。
Dataset / Evaluation
这是一篇优化理论论文，evaluation 主要服务于验证解析 rate，而不是证明经验泛化。实验包括 prescribed-angle subspace 构造、随机子空间数值验证、迭代次数比较、warm-start 策略比较、与 AAMR/GAP 的 rate crossover。覆盖范围足以验证核心 claim：给定主角度时，CRM 在 worst-case ray 上确实达到 ρ_V，随机方向不超过该 bound。
但实验并不承担真实应用验证。没有真实信号恢复、图像重建或大型约束系统上的 deployment；这没问题，因为论文的核心 claim 是 sharp theorem。需要注意的是，实验数字本身没有太多信息增量，主要是 sanity check 和可视化。真正 support claim 的是解析证明。
benchmark 没有明显 leakage 或 evaluation bias 问题，因为这里不是学习 benchmark；不过实验不能说明在非子空间、近似投影、噪声反射、多集合 convex feasibility 中同样保持 ρ_V。若读者关心 CRM 在实际非线性集合上的表现，本文 evaluation 支持力度有限。
Limitation
最大限制是 warm-start 到 V。ρ_V 不是 CRM 从任意 x∈R^n 的全局 sharp rate；直接迭代只能保证 c_F，实验也显示通常不会自然接近 ρ_V。这意味着方法的收益依赖状态空间选择：先投影到某个子空间是机制的一部分，不是无关 preprocessing。
第二，问题设定是两个线性子空间，且需要精确正交投影/反射。principal-angle frame、自伴 M、Kantorovich tightness 都严重依赖线性结构。推广到一般凸集、流形、锥或近似投影时，M 不存在或不自伴，核心证明链条会断。文中未充分说明这些 sharp insight 在非线性 CRM 变体中能保留多少。
第三，ρ_V 的“最优性”有严格限定：single-step、沿 T(v)−v、parameter-free。允许 memory 后 Chebyshev 更快；允许角度调参后 AAMR/GAP 在 θ_p≈π/2 且 θ_F 小时更快；允许 full Krylov/CG 理论上还可能更快。所以本文并没有证明 CRM 是整体最优，只是精确刻画了一类方法的边界。
第四，某些 extended family 的讨论增益归因相对分散。linear CDR family、symmetrization、composition、pre-projection 的分析说明了很多负结果，但与主 theorem 的关系更像边界地图。若作为算法论文看，这些部分的 practical payoff 不强；它们主要是理论定位。
Takeaway
- 1. CRM 的关键不是 circumcenter 本身“神奇”，而是在合适子空间上 circumcenter 强制实现了 exact adaptive relaxation；以后分析几何投影法时，应优先寻找这种隐藏的一维 linesearch / self-adjoint reduction。
- 2. Sharp rate 的正确对象不是 Friedrichs angle 单点，而是正谱区间 [sin^2 θ_F, sin^2 θ_p]。
- θ_p 的作用不可忽略；只用 c_F 会系统性低估 CRM-on-V。
- 3. ρ_V 是 parameter-free memoryless 方法的边界；要继续加速必须付出代价：memory、谱/角度估计、或更强的状态建模。
一句话总结
这篇论文把 CRM 在两个子空间上的真实加速机制钉死为“warm-start 后的隐式 adaptive relaxed MAP + Kantorovich-sharp 谱区间收缩”，从而把它定位为 parameter-free、memoryless 投影方法的 sharp frontier，而不是一个仅凭几何直觉更快的 circumcenter heuristic。