精读笔记

Problem Setting

论文处理的不是一个单一物理任务，而是一类共同结构：高维变量、线性或近线性物理前向模型、大量上下界/剂量/投影约束，以及一个并非必须严格最小化的软目标。典型例子是图像重建中的去噪/TV 偏好，或放疗计划中的 body dose 降低。

真正困难点在于这些问题经常处在 constrained optimization 与 feasibility seeking 的灰区：约束可能不一致，目标函数常常只是 preference 而非 hard requirement，严格最优解既可能不存在，也未必是实际需要的对象。传统优化把所有东西合成一个目标或约束系统，遇到 infeasibility 或调权问题会非常脆弱；单纯 feasibility seeking 又忽略了软质量指标，容易得到可行但不够好的解。

因此关键矛盾是：如何在不牺牲可行性/近可行性收敛的前提下，把软目标稳定地注入迭代过程。本文的回答是 superiorization：不是解优化问题，而是改变 feasibility algorithm 的轨迹，使其落到 objective value 更低的可行/近可行点。

Motivation

已有 superiorization 理论和应用并不少，尤其在 CT、radiotherapy、投影算法社区已经有长期积累。缺口不在概念，而在可用性和系统比较：很多工作是单任务、单实现、单参数配置，难以复现，也难以回答“换一个 feasibility operator 或 perturbation strategy 会怎样”。

作者的核心观察是，superiorization 很适合那些目标函数只是软偏好的物理问题，但其传播受限于工程基础设施。没有开放、模块化、支持 GPU 的工具，研究者很难把这种方法当作通用实验范式，而只能把它看成某个应用里的 heuristic。

所以本文的动机有两层：算法层面，展示 superiorization 在不一致约束和软正则化目标下的实用性；系统层面，把 feasibility operators、perturbations、CPU/GPU execution 解耦，降低新物理问题迁移成本。后者其实是本文最明确的贡献。

Core Idea

核心思想是把 constrained optimization 的单一目标化建模拆开：硬约束交给 projection-based feasibility seeking，软目标只通过 bounded perturbations 影响迭代路径。这样目标函数不需要被精确最小化，也不需要与约束权重合并成一个难调的 penalty objective；它只负责在每轮 feasibility update 之前轻推当前点向更低 objective value 的方向移动。

这个设计的理论直觉来自 bounded perturbation resilience：只要扰动方向有界、步长总和有限，某些 projection algorithms 的可行性收敛性质不会被破坏。于是 perturbation 变成一种受控 inductive bias，而不是主优化器。它改变的不是 feasible set 本身，而是 feasibility dynamics 在可行集/近可行区域中的吸引位置。

与常规 optimization 的本质差异是目标函数地位被降级：它不是要定义全局最优，而是选择“更好”的可行解。与普通正则化的差异是，正则项没有被并入主目标求解，而是作为 test-time trajectory shaping 插入迭代。这个视角对不一致约束尤其自然，因为算法可以追求 closest proximity，同时用软目标改善实际质量。

Method

方法关键机制可以压缩为四个部分。

第一，base algorithm 必须是对有界扰动稳定的 feasibility-seeking operator。它解决的是“扰动后是否还能收敛/接近约束”的问题。没有这一点，superiorization 就只是普通 gradient heuristic，无法区别于随便在投影之间插梯度步。

第二，perturbation 使用归一化梯度或次梯度，并配合递减 power-law 或 adaptive step-size。归一化的作用是控制方向序列有界，步长规则的作用是控制扰动总量。它带来的核心变化是：目标函数下降被限制在不会主导 feasibility dynamics 的范围内。

第三，restart perturbation 是本文实验中很重要但理论上微妙的机制。递减步长后期会让 perturbation 几乎失效，restart 试图恢复扰动强度，同时通过让指数参数整体增长来维持 boundedness。CT 实验显示 restart 对抑制 error curve 的 L-shape 和减少噪声过拟合有明显作用；这部分很可能是实际效果的重要来源之一。

第四，SupPy 的架构把 feasibility operator 与 perturbation strategy 解耦，并支持 NumPy/SciPy 与 CuPy 后端。这解决的是研究效率与 scaling 问题，不是新的数学机制。GPU 使大矩阵 simultaneous methods 可在实用时间内比较，尤其对 CT 和放疗这类矩阵操作密集任务很关键。

Key Insight / Why It Works

这篇论文最值得抓住的 insight 是：在许多物理问题里，“求最优”是过度建模，“求一个被软目标偏置的近可行点”反而更符合实际需求。superiorization 的有效性并不来自比优化器更强的全局搜索，而来自它把目标函数影响限制在一个不会破坏约束收敛的低能量通道里。

真正起作用的部分可能有三类。第一是 inductive bias：TV/L1/mean dose 这些扰动方向本身就是强先验，能把 feasibility solution 从噪声拟合或剂量扩散的方向拉回来。第二是 test-time compute：算法通过反复投影和扰动在推理阶段调整解，而不是依赖离线训练或固定闭式解。第三是 algorithm selection：EMR/CG 本身在 CT 上明显比 ART/EL 更合适，superiorization 的收益建立在 base operator 已经能给出好轨迹的前提上。

需要直说的是，论文里的很多增益不应归因于“superiorization 这个概念本身”单独带来的突破。CT 中 TV-superiorization 的改善很大程度上就是把已知 TV 先验以扰动形式注入；restart 和 early stopping 也承担了抑制噪声过拟合的作用。放疗中 body dose 降低则基本符合 mean-dose perturbation 的直接预期，并不说明它解决了更深的 plan optimality 问题。

