CVXPY 1.0中ZERO与NONNEGATIVE表达式差异问题咨询

这是CVXPY 1.0版本中一个典型的类型推断逻辑bug，和表达式运算顺序、常量类型的处理方式直接相关，我来拆解一下：

为什么两个表达式的标记不同？ #

• 第一个表达式：cv.square(a-b)*1/(1e-6)
  运算顺序是先执行 cv.square(a-b) * 1（整数1和表达式相乘，此时表达式的凸性、符号属性都保持不变），再执行除以1e-6（浮点数）。在CVXPY 1.0的旧版处理逻辑中，除以极小浮点数时，内部的数值检查模块错误地将整个表达式判定为CONSTANT ZERO——这完全是个误判，因为square(a-b)明明是依赖变量a、b的凸表达式。

• 第二个表达式：cv.square(a-b)*float(1/1e-6)
  这里1/1e-6先被Python解释器计算为浮点数1000000.0，再直接和square(a-b)相乘。CVXPY能正确识别这是正标量乘以凸非负表达式，所以保留了CONVEX, NONNEGATIVE的正确属性。

怎么解决这个问题？ #

• 直接合并常量：把乘除运算合并成一个浮点数常量，比如写成 cv.square(a-b) * 1e6，这样CVXPY能准确识别表达式属性。
• 用乘法代替除法：如果需要用倒数，预先计算好浮点数形式的倒数再相乘，避免分步的乘除操作。
• 升级CVXPY版本：这类旧版的类型推断bug在后续的1.1+版本中已经被修复，升级后能避免这类奇怪的问题。

内容的提问来源于stack exchange，提问作者Tony