At 23 Aug 2003 17:43:05 +0900,
R4000 2.2 wrote:
> 
> こんにちは、立花@鎌ヶ谷市です。
> 
> Fujii Hironori <fujii@chi.its.hiroshima-cu.ac.jp> writes:
> 
> > sqrt (x*x) == x が真となることを前提としたプログラムを見たのですが、
> > 真となるものですか。
> 
> 常には成り立たないと思います。

やはりそうでしたか。
どうもありがとうございます。

> > sqrt をニュートン法で誤差が十分小さくなったら終了とし、
> > 誤差がでる処理系があるのではないかと思いまして。
> > ちなみに、いくつか実験したら真でした。
> 
> 手元(gcc version 2.95.3, glibc 2.2 on Plamo Linux2.2.6)で
> Fujii さんのテストプログラムを動かしてみたら、

立花さんのテストをみて気づいたのですが、
123.456ってのは2進数では循環小数でまともに2乗できない数値でした。

---
藤井宏憲