[ruby-dev:5084] [PATCH] (Re: [PATCH] configure w/ isinf())

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

ruby-dev:5007 で扱った無限／非数の整数への変換
等の問題へのパッチの更新です。

これにより ruby-dev:5007 は obsolate になりました。

展開したての ruby-1.3.1-990210 に適用して下さい。

 * configure.in で isinf(), isnan() 及び finite() の有無を調べます。
   なかった場合の為 missing/{isinf,isnan,finite}.c を用意しました。
   これにより ruby インタプリタ及び拡張モジュールの中で、
   isinf(), isnan() 及び finite() が存在すると仮定したコーディングが
   可能になります。 
   configure.in missing/isinf.c missing/isnan.c missing/finite.c

 + ある種の最適化が isnan() の正しい生成を阻害する問題を回避しました。
   n != n から常に偽を生成する問題（n が NaN の場合に限り真が正しい）

 * 例外 DomainError の追加を行ないました。
   Float -> Integer への変換で以下の条件で無限ループに陥っていました。

     o 正の無限(Inf)を整数に変換しようとした時
     o 負の無限(-Inf)を整数に変換しようとした時
     o 非数(NaN)を整数に変換しようとした時

   このパッチで DomainError 例外を raise する様にしました。
   bignam.c numeric.c ruby.h

 * FreeBSD/i386 プラットホームでの、浮動小数点桁溢れでの
   SIGFPE の発生を抑止しました。
   numeric.c

このパッチの結果、無限ループに陥らずに、

     %  ruby -e '1e999.to_i'
     /var/tmp//rbY26973:1: warning: Float 1e999 out of range
     /var/tmp//rbY26973:1:in `to_i': Inifinity (DomainError)
	     from /var/tmp//rbY26973:1
     % ruby -e '-1e999.to_i'
     /var/tmp//rbV26974:1: warning: Float -1e999 out of range
     /var/tmp//rbV26974:1:in `to_i': -Inifinity (DomainError)
	     from /var/tmp//rbV26974:1
     % ruby -e '(0.0%0.0).to_i'
     /var/tmp//rbw26975:1:in `to_i': NaN (DomainError)
	     from /var/tmp//rbw26975:1
     % 

の様な、実行結果になります。

	えぐち

diff --exclude=configure -rudN ../ruby-1.3.1-990210/bignum.c ./bignum.c
--- ../ruby-1.3.1-990210/bignum.c	Wed Feb 10 17:44:23 1999
+++ ./bignum.c	Sun Feb 14 12:47:48 1999
@@ -403,6 +403,13 @@
     VALUE z;
     double u = (d < 0)?-d:d;
 
