APL中*○0j1与*(○0j1)输出差异原因及区别问询

核心区别 #

两个表达式的本质差异在于APL解释器的计算优化策略：

• *○0j1：被解释器识别为经典数学等式的特例，执行符号级精确计算，直接返回结果-1。
• *(○0j1)：按分步计算逻辑执行，引入浮点数精度误差，最终返回带极小虚部的近似值¯1j1e¯16。

详细解释 #

1. APL函数的基础行为

   • ○函数的单参数语义：○x返回x×π，因此○0j1的数学理论值是0jπ。
   • *函数的语义：*z表示自然指数函数e^z，对于复数z=a+bj，等价于e^a × (cosb + i sinb)。

2. *○0j1的精确计算逻辑
   APL解释器会对常见的复合表达式做优化，*○0j1对应数学中的欧拉公式特例e^(iπ) = -1，解释器直接匹配这个规则，跳过浮点数的近似计算，直接返回精确值-1。

3. *(○0j1)的近似计算逻辑
   这个表达式被拆分为两步执行：

   • 第一步：计算○0j1，由于π无法用浮点数精确表示，得到的是0j3.141592653589793（π的双精度浮点数近似值）。
   • 第二步：计算*作用于这个复数，即e^(0j3.141592653589793)。因为这个角度略小于精确的π，sin(3.141592653589793)会得到一个约1e-16的极小非零值，最终结果就是-1 + 1e-16j，即APL表示的¯1j1e¯16。

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