[ruby-dev:27601] Re: Matrix class is broken without mathn

まつもと ゆきひろです

In message "Re: [ruby-dev:27596] Re: Matrix class is broken without mathn"
    on Mon, 31 Oct 2005 10:47:27 +0900, Masahiro Sakai (酒井政裕) <[e-mail]> writes:

|matrix.rb は Numeric#/ ではなく Numeric#quo を使えば良いのでは
|ないでしょうか。
|Numeric#quo は正にこういう時のためのメソッドだと思います。

確かに。こんなパッチを作ってみましたが、正しく動作しているで
しょうか。元メールの問題は解決しているようですが。私はどうに
も算数が苦手で。

--- lib/matrix.rb	27 Oct 2005 08:18:37 -0000	1.22
+++ lib/matrix.rb	31 Oct 2005 03:55:37 -0000
@@ -604,3 +604,3 @@ class Matrix
         next if i == k
-        q = a[i][k] / akk
+        q = a[i][k].quo(akk)
         a[i][k] = 0
@@ -619,3 +619,3 @@ class Matrix
         |j|
-        a[k][j] /= akk
+        a[k][j] = a[k][j].quo(akk)
       end
@@ -623,3 +623,3 @@ class Matrix
         |j|
-        @rows[k][j] /= akk
+        @rows[k][j] = @rows[k][j].quo(akk)
       end
@@ -694,3 +694,3 @@ class Matrix
         |i|
-        q = a[i][k] / akk
+        q = a[i][k].quo(akk)
         (k + 1).upto(size) do
@@ -760,3 +760,3 @@ class Matrix
         |i|
-        q = a[i][k] / akk
+        q = a[i][k].quo(akk)
         (k + 1).upto(a_column_size - 1) do
@@ -927,3 +927,3 @@ class Matrix
         x, y = other.coerce(self)
-        x / y
+        x.quo(y)
       end