+    if (isinf(d)) {
+	rb_raise(rb_eDomainError, d < 0 ? "-Inifinity" : "Inifinity");
+    }
+    if (isnan(d)) {
+	rb_raise(rb_eDomainError, "NaN");
+    }
+
     while (!POSFIXABLE(u) || 0 != (long)u) {
 	u /= (double)(BIGRAD);
 	i++;
diff --exclude=configure -rudN ../ruby-1.3.1-990210/configure.in ./configure.in
--- ../ruby-1.3.1-990210/configure.in	Wed Feb 10 14:40:12 1999
+++ ./configure.in	Thu Feb 11 13:04:28 1999
@@ -187,7 +187,8 @@
 AC_FUNC_VFORK
 AC_FUNC_MEMCMP
 AC_REPLACE_FUNCS(dup2 memmove mkdir strcasecmp strerror strftime\
-		 strchr strstr strtoul strdup crypt flock vsnprintf)
+		 strchr strstr strtoul strdup crypt flock vsnprintf\
+		 isinf isnan finite)
 AC_CHECK_FUNCS(fmod killpg drand48 random wait4 waitpid syscall getcwd\
 	      truncate chsize times utimes fcntl lockf setitimer\
 	      setruid seteuid setreuid setrgid setegid setregid\
diff --exclude=configure -rudN ../ruby-1.3.1-990210/missing/finite.c ./missing/finite.c
--- ../ruby-1.3.1-990210/missing/finite.c	Thu Jan  1 09:00:00 1970
+++ ./missing/finite.c	Thu Feb 11 13:03:06 1999
@@ -0,0 +1,6 @@
+int
+finite(n)
+    double n;
+{
+    return !isnan(n) && !isinf(n);
+}
diff --exclude=configure -rudN ../ruby-1.3.1-990210/missing/isinf.c ./missing/isinf.c
--- ../ruby-1.3.1-990210/missing/isinf.c	Thu Jan  1 09:00:00 1970
+++ ./missing/isinf.c	Thu Feb 11 13:03:11 1999
@@ -0,0 +1,18 @@
+static double zero()	{ return 0.0; }
+static double one()	{ return 1.0; }
+static double inf()	{ return one() / zero(); }
+
+int
+isinf(n)
+    double n;
+{
+    static double pinf = 0.0;
+    static double ninf = 0.0;
+
+    if (pinf == 0.0) {
+	pinf = inf();
+	ninf = -pinf;
+    }
+    return memcmp(&n, &pinf, sizeof n) == 0
+	|| memcmp(&n, &ninf, sizeof n) == 0;
+}
diff --exclude=configure -rudN ../ruby-1.3.1-990210/missing/isnan.c ./missing/isnan.c
--- ../ruby-1.3.1-990210/missing/isnan.c	Thu Jan  1 09:00:00 1970
+++ ./missing/isnan.c	Sat Feb 13 09:30:11 1999
@@ -0,0 +1,18 @@
+static int double_ne();
+
+int
+isnan(n)
+  double n;
+{
+  return double_ne(n, n);
+}
+
+static
+int
+double_ne(n1, n2)
+  double n1, n2;
+{
+  return n1 != n2;
+}
+
+
diff --exclude=configure -rudN ../ruby-1.3.1-990210/numeric.c ./numeric.c
--- ../ruby-1.3.1-990210/numeric.c	Wed Feb 10 17:44:25 1999
+++ ./numeric.c	Sat Feb 13 10:33:53 1999
@@ -26,6 +26,7 @@
 VALUE rb_cFixnum;
 
 VALUE rb_eZeroDiv;
+VALUE rb_eDomainError;
 
 ID rb_frame_last_func();
 VALUE rb_float_new();
@@ -37,6 +38,12 @@
     rb_raise(rb_eZeroDiv, "divided by 0");
 }
 
+void
+rb_num_domain()
+{
+    rb_raise(rb_eDomainError, "Domain");
+}
+
 static VALUE
 num_coerce(x, y)
     VALUE x, y;
@@ -1405,12 +1412,13 @@
 {
 #ifdef __FreeBSD__
     /* allow divide by zero -- Inf */
-    fpsetmask(fpgetmask() & ~(FP_X_DZ|FP_X_INV));
+    fpsetmask(fpgetmask() & ~(FP_X_DZ|FP_X_INV|FP_X_OFL));
 #endif
     coerce = rb_intern("coerce");
     to_i = rb_intern("to_i");
 
     rb_eZeroDiv = rb_define_class("ZeroDivisionError", rb_eStandardError);
+    rb_eDomainError = rb_define_class("DomainError", rb_eStandardError);
     rb_cNumeric = rb_define_class("Numeric", rb_cObject);
 
     rb_include_module(rb_cNumeric, rb_mComparable);
diff --exclude=configure -rudN ../ruby-1.3.1-990210/ruby.h ./ruby.h
--- ../ruby-1.3.1-990210/ruby.h	Tue Feb  9 14:13:15 1999
+++ ./ruby.h	Sat Feb 13 10:09:11 1999
@@ -494,6 +494,7 @@
 EXTERN VALUE rb_eTypeError;
 EXTERN VALUE rb_eZeroDiv;
 EXTERN VALUE rb_eNotImpError;
+EXTERN VALUE rb_eDomainError;
 
 #if defined(__GNUC__) && __GNUC__ >= 2 && !defined(RUBY_NO_INLINE)
 extern __inline__ VALUE rb_class_of _((VALUE));