Talk:平方根倒数速算法

message: **:由英文GA翻译而得,很有趣的文章(可当都市传说(?)看),准备试着冲GA,欢迎顺便给些改进意见. - Dr. Cravix ♬La Pluie 2012年5月15日 (二) 15:14 (UTC)
**有意思的條目，不過看的有點累--Wolfch (留言) 2012年5月15日 (二) 15:36 (UTC)
**，可以申请GA。--CHEM.is.TRY 2012年5月16日 (三) 01:27 (UTC)
**：有趣的条目。 AlexHe34（留言） 2012年5月16日 (三) 09:46 (UTC)
**--Amazingloong ☏ 2012年5月16日 (三) 10:01 (UTC)
**--yrr933! (留言)|化學、科技提升計畫! 2012年5月16日 (三) 12:03 (UTC)
**--Alexchris（留言） 2012年5月17日 (四) 00:46 (UTC)
**--B2322858（留言） 2012年5月17日 (四) 02:42 (UTC)
**：摩卡·賀昇 2012年5月18日 (五) 08:28 (UTC)
**--Iflwlou [ M {  2012年5月18日 (五) 13:50 (UTC)

message: ，分类：电脑信息-算法，提名人： Dr. Cravix ♬La Pluie 2012年5月21日 (一) 08:24 (UTC)
replies: text: 投票期：2012年5月21日 (一) 08:24 (UTC) 至 2012年5月28日 (一) 08:24 (UTC)
*：提名人票。个人翻译自英文优良条目,dedicated to id Software.BTW,看了一下没有"算法"这个分类,如果没办法就直接和enwiki一样归进数学类吧.-- Dr. Cravix ♬La Pluie 2012年5月21日 (一) 08:24 (UTC)
**：顺便借问:根据FA的标准此条目若要入FA的话还需作何修改呢?(虽然最近也会比较忙而没空马上改,囧) - Dr. Cravix ♬La Pluie 2012年5月21日 (一) 23:51 (UTC)
*：深知编程算开方之麻烦没效率，然而居然会有这种算法……条目素质甚佳，叙述很清楚易懂，赞一个。关于没有算法分类我只想引用该条目代码概述第十行之注释--Hans Li（Li|Sn|ΔH） 2012年5月21日 (一) 08:57 (UTC)
*雷神之锤xD--铁铁的火大了（留言） 2012年5月21日 (一) 10:35 (UTC)
*--有意思，参考很多足以支持全文。Merphisto（留言） 2012年5月22日 (二) 00:49 (UTC)
*题外话，MathJax加载慢死了……--达师 － 218 － 372 2012年5月22日 (二) 06:44 (UTC)
**Firefox下速度正常,会不会是浏览器问题? - Dr. Cravix ♬La Pluie 2012年5月23日 (三) 03:30 (UTC)
*，比较专业，也很详细--Huandy618 (留言) 2012年5月22日 (二) 13:53 (UTC)
*--有意思，长见识了。Shypanda（留言） 2012年5月23日 (三) 02:27 (UTC)
*：注释部分参考文献不足。那些灰色背景的公式为什么要弄灰色背景？有一些多余的粗体。“魔术数字”部分有个公式显示错误---{百無一用是書生}- (☎) 2012年5月23日 (三) 02:30 (UTC)
**:嗯,谢谢建议与修正注释的确没写清参考,于是加了;灰色背景是标明界限利于与正文相区分;多余的粗体是指什么呢?全文除表格外只有三处粗体呀(刚才顺便修掉了一处);我这里公式显示正常,还是说我公式里写错什么了呢?如果是的话还请帮我修改一下吧 - Dr. Cravix ♬La Pluie 2012年5月23日 (三) 03:30 (UTC)
*，翻译通顺。PS:为毛俺的Fx就不解析MathJax……--铁铁的火大了（留言） 2012年5月23日 (三) 06:35 (UTC)
**只能说我的Arch下Firefox,Opera,Midori甚至虚拟机XP下的IE都可以解析是版本问题吗?还是如果可以的话请描述一下具体的问题吧,我试着检查一下~  - Dr. Cravix ♬La Pluie 2012年5月23日 (三) 09:02 (UTC)
***呃，全是灰色的字，写的全是线性排列的公式。Fx4和Fx12都不能解析，神奇……--铁铁的火大了（留言） 2012年5月23日 (三) 14:14 (UTC)
****这个...是所有公式都显示成原始代码吗?也看看Help:数学公式,样例都能正常显示吗?如果不行的话我的确搞不清楚是怎么回事了= = - Dr. Cravix ♬La Pluie 2012年5月23日 (三) 23:48 (UTC)
*****呃，我还是接着用图片吧，公式能显示出来就行了……--铁铁的火大了（留言） 2012年5月24日 (四) 04:48 (UTC)
******: 囧,说到图片我才明白是怎么回事,我都没注意我本来设置的就是显示为png= =刚才试了一下MathJax,也是可以的,但是流程是:页面载入完后(一定要载入完不能中断)先显示为灰色的原式(应该就是你看到的),然后左下角显示 载入MathJax(一个小框里面显示"Loading (MathJax各组件)"),载入完就能正常显示了~要不你再试试吧 - Dr. Cravix ♬La Pluie 2012年5月24日 (四) 05:31 (UTC)
*，内容有趣--CHEM.is.TRY 2012年5月23日 (三) 10:46 (UTC)
*，连我都能看懂，说明这语句够通俗了。-- 豆腐daveduv编写上海留言  2012年5月24日 (四) 15:59 (UTC)
*：內容豐富详细，參考資料足以支撐全文，段落大致上都有注腳，以支持票作獎勵。—ArikamaI 在沒有人有槍的國度裡，一把手槍的人就是國王（謝絕廢話|戰鬥記錄） 2012年5月25日 (五) 10:31 (UTC)
thread_title:  优良条目候选 

message: 原来的那堆又是rho又是M的证明，相当复杂，并且难以理解。还不如我给添加上去的那段解释，实际上等于后面
的详细解释，并较其严谨和易懂。其实关于尾数部分的选择，我记得当年也看到过一篇文章详细介绍其思路的，这里没有这部分的记录，较为可惜。当年我还仔细推敲过，这个尾数的选择算是非常巧妙。可惜我已经不记得了（包括来源），无法贡献了。如果有朋友记得，还请修改一下。—— Sumtec（赞美 骂街 讨论 察看贡献） 2012年6月8日 (五) 09:14 (UTC)
replies: text: 我觉得指数如何被-2除是很容易理解的事情(参看上面GAN存档的理由),浮点数部分也没必要举例(前面的说明已经很清楚了,再详细有点越俎代庖代替浮点条目之嫌...个人感觉吧),但这且按下不提,你这样直接插一段让我非常难办,因为内容明显是重复的,条理也被打乱了,现在一时间没有整理的头绪.这个条目毕竟已经是GA了,扩充内容固然好,但修改前还是想想怎么才能比较保持条理吧...现在只作了校对,剩下的后面再说. - Dr. Cravix ♬La Pluie 2012年6月8日 (五) 14:29 (UTC)
replies: text: :抱歉,但我不能不说你简直是乱搞,"魔术数字"那一部分的主题明显就不是尾数,编辑前请先仔细看一遍好吗?逗号也不能滥用,这么混乱再改下去连B级条目都算不上吧.抱歉我必须撤掉你加入的内容并存档于此,以后想清楚怎么整理好再加入吧. - Dr. Cravix ♬La Pluie 2012年6月8日 (五) 14:38 (UTC)

text: 于是存档如下,mathjax里乱糟糟的"<"与">"号也已校正.顺便对某IP用户说一句:本质上是br clear=all,可能会造成潜在的隐患(最直接的就是隐藏页面bug),请不要再度滥用这一模板,好心也会办坏事的. - Dr. Cravix ♬La Pluie 2012年6月8日 (五) 15:07 (UTC)
#以右图中的数字0.15625为例，其二进制表示法为0.00101000...（即，0.15625=0.125+0.03125=1/8+1/32）。根据浮点数的存储规则，“0.00101...”需要表示成1.01*10-3。其中有效数字0.00101用科学计数法表示为\scriptstyle 1.01 \times 2^{-3}（二进制），其首位数1不需要存储在浮点数当中，于是剩下的“0100...”即成为右图中红色部分的内容。而指数-3，则需要加上127之后得到124（二进制01111100），成为右图中的绿色部分。蓝色部分为原本有理数的符号位，大于等于0用0表示，小于0用1表示。
#运算解释
这一公式的第一步是进行了整数的右位移操作（高位向低位位移），其实际效果即是将指数部分（右上图绿色部分）除以2。同时，尾数部分也因为右位移操作而除以2，而指数部分最低位也会随之成为新尾数部分的最高位。需要注意的是，由于有效数字的最高位并不在存储结构中而使其不会受到右移操作的影响，因此这第一步的整体效果近似于指数部分\scriptstyle E除以2（注意，不是实际的指数\scriptstyle E-B）。此运算步骤所做的操作很容易理解，但是实际上至此仍难以理解其真正用意。
上述代码中魔术数字0x5f3759df的二进制展开形式为0 10111110 01101110101100111011111，即相当于\scriptstyle 1.43243014812469491421875 \times 2^{63}，这一个数值是整个算法当中最难以理解的部分。其中指数部分\scriptstyle E = 10111110（二进制） = 190（十进制），是整个运算的过程的关键。为了便于理解，在此回顾原始计算公式：
Y=X^{-\frac{1}{2}}
设
X=X_m \times 2^{X_e}，其中X_m为有效数字部分且满足1 \le X_m \lt 2，X_e = E-B为指数，易得：
Y=X^{-\frac{1}{2}}= ( X_m \times 2^{X_e} ) ^{-\frac{1}{2}} = X_m^{-\frac{1}{2}} \times 2^{-\frac{X_e}{2}}  
由于\scriptstyle1 \le X_m \lt 2，因此\scriptstyle 0.707106 \lt X_m^{-\tfrac{1}{2}} \le 1，易见此部分引起的变化较为有限。使近似值快速趋向结果值的关键之一，是后一部分\scriptstyle2^{-\tfrac{X_e}{2}}  的运算，在此也即指数的运算。由于浮点数当中的指数部分保存的是\scriptstyle E而不是\scriptstyle X_e=E-B，因此不能直接用\scriptstyle E_y=-\tfrac{E_x}{2}来直接求得\scriptstyle {X_e}，其中\scriptstyle E_y为结果\scriptstyle Y在浮点数存储结构当中的指数部分\scriptstyle E（注意，不是真正的指数\scriptstyle E-B），而正确的计算方法应当是：
E_y - B = -\frac{X_e}{2} = -\frac{E_x - B}{2}，因此
E_y = B - \frac{E_x - B}{2}
= B - \frac{E_x}{2} + \frac{B}{2}
= 127 + 63 - (E_x \gt\gt 1)
= 190 - (E_x \gt\gt 1)
魔术数字0x5f3759df当中指数部分\scriptstyle E的值190正是由此求得。相对应的，算法中所使用的这段代码正实现了上述指数部分运算。
thread_title:  关于那堆证明