我正在尝试为三角函数实现范围缩减。 我发现这篇论文http://www.computer.org/csdl/proceedings/pcspa/2010/4180/00/4180b048-abs.html讨论了使用64位整数运算。
提出的想法应该有效，但本文中的方程似乎存在一些问题。 这比fdlibm中实现的效率高吗？
答案 0 :(得分：0)
[编辑]抱歉，我没有将此视为2012年的问题而是2013年的问题。
这不适合作为评论，所以发布在这里。我知道它没有直接回答OP的帖子，但我强烈认为这是该问题的先决条件。
如果您想要执行完整的浮点范围缩减，请咨询K.C.网络上很容易找到Ng的“对于巨大的论据的赞美减少：从头到尾”。
显着的问题是，要对标准三角函数（例如正弦（x））进行范围缩小，其中x是弧度，必须进行涉及Pi的精确mod运算。 mod需要将4 / Pi扩展到足够的派系位以获得有意义的结果。本文详细介绍了这个过程以及需要走多远的过程。事实证明，它是潜在的100位，但不是数百万位。您可能已经意识到这个问题，但如果没有，那么您需要知道使用64位例程或其他任何内容来减少良好。
正如@ eh9所说，你能否在幕后链接上发布一些信息以及你遇到的问题？