CVXPY求解器在1.4.3版本后出现简单问题求解错误分析

问题背景

CVXPY是一个用于凸优化问题建模和求解的Python库。近期在版本升级过程中，用户发现了一个值得注意的数值求解问题：从1.4.3版本开始，多个求解器（PIQP、CLARABEL和OSQP）在处理特定形式的优化问题时会出现错误结果，而ECOS求解器则表现正常。

问题现象

考虑一个简单的投资组合优化问题，目标函数包含三部分：

1. 资产预期收益（线性项）
2. 风险惩罚项（二次项）
3. 交易成本（绝对值项）

当风险惩罚项采用元素级平方和形式（即0.02*(cp.multiply(w,ivol)**2).sum()）时，CVXPY 1.4.2版本各求解器结果一致，但1.4.3及以后版本中，除ECOS外其他求解器都给出了明显错误的解（约817 vs 正确的1384）。

技术分析

这个问题揭示了几个关键点：

1. 表达式转换差异：CVXPY在内部会将用户定义的优化问题转换为标准形式。1.4.3版本可能在表达式转换过程中对元素级平方运算的处理发生了变化。

2. 求解器敏感性：不同求解器对问题表述的敏感度不同。ECOS作为成熟的二阶锥求解器，可能对这类问题的数值稳定性更好。

3. 替代方案验证：当风险项改用标准的二次型形式（cp.quad_form(w,sigma)）时，所有求解器都能给出正确结果，这说明问题特定于元素级平方运算的处理。

解决方案

开发团队已经确认这是一个bug，并将在近期发布修复版本。对于遇到类似问题的用户，建议：

1. 暂时使用ECOS求解器
2. 或者将风险项改写为标准二次型形式
3. 等待官方发布的修复版本

对用户的建议

这类问题在数值优化中并不罕见，特别是当涉及不同形式的数学表达式转换时。建议用户：

1. 对于关键优化问题，使用多个求解器交叉验证结果
2. 注意记录使用的软件版本号
3. 关注表达式形式对求解精度的影响
4. 及时更新到稳定版本

数值优化问题的正确性验证是实际应用中不可忽视的环节，特别是在金融、工程等关键领域。这个案例也展示了开源社区快速响应和修复问题的优势。