diff options
| author | nobu <nobu@b2dd03c8-39d4-4d8f-98ff-823fe69b080e> | 2010-01-09 20:35:28 +0000 |
|---|---|---|
| committer | nobu <nobu@b2dd03c8-39d4-4d8f-98ff-823fe69b080e> | 2010-01-09 20:35:28 +0000 |
| commit | 2ae2d9dd4dae4089bc1c0a813cd520734ad33d28 (patch) | |
| tree | 0db059ee0330d04d475ab5d17f116a88d92171dd /lib | |
| parent | 6051cf527b70173a675a51084894be4dce56a0bd (diff) | |
* lib/matrix.rb (Matrix#inverse_from): use #quo. backported r9490.
* lib/matrix.rb (Matrix#determinant): ditto. [ruby-core:27507]
* lib/matrix.rb (Matrix#rank): ditto.
* lib/matrix.rb (Matrix::Scalar#initialize): ditto.
git-svn-id: svn+ssh://ci.ruby-lang.org/ruby/branches/ruby_1_8@26263 b2dd03c8-39d4-4d8f-98ff-823fe69b080e
Diffstat (limited to 'lib')
| -rw-r--r-- | lib/matrix.rb | 16 |
1 files changed, 8 insertions, 8 deletions
diff --git a/lib/matrix.rb b/lib/matrix.rb index 93ae3a9909..35ccb62625 100644 --- a/lib/matrix.rb +++ b/lib/matrix.rb @@ -549,7 +549,7 @@ class Matrix # # Returns the inverse of the matrix. - # Matrix[[1, 2], [2, 1]].inverse + # Matrix[[-1, -1], [0, -1]].inverse # => -1 1 # 0 -1 # @@ -585,7 +585,7 @@ class Matrix size.times do |i| next if i == k - q = a[i][k] / akk + q = a[i][k].quo(akk) a[i][k] = 0 (k + 1 ... size).each do |j| @@ -597,10 +597,10 @@ class Matrix end (k + 1 ... size).each do |j| - a[k][j] /= akk + a[k][j] = a[k][j].quo(akk) end size.times do |j| - @rows[k][j] /= akk + @rows[k][j] = @rows[k][j].quo(akk) end end self @@ -649,7 +649,7 @@ class Matrix # Returns the determinant of the matrix. If the matrix is not square, the # result is 0. # Matrix[[7,6], [3,9]].determinant - # => 63 + # => 45 # def determinant return 0 unless square? @@ -670,7 +670,7 @@ class Matrix end (k + 1 ... size).each do |i| - q = a[i][k] / akk + q = a[i][k].quo(akk) (k + 1 ... size).each do |j| a[i][j] -= a[k][j] * q end @@ -720,7 +720,7 @@ class Matrix end (k + 1 ... a_row_size).each do |i| - q = a[i][k] / akk + q = a[i][k].quo(akk) (k + 1... a_column_size).each do |j| a[i][j] -= a[k][j] * q end @@ -879,7 +879,7 @@ class Matrix self * other.inverse else x, y = other.coerce(self) - x / y + x.quo(y) end end |