Feature #10057 » add_test_for_Matrix#lup.patch
| test/matrix/test_matrix.rb | ||
|---|---|---|
|
assert_equal((v * d * v_inv).round(5), m)
|
||
|
end
|
||
|
def test_lup
|
||
|
m = Matrix[[1, 2], [3, 4]]
|
||
|
l, u, p = m.lup
|
||
|
assert(l.lower_triangular?)
|
||
|
assert(u.upper_triangular?)
|
||
|
assert(p.permutation?)
|
||
|
assert(l * u == p * m)
|
||
|
assert_equal(m.lup.solve([2, 5]), Vector[1, Rational(1,2)])
|
||
|
end
|
||
|
def test_row_vectors
|
||
|
assert_equal([Vector[1,2,3], Vector[4,5,6]], @m1.row_vectors)
|
||
|
end
|
||