草庐IT

c++ - 符合 IEEE-754 标准的四舍五入

coder 2024-01-31 原文

C 标准库提供了 round , lround , 和 llround C99 中的函数族。但是，这些函数不符合 IEEE-754，因为它们没有实现 IEEE 规定的半对偶的“银行家四舍五入”。如果小数部分恰好为 0.5，则半对偶四舍五入要求将结果四舍五入到最接近的偶数。如 cppreference.com 中所述，C99 标准改为要求离零一半。

1-3) Computes the nearest integer value to arg (in floating-point format), rounding halfway cases away from zero, regardless of the current rounding mode.

在 C 中实现舍入的常用方法是表达式 (int)(x + 0.5f)其中，尽管是incorrect在严格的 IEEE-754 数学中，通常由编译器翻译成正确的 cvtss2si操作说明。然而，这当然不是一个可移植的假设。

我如何实现一个函数，该函数将使用半对偶语义对任何浮点值进行舍入？如果可能，该函数应仅依赖于语言和标准库语义，以便它可以在非 IEEE 浮点类型上运行。如果这不可能，根据 IEEE-754 位表示定义的答案也是可以接受的。请根据 <limits.h> 描述任何常数或 <limits> .

最佳答案

The C standard library provides the round, lround, and llround family of functions in C99. However, these functions are not IEEE-754 compliant, because they do not implement the "banker's rounding" of half-to-even as mandated by IEEE...

谈论某个功能是否“符合 IEEE-754”是没有意义的。 IEEE-754 合规性要求一组具有已定义语义的数据类型操作可用。它不要求那些类型或操作具有特定名称，也不要求仅那些操作可用。一个实现可以提供它想要的任何附加功能并且仍然是兼容的。如果一个实现想要提供舍入到奇数、舍入随机、舍入零和不精确的陷阱，它可以这样做。

IEEE-754 对舍入的实际要求是提供以下六种操作:

convertToIntegerTiesToEven(x)

convertToIntegerTowardZero(x)

convertToIntegerTowardPositive(x)

convertToIntegerTowardNegative(x)

convertToIntegerTiesToAway(x)

convertToIntegerExact(x)

在 C 和 C++ 中，这些操作的最后五个绑定(bind)到 trunc、ceil、floor、round 和 rint 函数。 C11 和 C++14 第一个没有绑定(bind)，但 future 的修订将使用 roundeven。如您所见，round 实际上是必需的操作之一。

但是，roundeven 在当前的实现中不可用，这将我们带到您问题的下一部分:

The usual ad-hoc way to implement rounding in C is the expression (int)(x + 0.5f) which, despite being incorrect in strict IEEE-754 math, is usually translated by compilers into the correct cvtss2si instruction. However, this is certainly not a portable assumption.

该表达式的问题远远超出了“严格的 IEEE-754 数学”。它对负 x 完全不正确，对 nextDown(0.5) 给出错误答案，并将 2**23 binade 中的所有奇数整数转换为偶数整数。任何将它翻译成 cvtss2si 的编译器都非常非常糟糕。如果您有发生这种情况的示例，我很乐意看到。

How can I implement a function that will round any floating point value with half-to-even semantics?

正如 njuffa 在评论中指出的那样，您可以确保设置默认舍入模式并使用 rint(或 lrint，因为它听起来您实际上想要一个整数结果)，或者您可以通过调用 round 然后像 gnasher729 建议的那样修复中间情况来实现您自己的舍入函数。一旦采用了 C 的 n1778 绑定(bind)，您就可以使用 roundeven 或 fromfp 函数执行此操作，而无需控制舍入模式。

关于c++ - 符合 IEEE-754 标准的四舍五入，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/32746523/

四舍ampcodestrongblockquotec++cfloating-pointieee-754bankers-rounding

