精读笔记
Problem Setting
这篇论文实际处理的是非光滑凸变分问题的有限元离散求解：目标泛函形如 ∫φ(x,∇u)+∫ψ(x,u)，覆盖 TV、p-Dirichlet、障碍/梯度约束、弹塑性扭转等。关键矛盾是：这些问题的自然最优性条件是次微分包含关系，直接 Newton 化要么不可微，要么只在特定问题中可写成 active-set 结构；而连续函数空间中即使 pointwise prox 是 Lipschitz/半光滑，也会遇到 norm-gap，难以得到可用的 Newton derivative。以前方法往往在 primal 或 dual 单侧工作，局部很快但对网格、退化区、约束激活集和初始化敏感。本文真正要解决的是：能否构造一个统一的 primal-dual residual，使非光滑性以 prox 形式出现，并在离散层面给出良定性和局部超线性收敛。
Motivation
已有路线缺的是一个结构统一的 Newton 方程，而不是又一个针对 TV 或 p-Laplace 的 solver。Fenchel duality 已经告诉我们 primal 变量 u 和 flux/dual 变量 z 的最优关系，但传统做法通常把它转成 Euler-Lagrange 方程、active-set 条件或问题专用的互补系统。作者的观察是：次微分包含 z ∈ ∂φ(∇u) 与 div z ∈ ∂ψ(u) 可以通过 resolvent/prox 完全等价地改写为固定点方程；prox 是 1-Lipschitz 且 pointwise Newton differentiable，这给半光滑 Newton 一个统一入口。关键缺口在于如何让这个入口在有限元离散后仍保留 Fenchel 对偶结构、局部解耦 prox 结构和线性化可逆性。
Core Idea
核心思想不是“用 prox 做平滑”，而是用 prox 重新组织 primal-dual optimality system。原来的非光滑性位于 ∂φ、∂ψ；本文把它移到 prox_{γφ}、prox_{γψ} 中，然后对 residual F_h(z_h,u_h) = [∇u_h - prox_{γ1 G_h}(∇u_h+γ1 z_h), u_h - prox_{γ2 F_h}(u_h+γ2 div_h z_h)] 做半光滑 Newton。这样 Newton 看到的是一个 Lipschitz 固定点残差，而不是裸露的多值次微分。
本质区别在于它把 flux 作为独立未知量引入 Newton 迭代，而不是从 primal 解后处理得到。这个 primal-dual 对称性很重要：dual residual 不是附属诊断，而是参与 conditioning 和线性化。方法的 inductive bias 是 Fenchel pairing：∇u 与 z、u 与 div z 成对出现，prox 参数控制二者的耦合尺度。相比 primal-only Newton，它更像在求解一组“局部 constitutive law + 全局 equilibrium”的混合系统，因此在退化、饱和、约束激活区域更稳定。
Method
方法的关键机制可以压缩为四点。第一，连续 Fenchel 框架给出 primal/dual 能量和强对偶关系，保证最优性包含 z ∈ ∂φ(∇u)、div z ∈ ∂ψ(u) 不是人为引入的 mixed formulation，而是原问题的精确对偶结构。第二，离散化选用 V_h=P1 连续函数与 Y_h=P0 分片常向量场，并定义 div_h 为梯度的负伴随；这让 ∇V_h ⊂ Y_h，prox_G_h 在单元上解耦，prox_F_h 在节点或 L2 结构下解耦。第三，proximal residual F_h 将多值包含转成单值方程，Newton derivative 由 prox 的 Newton derivative 组成，谱界来自 prox 的 firm nonexpansiveness、强凸性或 Moreau envelope。第四，Newton 线性系统虽然是 primal-dual saddle-point 形式，但可通过 Schur complement 化成 primal 方向的椭圆问题，再回代 dual 方向；这解释了为什么实现上仍接近标准 FEM elliptic solve，而不是昂贵的全 saddle solve。
Key Insight / Why It Works
真正有效的原因是 prox residual 同时具备两个性质：它保留了原始非光滑问题的精确最优性，又把局部非光滑 constitutive relation 变成了单值、Lipschitz、局部可 Newton 化的映射。对 TV 和梯度约束这类问题，难点通常集中在 |∇u|=0、|z|=1 或饱和边界附近；prox 的 Newton derivative 在这些区域给出有界的 generalized slope，而不是让 Hessian 退化到不可控。
最核心贡献是 primal-dual invariance 与离散 Fenchel-compatible Newton residual。它不是简单的 scaling，也不是靠更多 test-time compute；它是更合适的变量组织和 residual 选择。引入 z_h 后，方法不再要求 primal 变量单独承载 flux 信息，这在 p 很接近 1、p 很大、TV 正则化或梯度约束下尤其关键。文中实验里 prox 方法与传统 Newton 的局部下降都很快，说明优势不在“Newton 阶数更高”，而在进入局部收敛域更稳定、线性化系统更少被网格/退化结构破坏。
也要直接说：部分收益可能来自工程层面的 residual 与 globalization 设计，包括 gradient-flow initialization、Armijo line search、Huber/Moreau 正则化和参数 ε∼h 或 h² 的选择。论文理论证明了局部超线性和若干全局良定性，但鲁棒性增益的细粒度归因并不完全清楚；尤其不同 γ、ε、mesh regime 下 conditioning 的系统分析还不够。
Relation To Prior Work
它最接近三条线：Chambolle-Pock/first-order primal-dual splitting 的 proximal optimality view，Hintermüller-Ito-Kunisch/active-set 型半光滑 Newton，以及 proximal Newton / Moreau-Yosida regularization。看似新的 prox 方程本身并不新，Moreau decomposition 和 resolvent 重写都是标准工具；新意在于把这些工具嵌入一般凸变分积分泛函的离散 Fenchel duality，并把它变成一个统一的 semi-smooth Newton residual。
和 primal-dual active set 的本质差异是：本文不是先识别 active set 再写互补条件，而是通过 prox derivative 自动编码 active/saturated/linear 区域。和 Chambolle-Pock 的差异是：Chambolle-Pock 把 prox fixed point 用于一阶迭代，本文把同一个 fixed point 作为 Newton 方程来解。和传统 primal Newton 的差异是：传统方法对 DI_h 或 D²I_h 动手，退化结构直接进入 Hessian；本文对 primal-dual prox residual 动手，把 constitutive nonlinearity 局部化并保留 flux 变量。实质创新在“统一离散 Fenchel-prox-Newton 闭环”，不是单个 prox 公式。
Dataset / Evaluation
评估不是数据集意义上的 ML benchmark，而是三个代表性 PDE/变分问题族：TV-minimization、p-Dirichlet、elasto-plastic torsion。覆盖面选择得比较有针对性：线性增长、p-增长极端情形、∞-增长/梯度约束，基本对应这类方法最容易出问题的三种结构。实验对比的是问题专用 primal 或 dual Newton，并在网格细化、mass lumping / non-mass lumping、两种 globalization 策略下观察 residual decay。
实验支持的 claim 是“prox-primal-dual Newton 在离散 refinement 下更鲁棒”，而不是“局部收敛显著更快”。文中没有大规模真实工程 deployment，也没有自适应网格、复杂几何、多物理耦合或非凸能量测试。benchmark 选择与理论假设高度一致，因此验证的是框架在设计域内有效；跨出 convex variational / Fenchel-clean setting 后是否仍有优势，文中未充分说明。
Limitation
最大限制是连续层面的问题没有真正解决。作者明确指出 continuous F 的 Newton differentiability 受 norm-gap 阻碍，因此主要理论落在有限维离散问题上；这对数值方法足够实用，但不能解读为无限维半光滑 Newton 理论的突破。第二，良定性证明依赖 prox derivative 的正定性、上界小于 identity、或 Moreau/强凸结构；TV 情况还需要额外论证，说明通用定理并不能直接覆盖所有重要边界情形。第三，方法把一部分困难转移到 prox evaluation、Newton derivative selection 和参数耦合上：p-Dirichlet 的 prox 需要逐单元解标量方程，γ 和 ε 的选择会影响 conditioning。第四，实验增益的归因不完全干净：gradient-flow initialization、line search、dual variable 独立更新、正则化尺度都可能贡献鲁棒性，文中没有系统 ablation。第五，所有问题都是凸的，且离散 Fenchel duality 很干净；对非凸材料模型、不可分 integrand、强耦合系统，prox 的局部解耦优势可能消失。
Takeaway
- 最值得记住的不是某个具体 TV 或 p-Laplace 公式，而是：把非光滑变分问题的 Newton residual 设计在 Fenchel primal-dual prox optimality 上，比直接 Newton 化 primal Euler equation 更稳。
- 这个思路可以迁移到其他 convex PDE-constrained 或 variational inequality 问题：先找对偶 pairing，再用 resolvent 把次微分封装成 prox，再在离散层面构造 Newton 方程。
- 未来真正值得做的是参数/conditioning 理论、自适应网格下的可靠 stopping 与 line search、以及突破 norm-gap 的无限维分析。
- 它推动的是“proximal fixed-point view 从一阶算法工具变成二阶 Newton residual 设计原则”。
一句话总结
这篇论文把 Fenchel 对偶、proximal optimality 和半光滑 Newton 组合成一个统一的离散 primal-dual 求解范式，其贡献主要是 residual 设计和良定性机制，而不是新的 prox 算子或更快的局部 Newton 阶数。