精读笔记
Problem Setting
这篇论文实际解决的是finite-cardinality set-valued objective map的无约束set optimization，在lower set preorder下寻找stationary point。目标F(x)不是单个向量，而是一组向量值函数的像集；排序由锥K诱导，且K只要求闭、凸、尖、solid，不要求polyhedral或有限生成。
真正困难点在于：set-valued目标的“下降”不是单一梯度方向问题。每个x处只有minimal elements真正参与下集序比较，而minimal elements本身会随x变化；同一个minimal value还可能由多个f^j产生，导致active index组合P_x。算法既要处理这种活动集合切换，又要给出类似CG的方向递推和Wolfe线搜索。
以前方法主要卡在两个地方：一是有限生成锥假设使scalarization和Wolfe分析比较容易，但排除了重要非polyhedral cone；二是部分set optimization一阶方法依赖最优点正则性，理论适用范围窄。这里的关键矛盾是：想保留CG方法的低内存和加速特性，同时又不能破坏锥序意义下的充分下降。
Motivation
已有路线不够的地方不是缺一个新的beta公式，而是缺一个能在set preorder、非有限生成锥、无最优点正则性三者同时存在时闭合的CG框架。Steepest descent在理论上直接，但每步只用当前最陡方向，数值上通常保守；FR/CD/DY/PRP/HS等set-valued CG已有推进，但HZ这一类带保证下降的参数尚未被完整搬到该setting。
作者的核心观察是：HZ在scalar/vector optimization中的优势来自把共轭参数设计成“既利用前一方向，又显式压制破坏下降的曲率项”。如果能把vector inner product中的梯度差、方向导数等对象替换为Drummond-Svaiter scalarized max descent quantity，就可能在set-valued目标下复用HZ的保证下降结构。
关键缺口是理论层面的：如何在不枚举有限生成锥极射线的情况下定义Wolfe条件、证明步长存在，并让HZ方向满足K-descent。
Core Idea
论文真正核心的方法思想是：不把set-valued objective先压成一个全局固定的scalar objective，而是在每个迭代点只看当前minimal image set，并对这些minimal elements的候选active indices构造局部下降模型。V_x(a,d)=max_j phi(grad f^{a_j}(x)^T d)+1/2||d||^2给出一个“锥序意义下的最陡下降代理”；u_k是这个代理的唯一最小方向，||u_k||同时作为stationarity residual。
在这个局部代理之上，作者引入HZ型非线性CG方向d_k=u_k+beta_k d_{k-1}。本质变化是：CG不再依赖单个目标梯度，而依赖当前active minimal family的scalarized directional derivatives。HZ参数中的负修正项扮演了防止共轭项破坏下降的角色；beta截断为非负，使方向在不利情况下退化回steepest-like direction。
与prior的本质区别在于，论文把“非有限生成锥”通过C={w in K*: w^T e=1}和DS scalarization吸收进phi，而不是依赖有限极射线枚举；同时把HZ下降保证嵌入set-valued minimal-index结构中。它的generalization主要来自锥表示方式更一般，而不是来自更强的搜索策略。
Method
1. 当前minimal set与partition：它解决的是set-valued目标中哪些分量真正参与lower-order比较的问题。只看Min(F(x),K)避免把被支配元素误纳入下降条件；partition P_x处理多个函数产生同一minimal image value的歧义。核心变化是把全局p个函数缩成当前有效的minimal representatives。
2. Drummond-Svaiter scalarization phi：它解决的是一般锥序下“向负内点方向下降”如何被标量判定的问题。phi(y)<0等价于y属于-int(K)，因此max_j phi(grad f^{a_j}(x)^T d)<0就是K-descent的可计算判据。它是整篇理论能脱离有限生成锥的关键接口。
3. V_x子问题产生u_k：它解决的是没有单一梯度时如何定义stationarity residual和基础下降方向的问题。max-linear scalarized directional derivative加上1/2||d||^2形成强凸模型，保证u_k唯一；非零u_k对应非stationarity，||u_k||趋零对应渐近stationarity。
4. Wolfe条件的set-valued版本：它解决的是步长既要让集合值在lower preorder下下降，又要控制方向导数变化的问题。Armijo-like条件写成F(x+alpha d)相对当前minimal values的下集改善，curvature条件写成scalarized max directional derivative控制。
5. HZ beta与非负截断：它解决的是如何在保持K-descent的同时利用前一方向。HZ公式中的负二次项惩罚梯度变化，证明中正是这一项给出beta * Hbar_2 <= ||u_k||^2/(4mu)，从而推出充分下降。beta<0时截断为0，本质上是把不可靠的共轭记忆丢掉。
Key Insight / Why It Works
最重要的insight是：set-valued optimization里的下降可以被压缩成“当前minimal active family上的最大scalarized方向导数”。一旦这个量为负，就不仅是某个标量目标下降，而是对lower set preorder有意义的K-descent。DS scalarization在这里不是普通工具，而是把锥序几何变成Wolfe/CG可操作对象的桥。
方法成立的核心是三层对齐：stationarity residual ||u_k||、K-descent条件、Wolfe curvature条件都使用同一个max_j phi(grad f^{a_j}(x)^T d)。这避免了很多set-valued方法中常见的目标不一致：一个量用于求方向，另一个量用于线搜索，再另一个量用于收敛证明。这里的信息流比较干净。
HZ部分的真正作用是保证下降，而不只是“更快”。证明中关键不是共轭性本身，而是HZ参数的负修正项足以抵消beta d_{k-1}可能带来的正方向导数贡献。换言之，核心贡献更接近better inductive bias for descent preservation，而不是scaling或test-time compute。
但数值增益未必完全来自HZ理论。实验里每步都要求解Step 3子问题，line search和MATLAB solver也会显著影响时间；HZ相对PRP/HS的改善可能有一部分只是特定测试分布、参数设置和solver行为的结果。文中没有做消融来分离这些因素，增益来源不清。
非有限生成锥的推广是实质理论贡献，但实现层面并没有免费午餐。它把“有限极射线枚举”替换为“能计算phi和解相应max-scalarized子问题”。对于K3这种二阶锥型例子phi有闭式，效果自然；对一般锥，计算负担可能成为主要瓶颈。
Relation To Prior Work
最接近的路线有三条：Drummond-Svaiter的vector optimization steepest descent，Pérez-Prudente/Gonçalves-Prudente一系的vector-valued nonlinear CG，和Bouza/Kumar/Ghosh等finite-cardinality set optimization一阶方法。
相对vector optimization CG，这篇的新增信息是把active minimal set和partition引入CG方向构造，使方法不再只处理单个vector objective f(x)，而处理F(x)={f^j(x)}的set comparison。vector optimization对应p=1，是它的特例。
相对Bouza等steepest descent，这篇不是重新定义stationarity，而是在同一类stationarity residual上加入HZ记忆项，试图获得更好的数值效率。理论难点也从“最陡方向存在”转为“递推方向仍然下降”。
相对Kumar et al.和Ghosh et al.的set-valued CG，这篇看似只是补上HZ variant，但HZ不是纯参数替换。HZ的下降证明依赖特定负修正结构，与PRP/HS这类参数的风险控制不同。实质创新在于把HZ保证下降机制移植到set-valued、非有限生成锥、无正则性假设的框架中。
不过，很多构件并非全新：minimal set/partition、V_x stationarity function、DS scalarization、Wolfe/Zoutendijk套路都来自已有工作。论文的贡献更像一次理论拼接与闭合，而不是提出新的set optimization建模范式。
Dataset / Evaluation
实验覆盖了一组标准multiobjective test functions转成set-valued objectives，以及作者构造的非有限生成锥例子。覆盖面包括不同m、n、p、有限生成锥K1/K2和非有限生成锥K3，能说明算法至少在典型合成测试上可运行。
但evaluation主要验证的是“该算法能跑，并且在这些合成问题上常比PRP/HS更快或迭代更少”。它没有真正验证大规模set-valued optimization、复杂非polyhedral cone、真实应用场景或高p情况下的可扩展性。所谓practical effectiveness成立范围有限。
比较对象只选PRP和HS，理由是文献中它们优于FR/CD/DY。这是可以接受的，但没有消融HZ负修正、beta截断、line search参数、Step 3 solver选择，因此无法归因。性能profile支持相对优势，但不支持“HZ机制本身是唯一增益来源”。
实验没有真实世界数据，也没有部署型任务。对一篇数学优化算法论文这不致命，但如果claim扩展到finance/robust optimization等应用，当前证据明显不足。
Limitation
第一，收敛目标只是liminf ||u_k||=0，即渐近stationarity，不是弱极小解收敛，也不是全序列收敛。对于set optimization来说，stationary point和weakly minimal point之间差距可能很大；论文承认stationarity较弱，但实验叙述容易让人误以为得到了解。
第二，复杂度被转移到每步子问题。需要计算Min(F(x),K)、构造partition P_x、求解V_x上的min-max强凸问题。p大、minimal elements多、重复minimal value多时，P_x可能成为实际瓶颈。文中未充分说明这部分如何scale。
第三，非有限生成锥的理论一般性依赖C={w in K*: w^T e=1}的紧性和phi的良好性质，但算法实现需要phi可计算。对于一般闭凸锥，phi本身就是一个support function优化问题；若没有闭式或高效oracle，所谓general cone advantage会变成内层优化负担。
第四，Assumption 4.1/4.2/4.3并不弱。bounded lower level set、Jacobian Lipschitz、以及set sequence存在共同lower bound，这些条件在非凸set-valued问题中可能不容易验证。理论不要求最优点正则性，但仍依赖全局有界性和光滑性结构。
第五，Theorem 4.3的证明写得过于简略，直接称“exactly as in [12, Theorem 2]”。考虑到这里有set-valued minimal index切换、P_x变化和scalarization max结构，文中未充分说明这些差异为何不会影响最后矛盾证明。
第六，数值优势的归因不清。可能主要来自HZ参数，也可能来自特定测试问题、line search策略、solver设置或beta截断导致的稳定性。没有复杂度报告、失败率分析或不同p/m/n scaling曲线。
Takeaway
- 1. 最值得迁移的是“用当前minimal active family上的max scalarized directional derivative统一方向、线搜索和收敛证明”。
- 这个模式比固定标量化更适合set-valued/vector-valued优化。
- 2. 非有限生成锥不一定要靠有限生成近似；用归一化dual generator C和support-function式scalarization可以在理论上绕开polyhedral限制。
- 但工程实现必须配套高效phi oracle。
一句话总结
这篇论文把HZ共轭梯度的保证下降机制移植到finite-cardinality set optimization，并通过Drummond-Svaiter scalarization绕开有限生成锥假设，是一篇偏理论闭合与方法扩展的set-valued first-order optimization工作。