▼ げんたさん
> > どう頑張っても単精度変数に一度落としてからでないと比較できないなんてアルゴリズム上の問題じゃないよなあ。
> そういうものだと思うしか無いと思います。今回のように無理矢理精度を落としたいという方がまれだと思いますが。
(途中省略)
> > 困ったモンです。>MS それに引き換えgccはフリーなのにさすがだな。
> うーん、どうだか。1.0+a==1.0というのは明らかに間違いですからaの精度以上に正しい結果(1.0+a>1.0)を出す方が親切と考える人がいてもおかしくない。それに、いちいち精度を変換して比較するのはそれだけ遅くなっているとも考えられますし。

そういうことではないんです。主旨は「計算には精度というものがあり、その範囲でしか計算結果を信用してはいけない」という教育的内容なわけです。だから内部的には1.0+a=1.0ということも起こり得るという意味で誤りではありません。精度よく計算したい場合には、あるいは予期しない結果になっていた場合には精度の影響が考えられるので高精度で計算しなさいという主旨。
Fortranでは倍精度の他に四倍精度もほぼ標準的にあります。(JIS規格かまでは知りませんが)
一方Cは長精度は処理系によってはありますが、長さまで規定されてないし、1.5倍とかいい加減です。
言語の標準精度はだいたい単精度で、Cのみ倍精度とちょっと変わっている、と思っているのですが、最近のPascal(Delphi)やBasic(VB)、javaなどはどうなんでしょ？
別に倍精度と長精度の比較でもいいんですが、標準的ではないのと処理系によって異なるというのを嫌ったわけです。結果的にはそれも(プログラムの書き方により)処理系によってしまったわけですが。長精度でやっても今度は長精度で比較してほしいところを倍精度で比較されてしまうわけだから(ですよね？)やっぱり良くない仕様だと思います。
単精度は電卓にも劣る精度だとかわざわざ精度を落として計算するのはあほらしい、とか今では単精度の方がむしろ時間がかかる、というのは精度ということを知ることとは別の次元の、むしろもっと泥臭い話です。コンパイラの作りによるわけだから。
精度という考えはどの計算用言語にもある話です(よね？)。結局、C言語では単精度を使うといろいろ予期しないが起こり得るし、使うメリットもほとんどないので倍精度で計算しなさい、というあまり明解でない回答になってしまいます。