1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
|
$:.push File.dirname(__FILE__) + "/../.."
require "lapack_test"
class GbsvTest < Test::Unit::TestCase
include LapackTest
def setup
@ab = Hash.new
@b = Hash.new
@b_exp = Hash.new
@ipiv_exp = Hash.new
@kl = 1
@ku = 2
@ab[:r] = NMatrix[[-0.23, 2.54, -3.66, 0.00],
[-6.98, 2.46, -2.73, -2.13],
[ 0.00, 2.56, 2.46, 4.07],
[ 0.00, 0.00, -4.78, -3.82]].to_lb(@kl, @ku, @kl)
@b[:r] = NVector[[4.42, 27.13, -6.14, 10.50]]
@b_exp[:r] = NArray[[-2.0, 3.0, 1.0, -4.0]]
@ipiv_exp[:r] = NArray[2, 3, 3, 4]
@ab[:c] = NMatrix[[-1.65+2.26*I, -2.05-0.85*I, 0.97-2.84*I, 0.00+0.00*I],
[ 0.00+6.30*I, -1.48-1.75*I, -3.99+4.01*I, 0.59-0.48*I],
[ 0.00+0.00*I, -0.77+2.83*I, -1.06+1.94*I, 3.33-1.04*I],
[ 0.00+0.00*I, 0.00+0.00*I, 4.48-1.09*I, -0.46-1.72*I]].to_lb(@kl, @ku, @kl)
@b[:c] = NVector[[-1.06+21.50*I, -22.72-53.90*I, 28.24-38.60*I, -34.56+16.73*I]]
@b_exp[:c] = NArray[[-3.0+2.0*I, 1.0-7.0*I, -5.0+4.0*I, 6.0-8.0*I]]
@ipiv_exp[:c] = NArray[2, 3, 3, 4]
end
%w(s d c z).each do |x|
method = "#{x}gbsv"
rc = LapackTest.get_rc(x)
define_method("test_#{method}") do
ipiv, info, ab, b = NumRu::Lapack.send(method, @kl, @ku, @ab[rc], @b[rc])
assert_equal 0, info
assert_narray @b_exp[rc], b
assert_equal @ipiv_exp[rc], ipiv
end
end
end
|
