1
0
Fork 0
mirror of https://github.com/ruby/ruby.git synced 2022-11-09 12:17:21 -05:00
git-svn-id: svn+ssh://ci.ruby-lang.org/ruby/trunk@922 b2dd03c8-39d4-4d8f-98ff-823fe69b080e
This commit is contained in:
matz 2000-09-01 03:31:05 +00:00
parent 6001314def
commit 2252d4ce02
7 changed files with 44 additions and 31 deletions

View file

@ -612,8 +612,12 @@ class Matrix
def rank
if column_size > row_size
a = transpose.to_a
a_column_size = row_size
a_row_size = column_size
else
a = to_a
a_column_size = column_size
a_row_size = row_size
end
rank = 0
k = 0
@ -622,7 +626,7 @@ class Matrix
i = k
exists = true
begin
if (i += 1) > column_size - 1
if (i += 1) > a_column_size - 1
exists = false
break
end
@ -634,13 +638,13 @@ class Matrix
i = k
exists = true
begin
if (i += 1) > row_size - 1
if (i += 1) > a_row_size - 1
exists = false
break
end
end while a[k][i] == 0
if exists
k.upto(column_size - 1) do
k.upto(a_column_size - 1) do
|j|
a[j][k], a[j][i] = a[j][i], a[j][k]
end
@ -650,16 +654,16 @@ class Matrix
end
end
end
(k + 1).upto(row_size - 1) do
(k + 1).upto(a_row_size - 1) do
|i|
q = a[i][k] / akk
(k + 1).upto(column_size - 1) do
(k + 1).upto(a_column_size - 1) do
|j|
a[i][j] -= a[k][j] * q
end
end
rank += 1
end while (k += 1) <= column_size - 1
end while (k += 1) <= a_column_size - 1
return rank
end