草庐IT

c++ - 浮点除法的软件实现，舍入问题

coder 2023-11-15 原文

作为一个学习项目，我正在使用 C++ 在软件中实现浮点运算(add、sub、mul、div)。目标是更加熟悉浮点行为的底层细节。

我正在尝试将我的处理器操作与精确位匹配，即 IEEE 754 标准。到目前为止，它一直运行良好，add、sub 和 mult 表现完美，我对它进行了大约 1.1 亿次随机操作测试，并得到了与处理器在硬件中所做的完全相同的结果。 (虽然没有考虑边缘情况、溢出等)。

之后，我开始移动到最后一个操作，除法。它工作正常并达到了想要的结果，但有时我会把最后一个尾数弄错，而不是四舍五入。我很难理解为什么。 我一直在使用的主要引用是来自 John Farrier 的精彩演讲(时间戳位于显示如何舍入的位置):

https://youtu.be/k12BJGSc2Nc?t=1153

四舍五入对所有操作都非常有效，但给我的部门带来了麻烦。 让我给你一个具体的例子。 我正在尝试将 645.68011474609375 除以 493.20962524414063

我得到的最终结果是:

我的:0-01111111-01001111001000111100000

c++_:0-01111111-01001111001000111100001

如您所见，除最后一位外，所有内容都匹配。我计算除法的方式是基于这个视频: https://www.youtube.com/watch?v=fi8A4zz1d-s

在此之后，我计算了尾数精度 24 的 28 位(隐藏 1 + 23 尾数)和守卫的 3 位，圆形粘性加上可能移位的额外一位。 使用视频的算法，我最多可以获得 1 的归一化偏移，这就是为什么我在末尾有一个额外的位以防在归一化中移入，所以在舍入中可用。现在这是我从除法算法中得到的结果:

 010100111100100011110000 0100
 ------------------------ ----
 ^                        grs^
 |__ to be normalized        |____ extra bit

如您所见，我在第 24 个位置得到了 0，因此我需要向左移动一位以获得正确的归一化。 这意味着我将得到:

10100111100100011110000 100

根据 John Farrier 的视频，在 100 grs 位的情况下，我仅在尾数的 LSB 为 1 时才进行归一化。在我的情况下为零，这就是为什么我没有对我的结果进行四舍五入。

我有点迷茫的原因是我确信我的算法正在计算正确的尾数，我已经用在线计算器仔细检查过它，舍入策略适用于所有其他操作。此外，以这种方式计算会触发归一化，最终产生正确的指数。

我错过了什么吗？某个地方的小细节？

让我感到奇怪的一件事是粘滞位，在加法和乘法中你会得到不同程度的移位，这导致粘滞位触发的可能性更高，在这种情况下，我只移动一个使粘性位不是真正粘性的最大值。

我确实希望我提供了足够的细节以使我的问题得到理解。在这里，您可以在底部找到我的除法实现，其中充满了我用于调试的打印，但应该让您了解我在做什么，代码从第 374 行开始:

https://gist.github.com/giordi91/1388504fadcf94b3f6f42103dfd1f938

PS:与此同时，我正在浏览“科学家应该知道的关于 float 的一切”，看看我是否遗漏了什么。

最佳答案

您从除法算法中得到的结果不充分。你显示:

 010100111100100011110000 0100
 ------------------------ ----
 ^                        grs^
 |__ to be normalized        |____ extra bit

数学上精确的商继续:

 010100111100100011110000 0100 110000111100100100011110…

因此，四舍五入处的余数超过 ½ ULP，因此应四舍五入。我没有详细研究您的代码，但看起来您可能刚刚计算了一个或两个有效数1。您实际上需要知道残差是非零的，而不仅仅是它的下一位或两位是否为零。如果精确数学结果中位于或超出该位置的位中的任何非零，则最终粘性位应为 1。

脚注

1 “有效数”是首选术语。 “尾数”是对数的小数部分的旧术语。浮点值的有效数字是线性的。尾数是对数。

关于c++ - 浮点除法的软件实现，舍入问题，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/47909541/

amp43四舍strong的c++algorithmfloating-point

