[论文解读] A Robust Complex Division in Scilab

Michaël Baudin, Robert L. Smith|arXiv (Cornell University)|Oct 16, 2012

Microwave Engineering and Waveguides参考文献 5被引用 5

一句话总结

本文提出了一种简化且鲁棒的双精度浮点数复数除法算法，通过避免中间计算中的溢出和下溢，改进了Smith于1962年提出的方法。作者证明了该算法的鲁棒性，并通过数值模拟表明，其可靠性显著优于Smith的方法，且与Stewart更复杂的算法相当，失败率比简单实现低多达四个数量级。

ABSTRACT

The most widely used algorithm for floating point complex division, known as Smith's method, may fail more often than expected. This document presents two improved complex division algorithms. We present a proof of the robustness of the first improved algorithm. Numerical simulations show that this algorithm performs well in practice and is significantly more robust than other known implementations. By combining additionnal scaling methods with this first algorithm, we were able to create a second algorithm, which rarely fails.

研究动机与目标

• 解决Smith复数除法算法因中间项如$c^2 + d^2$导致溢出和下溢而频繁失效的问题。
• 开发一种比Stewart的复数乘积计算方法更简单但同样鲁棒的替代方案，后者需要测试三种不同的三元乘积计算方式。
• 提供一个可证明鲁棒的算法，即使在接近下溢阈值时也能保持高精度。
• 通过大量随机测试，评估并比较所提算法与现有方法在数值鲁棒性和性能方面的表现。

提出的方法

• 改进的算法采用一种修改后的公式，通过重新组织复数除法表达式，避免了有问题的$c^2 + d^2$项，从而降低中间结果溢出和下溢的风险。
• 根据大小选择三元乘积的最佳乘法顺序，消除其中一种不必要的排列方式，从而简化实现。
• 算法检测下溢状态，并切换到更安全的表达形式，例如$e = a*t + (b*t)*r$，其中两项可能各自下溢，但其和仍保持准确。
• 第二种更鲁棒的变体结合了缩放技术与改进算法，用于处理结果接近下溢阈值$\alpha$的情况。
• 作者在ANSI C和Scilab中实现了所有算法，并使用包含1,574,802个随机复数除法的测试集来衡量性能与精度。
• 通过随机数值实验评估鲁棒性，重点关注失败率和结果的数字精度。

实验结果

研究问题

• RQ1 是否可以设计一种更简单的复数除法算法，在无需计算三元乘积所有三种方式的前提下，实现与Stewart方法相当的鲁棒性？
• RQ2 Smith方法在实际中失败的频率如何？这些失败是否可以系统性地识别并缓解？
• RQ3 在接近下溢阈值时，缩放和表达式重组能在多大程度上提升复数除法的精度？
• RQ4 增加鲁棒性带来的性能开销有多大？在实际应用中是否可接受？

主要发现

• 与简单实现相比，改进算法的失败率降低了多达四个数量级，与其它已知实现相比降低了两个数量级。
• 改进算法的鲁棒性与Stewart方法相当，其失败率显著低于Smith、Priest、Li等人以及C99的算法。
• 该算法的鲁棒变体在大多数情况下可产生超过52位有效数字的实部和虚部，即使结果接近下溢阈值$\alpha$。
• 性能测量显示，编译后的改进算法运行速度为5.0 MCDPS，略快于Scilab原生的4.7 MCDPS；而鲁棒版本为3.2 MCDPS，反映出其更高的复杂度。
• 算法能够检测中间项的下溢，并切换到更安全的表达式，即使各项本身为 subnormal 也能保持精度。
• 数值实验确认，该算法仅在极少数情况下失败，即当精确结果极其接近下溢阈值时，即便如此，其模长仍能被正确保留。