有关c++ - 符合 IEEE-754 标准的四舍五入的更多相关文章

  1. ruby-on-rails - 如何优雅地重启 thin + nginx？ - 2

    我的瘦服务器配置了nginx，我的ROR应用程序正在它们上运行。在我发布代码更新时运行thinrestart会给我的应用程序带来一些停机时间。我试图弄清楚如何优雅地重启正在运行的Thin实例，但找不到好的解决方案。有没有人能做到这一点？ 最佳答案 #Restartjustthethinserverdescribedbythatconfigsudothin-C/etc/thin/mysite.ymlrestartNginx将继续运行并代理请求。如果您将Nginx设置为使用多个上游服务器，例如server{listen80;server

  2. ruby - 将 spawn() 的标准输出/标准错误重定向到 Ruby 中的字符串 - 2

    我想使用spawn(针对多个并发子进程)在Ruby中执行一个外部进程，并将标准输出或标准错误收集到一个字符串中，其方式类似于使用Python的子进程Popen.communicate()可以完成的操作。我尝试将:out/:err重定向到一个新的StringIO对象，但这会生成一个ArgumentError，并且临时重新定义$stdxxx会混淆子进程的输出。 最佳答案 如果你不喜欢popen，这是我的方法:r,w=IO.pipepid=Process.spawn(command,:out=>w,:err=>[:child,:out])

  3. ruby-on-rails - 标准化文件名的字符串，删除重音和特殊字符 - 2

    我正在尝试找到一种方法来规范化字符串以将其作为文件名传递。到目前为止我有这个:my_string.mb_chars.normalize(:kd).gsub(/[^\x00-\x7F]/n,'').downcase.gsub(/[^a-z]/,'_')但第一个问题:-字符。我猜这个方法还有更多问题。我不控制名称，名称字符串可以有重音符、空格和特殊字符。我想删除所有这些，用相应的字母('é'=>'e')替换重音符号，并将其余的替换为'_'字符。名字是这样的:“Prélèvements-常规”“健康证”...我希望它们像一个没有空格/特殊字符的文件名:“prelevements_routin

  4. ruby - 使用 `+=` 和 `send` 方法 - 2

    如何将send与+=一起使用？a=20;a.send"+=",10undefinedmethod`+='for20:Fixnuma=20;a+=10=>30 最佳答案 恐怕你不能。+=不是方法，而是语法糖。参见http://www.ruby-doc.org/docs/ProgrammingRuby/html/tut_expressions.html它说Incommonwithmanyotherlanguages,Rubyhasasyntacticshortcut:a=a+2maybewrittenasa+=2.你能做的最好的事情是:

  5. ruby - 如何计算 Liquid 中的变量 +1 - 2

    我对如何计算通过{%assignvar=0%}赋值的变量加一完全感到困惑。这应该是最简单的任务。到目前为止，这是我尝试过的:{%assignamount=0%}{%forvariantinproduct.variants%}{%assignamount=amount+1%}{%endfor%}Amount:{{amount}}结果总是0。也许我忽略了一些明显的东西。也许有更好的方法。我想要存档的只是获取运行的迭代次数。 最佳答案 因为{{incrementamount}}将输出您的变量值并且不会影响{%assign%}定义的变量，我

  6. arrays - Ruby 数组 += vs 推送 - 2

    我有一个数组数组，想将元素附加到子数组。+=做我想做的，但我想了解为什么push不做。我期望的行为(并与+=一起工作):b=Array.new(3,[])b[0]+=["apple"]b[1]+=["orange"]b[2]+=["frog"]b=>[["苹果"],["橙子"],["Frog"]]通过推送，我将推送的元素附加到每个子数组(为什么？):a=Array.new(3,[])a[0].push("apple")a[1].push("orange")a[2].push("frog")a=>[[“苹果”、“橙子”、“Frog”]、[“苹果”、“橙子”、“Frog”]、[“苹果”、“

  7. Ruby:标准递归模式 - 2

    我经常迷上ruby​​的一件事是递归模式。例如，假设我有一个数组，它可能包含无限深度的数组作为元素。所以，例如:my_array=[1,[2,3,[4,5,[6,7]]]]我想创建一个方法，可以将数组展平为[1,2,3,4,5,6,7]。我知道.flatten可以完成这项工作，但这个问题是作为我经常遇到的递归问题的一个例子-因此我试图找到一个更可重用的解决方案。简而言之-我猜这种事情有一个标准模式，但我想不出任何特别优雅的东西。任何想法表示赞赏 最佳答案 递归是一种方法，它不依赖于语言。您在编写算法时要考虑两种情况:再次调用函数的情

  8. ruby-on-rails - 使用 Ruby 标准 Logger 每天只创建一个日志 - 2

    我正在使用ruby​​标准记录器，我想要每天轮换一次，所以在我的代码中我有:Logger.new("#{$ROOT_PATH}/log/errors.log",'daily')它运行完美，但它创建了两个文件errors.log.20130217和errors.log.20130217.1。如何强制它每天只创建一个文件？ 最佳答案 您的代码对于长时间运行的应用程序是正确的。发生的事情是您在给定的一天多次运行代码。第一次运行时，Ruby会创建一个日志文件“errors.log”。当日期改变时，Ruby将文件重命名为“errors.log

  9. += 的 Ruby 方法 - 2

    有没有办法让Ruby能够做这样的事情？classPlane@moved=0@x=0defx+=(v)#thisiserror@x+=v@moved+=1enddefto_s"moved#{@moved}times,currentxis#{@x}"endendplane=Plane.newplane.x+=5plane.x+=10putsplane.to_s#moved2times,currentxis15 最佳答案 您不能在Ruby中覆盖复合赋值运算符。任务在内部处理。您应该覆盖+，而不是+=。plane.a+=b与plane.a=

  10. ruby - Sinatra + Heroku + Datamapper 使用 dm-sqlite-adapter 部署问题 - 2

    出于某种原因，heroku尝试要求dm-sqlite-adapter，即使它应该在这里使用Postgres。请注意，这发生在我打开任何URL时-而不是在gitpush本身期间。我构建了一个默认的Facebook应用程序。gem文件:source:gemcuttergem"foreman"gem"sinatra"gem"mogli"gem"json"gem"httparty"gem"thin"gem"data_mapper"gem"heroku"group:productiondogem"pg"gem"dm-postgres-adapter"endgroup:development,:t

随机推荐