为何Numpy complex128做NaN除法告警，而内置complex不会？

你遇到的这个差异，本质是Numpy对实数和复数NaN运算的警告判定逻辑不同，以及和Python内置类型的实现差异导致的：

1. Python内置complex的处理逻辑
   Python内置的复数类型在处理NaN/(1j*NaN)这类运算时，会直接返回nan+nanj且不会触发警告。这是因为Python的复数运算逻辑将NaN参与的运算结果视为合法的NaN值，不会额外触发错误提示。

2. Numpy float64的处理逻辑
   对于Numpy的float64类型，NaN/NaN的结果是NaN，这属于“安静NaN（quiet NaN）”的生成操作，Numpy将这类操作视为合法的数值运算，因此不会触发警告。

3. Numpy complex128的特殊处理
   Numpy的complex128运算底层基于C语言的复数运算实现，处理复数除法时会分别计算结果的实部和虚部。当实部或虚部出现NaN/NaN的情况时，Numpy的错误检查机制会将其判定为“无效操作（invalid operation）”，从而抛出RuntimeWarning。尽管最终结果同样是nan+nanj，但Numpy在复数类型的运算校验上比实数类型更严格，这种警告触发逻辑是Numpy复数运算模块的特定实现导致的。

需要注意的是，这种行为可能和你使用的Numpy 1.24.2版本相关，后续版本可能对这类警告逻辑做了调整，但核心原因还是复数与实数运算在Numpy内部的错误判定规则不一致。

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