GPU 加速主要是 engineering / scaling 贡献。它使 simultaneous projection methods 的比较变得现实，但没有改变算法的统计性质或数学保证。SupPy 的价值在于把已有 superiorization 机制做成可组合系统，而不是提出新的 superiorization 原理。

文中未充分说明的是，在 infeasible case 下，projection algorithm 收敛到的 closest proximity point 是否与任务质量有稳定对应。尤其 radiotherapy 中，临床可用性来自剂量指标和视觉 DVH 比较，但这不等于证明 superiorization 在复杂约束冲突下有系统可靠性。

Relation To Prior Work

本文属于 projection methods、superiorization、inverse problems / radiotherapy planning 交叉谱系。最接近的是 Censor 等人的 bounded perturbation resilience / superiorization 理论、Penfold/Sidky 系列 TV-superiorization CT 重建、Nikazad 的 adaptive step-size，以及 radiotherapy 中 split feasibility + superiorization 的既有工作。

真正不同点不是算法公式，而是系统化组织：把 feasibility algorithms、perturbation strategies、split feasibility、CPU/GPU 后端做成同一工具箱，并用三个物理场景展示可迁移性。换句话说，它新增的信息主要是“平台化后 superiorization 是否能作为通用实验范式工作”，而不是“一个新的 superiorization 方法是否理论上更强”。

看似新的部分，例如 TV/L1/L2 perturbation、power-law step-size、restart、adaptive perturbation、CQ/split feasibility、EMR/CG/DROP，本质上都是已有思想的重组和实现集成。实质创新在工程抽象和跨应用验证：研究者可以在相同框架下比较 base operator 与 perturbation，而不是为每个应用重写算法。

与 proximal splitting、ADMM、projected gradient、penalty/barrier optimization 相比，superiorization 的差别在于不追求某个明确 composite objective 的最优性。它牺牲 optimality claim，换取在 inconsistent / soft objective 场景下的鲁棒工作流。这是定位上的差异，不是一定意义上的性能优势。

Dataset / Evaluation

评估覆盖 seismic tomography、LoDoPaB-CT 和 CORT radiotherapy，任务跨度较大，能支持“SupPy 可迁移到多类物理约束问题”的 claim。尤其 CT 和放疗都不是玩具问题，矩阵规模和实际物理语境足够说明 GPU 后端和模块化实现有价值。

但这些实验更像 capability demonstration，而不是严格 superiority benchmark。seismic 中提升较小，L2 甚至无收益；CT 中有效性主要集中在 EMR/CG + TV 组合，ART/EL 表现较差；放疗中展示了 infeasible constraints 下可产生合理计划，但缺少与成熟优化规划器在同一约束/目标设定下的公平系统对比。

评估确实支持两个结论：一是 superiorization 通常能在 feasibility seeking 基础上改善某些软指标；二是 GPU SupPy 能显著降低 simultaneous methods 的运行成本。它不充分支持更强结论，例如 superiorization 普遍优于 constrained optimization、或在临床复杂 planning 中更可靠。

另外，早停准则、restart、扰动强度对结果影响很大，文中虽然做了部分比较，但没有完全拆解增益来源。benchmark 没有明显 leakage 问题，因为这不是 learning-based 方法；但存在 evaluation bias：选择的软目标与报告指标高度一致，因此改善在一定程度上是预期内的。

Limitation

第一，方法成立依赖 base feasibility algorithm 的性质。若底层 operator 对不一致约束不稳定，或对有界扰动不具 resilience，superiorization 没有独立保证。CT 中 ART/EL 的失败已经说明：扰动不能拯救错误的底层动力系统。

第二，增益归因不清。TV/L1/mean-dose 本身就是强先验，restart 和 early stopping 也在实际效果中扮演重要角色。论文没有充分区分“superiorization 框架带来的收益”和“经典正则化/调参/早停带来的收益”。

第三，optimality 上限是内生的。superiorization 明确不追求最优，因此当任务真的需要最优权衡、Pareto frontier 或严格临床 objective trade-off 时，它只能给 heuristic solution。它适合“够好且稳健”的场景，不适合需要可证明最优质量的场景。

第四，scalability 仍受矩阵存储和扰动计算限制。GPU 加速解决了一部分密集矩阵操作问题，但 CT 和 radiotherapy 的真实大规模部署可能需要 matrix-free operator、专用 CUDA kernel 和更细的内存管理。作者也承认这仍是未来工作。

第五，infeasible constraints 下的理论-实践连接不足。closest proximity 的数学定义未必等价于物理或临床可接受性；特别是 DVC 约束收敛慢，且非线性/排序式投影可能使实际行为更依赖实现细节。文中未充分说明这些情况下的稳定性边界。

Takeaway

• 1. 这篇论文真正推动的是 superiorization 的可实验化和平台化，而不是提出新的 superiorization 理论。
• 2. 对物理反问题而言，把软正则化目标作为 bounded trajectory bias，而不是硬塞进 constrained optimization objective，是一个值得迁移的建模视角；尤其适合 infeasible 或 soft-preference-heavy 的场景。
• 3. 未来真正值得做的是自动选择 perturbation / restart / stopping rule，并把增益归因从 heuristic 调参中拆出来；否则 superiorization 仍会停留在“有效但靠经验”的方法族。
• 4. SupPy 的最大价值可能是成为比较 projection algorithms 与 perturbation strategies 的实验基座；如果后续加入 matrix-free、大规模 CUDA kernel 和非次梯度扰动，它的研究价值会明显上升。

一句话总结

这篇论文把已有 superiorization 思想工程化为 GPU 友好的通用物理约束求解平台，真正贡献在于将“可行性收敛 + 软目标轨迹偏置”变成可复用实验范式，而不是在优化理论上给出新的突破。