我想知道这是不是真的:当我对一个平方整数求平方根时,就像在
f = Math.sqrt(123*123)
我将得到一个非常接近 123 的 float 。由于浮点表示精度,这可能类似于 122.99999999999999999999 或 123.000000000000000000001。
因为 floor(122.999999999999999999) 是 122,我应该得到 122 而不是 123。所以我希望 floor(sqrt(i*i)) == i-1 在大约 50% 的情况下。奇怪的是,对于我测试过的所有数字,floor(sqrt(i*i) == i。这是一个用于测试前 1 亿个数字的小 ruby 脚本:
100_000_000.times do |i|
puts i if Math.sqrt(i*i).floor != i
end
上面的脚本从不打印任何东西。为什么?
更新:感谢您的快速回复,这似乎是解决方案:根据 wikipedia
Any integer with absolute value less than or equal to 2^24 can be exactly represented in the single precision format, and any integer with absolute value less than or equal to 2^53 can be exactly represented in the double precision format.
Math.sqrt(i*i) 从 i=9007199254740993 开始按照我的预期运行,即 2^53 + 1。
最佳答案
这是你困惑的本质:
Due to floating point representation precision, this could be something like 122.99999999999999999999 or 123.000000000000000000001.
这是错误的。在符合 IEEE-754 标准的系统上,它总是恰好是 123,这几乎是现代的所有系统。浮点运算没有“随机误差”或“噪声”。它具有精确的、确定性的舍入,许多简单的计算(比如这个)根本不会产生任何舍入。
123 可以精确表示为 float ,123*123 也是如此(所有 中等大小的整数也是如此)。因此当您将 123*123 转换为浮点类型时不会发生舍入错误。结果完全 15129。
根据 IEEE-754 标准,平方根是正确舍入的运算。这意味着如果有一个确切的答案,则需要平方根函数来产生它。由于您要对 exactly 15129 求平方根,即 exactly 123,因此 exactly 从平方根函数得到的结果。不进行舍入或近似。
现在,对于多大的整数,这将是正确的?
double 可以准确表示 2^53 以内的所有整数。因此,只要 i*i 小于 2^53,您的计算中就不会发生舍入,因此结果将是准确的。这意味着对于所有小于 94906265 的 i,我们知道计算是准确的。
但是你试过 i 比那个大!发生了什么事?
对于您尝试过的最大 i,i*i 仅略大于 2^53 (1.1102... * 2^53,实际上)。因为从整数到 double 的转换(或 double 乘法)也是正确舍入操作,i*i 将是最接近 i 精确平方的可表示值。在这种情况下,由于 i*i 是 54 位宽,舍入将发生在最低位。因此我们知道:
i*i as a double = the exact value of i*i + rounding
其中 rounding 是 -1,0 或 1。如果四舍五入为零,则平方是精确的,因此平方根是精确的,所以我们已经知道您得到了正确答案。让我们忽略这种情况。
现在我们要计算 i*i +/- 1 的平方根。使用泰勒级数展开,此平方根的无限精确(未四舍五入)值为:
i * (1 +/- 1/(2i^2) + O(1/i^4))
如果你以前没有做过任何浮点错误分析,现在这有点繁琐,但如果你使用 i^2 > 2^53 这个事实,你可以看到的:
1/(2i^2) + O(1/i^4)
term 小于 2^-54,这意味着(因为平方根被正确舍入,因此它的舍入误差必须小于 2^54),舍入 sqrt 的结果函数正是 i。
事实证明(通过类似的分析),对于任何可精确表示的 float x,sqrt(x*x) 正好是 x(假设 x*x 的中间计算没有不会溢出或下溢),因此您遇到此类计算舍入的唯一方法是在 x 本身的表示中,这就是为什么您看到它从 2^53 开始+ 1(不可表示的最小整数)。
关于ruby - 为什么 Math.sqrt(i*i).floor == i?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/2060486/
我正在学习如何使用Nokogiri,根据这段代码我遇到了一些问题:require'rubygems'require'mechanize'post_agent=WWW::Mechanize.newpost_page=post_agent.get('http://www.vbulletin.org/forum/showthread.php?t=230708')puts"\nabsolutepathwithtbodygivesnil"putspost_page.parser.xpath('/html/body/div/div/div/div/div/table/tbody/tr/td/div
总的来说,我对ruby还比较陌生,我正在为我正在创建的对象编写一些rspec测试用例。许多测试用例都非常基础,我只是想确保正确填充和返回值。我想知道是否有办法使用循环结构来执行此操作。不必为我要测试的每个方法都设置一个assertEquals。例如:describeitem,"TestingtheItem"doit"willhaveanullvaluetostart"doitem=Item.new#HereIcoulddotheitem.name.shouldbe_nil#thenIcoulddoitem.category.shouldbe_nilendend但我想要一些方法来使用
我有一个字符串input="maybe(thisis|thatwas)some((nice|ugly)(day|night)|(strange(weather|time)))"Ruby中解析该字符串的最佳方法是什么?我的意思是脚本应该能够像这样构建句子:maybethisissomeuglynightmaybethatwassomenicenightmaybethiswassomestrangetime等等,你明白了......我应该一个字符一个字符地读取字符串并构建一个带有堆栈的状态机来存储括号值以供以后计算,还是有更好的方法?也许为此目的准备了一个开箱即用的库?
我有一个Ruby程序,它使用rubyzip压缩XML文件的目录树。gem。我的问题是文件开始变得很重,我想提高压缩级别,因为压缩时间不是问题。我在rubyzipdocumentation中找不到一种为创建的ZIP文件指定压缩级别的方法。有人知道如何更改此设置吗?是否有另一个允许指定压缩级别的Ruby库? 最佳答案 这是我通过查看rubyzip内部创建的代码。level=Zlib::BEST_COMPRESSIONZip::ZipOutputStream.open(zip_file)do|zip|Dir.glob("**/*")d
类classAprivatedeffooputs:fooendpublicdefbarputs:barendprivatedefzimputs:zimendprotecteddefdibputs:dibendendA的实例a=A.new测试a.foorescueputs:faila.barrescueputs:faila.zimrescueputs:faila.dibrescueputs:faila.gazrescueputs:fail测试输出failbarfailfailfail.发送测试[:foo,:bar,:zim,:dib,:gaz].each{|m|a.send(m)resc
很好奇,就使用rubyonrails自动化单元测试而言,你们正在做什么?您是否创建了一个脚本来在cron中运行rake作业并将结果邮寄给您?git中的预提交Hook?只是手动调用?我完全理解测试,但想知道在错误发生之前捕获错误的最佳实践是什么。让我们理所当然地认为测试本身是完美无缺的,并且可以正常工作。下一步是什么以确保他们在正确的时间将可能有害的结果传达给您? 最佳答案 不确定您到底想听什么,但是有几个级别的自动代码库控制:在处理某项功能时,您可以使用类似autotest的内容获得关于哪些有效,哪些无效的即时反馈。要确保您的提
这似乎应该有一个直截了当的答案,但在Google上花了很多时间,所以我找不到它。这可能是缺少正确关键字的情况。在我的RoR应用程序中,我有几个模型共享一种特定类型的字符串属性,该属性具有特殊验证和其他功能。我能想到的最接近的类似示例是表示URL的字符串。这会导致模型中出现大量重复(甚至单元测试中会出现更多重复),但我不确定如何让它更DRY。我能想到几个可能的方向...按照“validates_url_format_of”插件,但这只会让验证干给这个特殊的字符串它自己的模型,但这看起来很像重溶液为这个特殊的字符串创建一个ruby类,但是我如何得到ActiveRecord关联这个类模型
假设我做了一个模块如下:m=Module.newdoclassCendend三个问题:除了对m的引用之外,还有什么方法可以访问C和m中的其他内容?我可以在创建匿名模块后为其命名吗(就像我输入“module...”一样)?如何在使用完匿名模块后将其删除,使其定义的常量不再存在? 最佳答案 三个答案:是的,使用ObjectSpace.此代码使c引用你的类(class)C不引用m:c=nilObjectSpace.each_object{|obj|c=objif(Class===objandobj.name=~/::C$/)}当然这取决于
我试图在一个项目中使用rake,如果我把所有东西都放到Rakefile中,它会很大并且很难读取/找到东西,所以我试着将每个命名空间放在lib/rake中它自己的文件中,我添加了这个到我的rake文件的顶部:Dir['#{File.dirname(__FILE__)}/lib/rake/*.rake'].map{|f|requiref}它加载文件没问题,但没有任务。我现在只有一个.rake文件作为测试,名为“servers.rake”,它看起来像这样:namespace:serverdotask:testdoputs"test"endend所以当我运行rakeserver:testid时
出于纯粹的兴趣,我很好奇如何按顺序创建PI,而不是在过程结果之后生成数字,而是让数字在过程本身生成时显示。如果是这种情况,那么数字可以自行产生,我可以对以前看到的数字实现垃圾收集,从而创建一个无限系列。结果只是在Pi系列之后每秒生成一个数字。这是我通过互联网筛选的结果:这是流行的计算机友好算法,类机器算法:defarccot(x,unity)xpow=unity/xn=1sign=1sum=0loopdoterm=xpow/nbreakifterm==0sum+=sign*(xpow/n)xpow/=x*xn+=2sign=-signendsumenddefcalc_pi(digits