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
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
|
#!/usr/bin/env ruby
require("gsl")
require("test/unit")
class VectorComplexTest < Test::Unit::TestCase
def test_vector_complex_get
v = GSL::Vector::Complex.indgen(5)
assert_equal(GSL::Vector::Complex[[3,0],[1,0],[2,0]], v.get([3, 1, 2]))
end
def test_vector_complex_addsub
a = GSL::Vector::Complex[[-2, 5], [ 4, -1]]
b = GSL::Vector::Complex[[10, 30], [20, 40]]
c = GSL::Vector::Complex[[ 8, 35], [24, 39]]
d = GSL::Vector::Complex[[12, 25], [16, 41]]
assert_equal(c, a+b)
assert_equal(d, b-a)
end
def test_vector_complex_collect
v = GSL::Vector::Complex.indgen(5)
u = GSL::Vector::Complex[[0,0], [1,0], [4,0], [9,0], [16,0]]
w = v.collect { |val| val*val }
assert_equal(u, w)
end
def test_vector_complex_subvector
v = GSL::Vector::Complex.indgen(12)
# args = []
vv = v.subvector
assert_not_equal(v.object_id, vv.object_id)
assert_equal(v.subvector, v)
# args = [Fixnum]
vv = v.subvector(3)
assert_equal([0, 0, 1, 0, 2, 0], vv.to_a)
assert_nothing_raised("subvector(-1)") {v.subvector(-1)}
vv = v.subvector(-1)
assert_equal([11, 0], vv.to_a)
vv = v.subvector(-2)
assert_equal([10, 0, 11, 0], vv.to_a)
assert_raise(RangeError) {v.subvector(-13)}
# args = [Fixnum, Fixnum]
vv = v.subvector(2, 3)
assert_equal([2, 0, 3, 0, 4, 0], vv.to_a)
vv = v.subvector(-4, 3)
assert_equal([8, 0, 9, 0, 10, 0], vv.to_a)
assert_nothing_raised("subvector(-4, -3)") {v.subvector(-4, -3)}
vv = v.subvector(-4, -3)
assert_equal([8, 0, 7, 0, 6, 0], vv.to_a)
assert_raise(GSL::ERROR::EINVAL) {v.subvector(-11, -3)}
# args = [Fixnum, Fixnum, Fixnum]
vv = v.subvector(1, 3, 4)
assert_equal([1, 0, 4, 0, 7, 0, 10, 0], vv.to_a)
# args = [Range]
tests = {
# ( range ) => [0, 0, 1, 0, 2, 0, 3, 0, 4, 0, 5, 0, 6, 0, 7, 0, 8, 0, 9, 0, 10, 0, 11, 0]
( 1.. 3) => [ 1, 0, 2, 0, 3, 0 ],
( 1... 3) => [ 1, 0, 2, 0, ],
( 3.. 1) => [ 3, 0, 2, 0, 1, 0 ],
( 3... 1) => [ 3, 0, 2, 0 ],
(-7.. 9) => [ 5, 0, 6, 0, 7, 0, 8, 0, 9, 0 ],
(-7... 9) => [ 5, 0, 6, 0, 7, 0, 8, 0 ],
( 4.. -3) => [ 4, 0, 5, 0, 6, 0, 7, 0, 8, 0, 9, 0 ],
( 4...-3) => [ 4, 0, 5, 0, 6, 0, 7, 0, 8, 0 ],
( 2.. -2) => [ 2, 0, 3, 0, 4, 0, 5, 0, 6, 0, 7, 0, 8, 0, 9, 0, 10, 0 ],
( 2...-2) => [ 2, 0, 3, 0, 4, 0, 5, 0, 6, 0, 7, 0, 8, 0, 9, 0 ],
(-2.. 2) => [ 10, 0, 9, 0, 8, 0, 7, 0, 6, 0, 5, 0, 4, 0, 3, 0, 2, 0 ],
(-2... 2) => [ 10, 0, 9, 0, 8, 0, 7, 0, 6, 0, 5, 0, 4, 0, 3, 0 ],
(-3.. -1) => [ 9, 0, 10, 0, 11, 0],
(-3...-1) => [ 9, 0, 10, 0 ],
(-1.. -3) => [ 11, 0, 10, 0, 9, 0],
(-1...-3) => [ 11, 0, 10, 0 ],
# Add more test cases here...
}
tests.each do |r, x|
assert_nothing_raised("subvector(#{r})") {v.subvector(r)}
assert_equal(x, v.subvector(r).to_a, "subvector(#{r})")
end
# args = [Range, Fixnum]
tests = {
# [( range ), s] => [0, 0, 1, 0, 2, 0, 3, 0, 4, 0, 5, 0, 6, 0, 7, 0, 8, 0, 9, 0, 10, 0, 11, 0]
[( 1.. 6), 2] => [ 1, 0, 3, 0, 5, 0 ],
[( 1... 6), 2] => [ 1, 0, 3, 0, 5, 0 ],
[( 0.. 6), 3] => [0, 0, 3, 0, 6, 0 ],
[( 0... 6), 3] => [0, 0, 3, 0 ],
# Add more test cases here...
}
tests.each do |(r,s), x|
assert_nothing_raised("subvector(#{r},#{s})") {v.subvector(r)}
assert_equal(x, v.subvector(r,s).to_a, "subvector(#{r},#{s})")
end
end
end
|
