[ruby-dev:4805] Re: compare to NaN (Re: Re: call for finite/isnan testers)

えぐち＠エスアンドイー です。

>>> On Sun, 7 Feb 1999 11:04:31 +0900, [e-mail] (GOTO Kentaro) said:

gotoken> ごとけんです
gotoken> 
gotoken> In message "[ruby-dev:4793] Re: compare to NaN (Re:  Re: call for finite/isnan testers)"
gotoken>     on 99/02/07, EGUCHI Osamu <[e-mail]> writes:
gotoken> 
gotoken> >えぐち＠エスアンドイー です。
gotoken> (中略)
gotoken> >gotoken> もう一つの考え方は、NaN, +-Inf それぞれ一つで、
gotoken> >gotoken> 一括して Numeric の定数だとする考え方です。
gotoken> >
gotoken> >これに賛成です。
gotoken>
gotoken> ややこしいので、これがよさそうですね。

ところで
  Numeric::INFINITY.type
  => Float
で良いでしょうか？

gotoken> >gotoken> >  isinf(),isnan()の移植性をどうするか．ごとけんさんの調査待
gotoken> >gotoken> >  ちか．
gotoken> >gotoken> 
gotoken> >gotoken> isnan() 相当はどこにでもあって、finite() もしくは isinf() 
gotoken> >gotoken> ならどこにでもある模様です。 
gotoken> >gotoken> BeOS, NeXT, それから 64 bit の Alpha とかはよく分かりませんが。
gotoken> >
gotoken> >ビット単位の比較でよいので自家製と言う選択もありですが、、
gotoken> >まずは、その数の生成がポータブルに出来るかですね。
gotoken> 
gotoken> そうなんですよね。ビットの並びが決まっていて、しかも
gotoken> IEEE 754 準拠なら、
gotoken> 
gotoken> #if sun || decosf1 || ...
gotoken> #define be_NaN(X)  (((union { double d; struct { unsigned :1, e:11; } s; } \
gotoken>                         *)&X)->s.e == 0x7ff)
gotoken> #elif bsd
gotoken> #define be_NaN(X) ...
gotoken> #endif
gotoken> 
gotoken> とかできるんですが。

『ビット単位の比較』と言う言葉が非常に曖昧でした。
私なりのニュアンスの isnan() の実装を３つ考えてみました。

  triton2.eguchi:72% cat nan.c
  #include <stdio.h>
  #include <stdlib.h>
  #include <floatingpoint.h>
  extern double double_zero;
  double make_nan()
  {
      return double_zero / double_zero;
  }
  double double_zero = 0.0;

  int is_nan1(double n)
  {
      double nan = make_nan();
      return memcmp(&n, &nan, sizeof n) == 0;
  }

  int is_nan2(double n)
  {
      double nan = make_nan();
      return *(long long*)&n == *(long long*)&nan;;
  }

  int is_nan3(double n)
  {
      return n != n;
  }

  main()
  {
      double nan;

      fpsetmask(fpgetmask() & ~(FP_X_DZ|FP_X_INV));
      nan = make_nan();
      printf("NaN = %qX\n", *(long long*)&nan);
      printf("is_nan1(0.0 / 0.0)      => %d\n", is_nan1(0.0 / 0.0));
      printf("is_nan2(0.0 / 0.0)      => %d\n", is_nan2(0.0 / 0.0));
      printf("is_nan3(0.0 / 0.0)      => %d\n", is_nan3(0.0 / 0.0));
      printf("0.0 / 0.0 == 0.0 / 0.0  => %d\n", 0.0 / 0.0 == 0.0 / 0.0);
  }
  triton2.eguchi:73% ./nan
  NaN = FFF8000000000000
  is_nan1(0.0 / 0.0)      => 1
  is_nan2(0.0 / 0.0)      => 1
  is_nan3(0.0 / 0.0)      => 1
  0.0 / 0.0 == 0.0 / 0.0  => 0
  triton2.eguchi:74% 

移植性から言うと、 double nan を必要としない第３案でしょうか？

	えぐち