有关c++ - 浮点除法的软件实现，舍入问题的更多相关文章

  1. ruby - 在 64 位 Snow Leopard 上使用 rvm、postgres 9.0、ruby 1.9.2-p136 安装 pg gem 时出现问题 - 2

    我想为Heroku构建一个Rails3应用程序。他们使用Postgres作为他们的数据库，所以我通过MacPorts安装了postgres9.0。现在我需要一个postgresgem并且共识是出于性能原因你想要pggem。但是我对我得到的错误感到非常困惑当我尝试在rvm下通过geminstall安装pg时。我已经非常明确地指定了所有postgres目录的位置可以找到但仍然无法完成安装:$envARCHFLAGS='-archx86_64'geminstallpg--\--with-pg-config=/opt/local/var/db/postgresql90/defaultdb/po

  2. ruby - 通过 rvm 升级 ruby​​gems 的问题 - 2

    尝试通过RVM将RubyGems升级到版本1.8.10并出现此错误:$rvmrubygemslatestRemovingoldRubygemsfiles...Installingrubygems-1.8.10forruby-1.9.2-p180...ERROR:Errorrunning'GEM_PATH="/Users/foo/.rvm/gems/ruby-1.9.2-p180:/Users/foo/.rvm/gems/ruby-1.9.2-p180@global:/Users/foo/.rvm/gems/ruby-1.9.2-p180:/Users/foo/.rvm/gems/rub

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

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

  4. ruby - 通过 RVM (OSX Mountain Lion) 安装 Ruby 2.0.0-p247 时遇到问题 - 2

    我的最终目标是安装当前版本的RubyonRails。我在OSXMountainLion上运行。到目前为止，这是我的过程:已安装的RVM$\curl-Lhttps://get.rvm.io|bash-sstable检查已知(我假设已批准)安装$rvmlistknown我看到当前的稳定版本可用[ruby-]2.0.0[-p247]输入命令安装$rvminstall2.0.0-p247注意:我也试过这些安装命令$rvminstallruby-2.0.0-p247$rvminstallruby=2.0.0-p247我很快就无处可去了。结果:$rvminstall2.0.0-p247Search

  5. ruby - Fast-stemmer 安装问题 - 2

    由于fast-stemmer的问题，我很难安装我想要的任何ruby​​gem。我把我得到的错误放在下面。Buildingnativeextensions.Thiscouldtakeawhile...ERROR:Errorinstallingfast-stemmer:ERROR:Failedtobuildgemnativeextension./System/Library/Frameworks/Ruby.framework/Versions/2.0/usr/bin/rubyextconf.rbcreatingMakefilemake"DESTDIR="cleanmake"DESTDIR=

  6. ruby - 如何根据特征实现 FactoryGirl 的条件行为 - 2

    我有一个用户工厂。我希望默认情况下确认用户。但是鉴于unconfirmed特征，我不希望它们被确认。虽然我有一个基于实现细节而不是抽象的工作实现，但我想知道如何正确地做到这一点。factory:userdoafter(:create)do|user,evaluator|#unwantedimplementationdetailshereunlessFactoryGirl.factories[:user].defined_traits.map(&:name).include?(:unconfirmed)user.confirm!endendtrait:unconfirmeddoenden

  7. ruby - 安装 Ruby 时遇到问题(无法下载资源 "readline--patch") - 2

    当我尝试安装Ruby时遇到此错误。我试过查看this和this但无济于事➜~brewinstallrubyWarning:YouareusingOSX10.12.Wedonotprovidesupportforthispre-releaseversion.Youmayencounterbuildfailuresorotherbreakages.Pleasecreatepull-requestsinsteadoffilingissues.==>Installingdependenciesforruby:readline,libyaml,makedepend==>Installingrub

  8. java - 从 JRuby 调用 Java 类的问题 - 2

    我正在尝试使用boilerpipe来自JRuby。我看过guide从JRuby调用Java，并成功地将它与另一个Java包一起使用，但无法弄清楚为什么同样的东西不能用于boilerpipe。我正在尝试基本上从JRuby中执行与此Java等效的操作:URLurl=newURL("http://www.example.com/some-location/index.html");Stringtext=ArticleExtractor.INSTANCE.getText(url);在JRuby中试过这个:require'java'url=java.net.URL.new("http://www

  9. ruby-on-rails - 简单的 Ruby on Rails 问题——如何将评论附加到用户和文章？ - 2

    我意识到这可能是一个非常基本的问题，但我现在已经花了几天时间回过头来解决这个问题，但出于某种原因，Google就是没有帮助我。(我认为部分问题在于我是一个初学者，我不知道该问什么......)我也看过O'Reilly的RubyCookbook和RailsAPI，但我仍然停留在这个问题上.我找到了一些关于多态关系的信息，但它似乎不是我需要的(尽管如果我错了请告诉我)。我正在尝试调整MichaelHartl'stutorial创建一个包含用户、文章和评论的博客应用程序(不使用脚手架)。我希望评论既属于用户又属于文章。我的主要问题是:我不知道如何将当前文章的ID放入评论Controller。

  10. 【高数】用拉格朗日中值定理解决极限问题 - 2

    首先回顾一下拉格朗日定理的内容：函数f(x)是在闭区间[a,b]上连续、开区间(a,b)上可导的函数，那么至少存在一个，使得:通过这个表达式我们可以知道，f(x)是函数的主体，a和b可以看作是主体函数f(x)中所取的两个值。那么可以有，  也就意味着我们可以用来替换 这种替换可以用在求某些多项式差的极限中。方法： 外层函数f(x)是一致的，并且h(x)和g(x)是等价无穷小。此时，利用拉格朗日定理，将原式替换为 ，再进行求解，往往会省去复合函数求极限的很多麻烦。使用要注意：1.要先找到主体函数f(x)，即外层函数必须相同。2.f(x)找到后，复合部分是等价无穷小。3.要满足作差的形式。如果是加

随机推荐