rubysl-matrix 1.0.0
Sign up to get free protection for your applications and to get access to all the features.
- checksums.yaml +7 -0
- data/.gitignore +17 -0
- data/.travis.yml +8 -0
- data/Gemfile +4 -0
- data/LICENSE +25 -0
- data/README.md +29 -0
- data/Rakefile +1 -0
- data/lib/matrix.rb +1 -0
- data/lib/rubysl/matrix.rb +2 -0
- data/lib/rubysl/matrix/matrix.rb +1537 -0
- data/lib/rubysl/matrix/version.rb +5 -0
- data/rubysl-matrix.gemspec +27 -0
- data/spec/I_spec.rb +5 -0
- data/spec/build_spec.rb +75 -0
- data/spec/clone_spec.rb +26 -0
- data/spec/coerce_spec.rb +11 -0
- data/spec/collect_spec.rb +5 -0
- data/spec/column_size_spec.rb +14 -0
- data/spec/column_spec.rb +36 -0
- data/spec/column_vector_spec.rb +36 -0
- data/spec/column_vectors_spec.rb +25 -0
- data/spec/columns_spec.rb +45 -0
- data/spec/conj_spec.rb +7 -0
- data/spec/conjugate_spec.rb +7 -0
- data/spec/constructor_spec.rb +68 -0
- data/spec/det_spec.rb +6 -0
- data/spec/determinant_spec.rb +6 -0
- data/spec/diagonal_spec.rb +73 -0
- data/spec/divide_spec.rb +60 -0
- data/spec/each_spec.rb +77 -0
- data/spec/each_with_index_spec.rb +84 -0
- data/spec/eigenvalue_decomposition/eigenvalue_matrix_spec.rb +10 -0
- data/spec/eigenvalue_decomposition/eigenvalues_spec.rb +23 -0
- data/spec/eigenvalue_decomposition/eigenvector_matrix_spec.rb +23 -0
- data/spec/eigenvalue_decomposition/eigenvectors_spec.rb +25 -0
- data/spec/eigenvalue_decomposition/initialize_spec.rb +27 -0
- data/spec/eigenvalue_decomposition/to_a_spec.rb +19 -0
- data/spec/element_reference_spec.rb +24 -0
- data/spec/empty_spec.rb +69 -0
- data/spec/eql_spec.rb +12 -0
- data/spec/equal_value_spec.rb +10 -0
- data/spec/exponent_spec.rb +66 -0
- data/spec/find_index_spec.rb +147 -0
- data/spec/fixtures/classes.rb +7 -0
- data/spec/hash_spec.rb +14 -0
- data/spec/hermitian_spec.rb +37 -0
- data/spec/identity_spec.rb +5 -0
- data/spec/imag_spec.rb +7 -0
- data/spec/imaginary_spec.rb +7 -0
- data/spec/inspect_spec.rb +30 -0
- data/spec/inv_spec.rb +6 -0
- data/spec/inverse_from_spec.rb +5 -0
- data/spec/inverse_spec.rb +6 -0
- data/spec/lower_triangular_spec.rb +25 -0
- data/spec/lup_decomposition/determinant_spec.rb +24 -0
- data/spec/lup_decomposition/initialize_spec.rb +14 -0
- data/spec/lup_decomposition/l_spec.rb +19 -0
- data/spec/lup_decomposition/p_spec.rb +19 -0
- data/spec/lup_decomposition/solve_spec.rb +54 -0
- data/spec/lup_decomposition/to_a_spec.rb +36 -0
- data/spec/lup_decomposition/u_spec.rb +19 -0
- data/spec/map_spec.rb +5 -0
- data/spec/minor_spec.rb +90 -0
- data/spec/minus_spec.rb +45 -0
- data/spec/multiply_spec.rb +73 -0
- data/spec/new_spec.rb +7 -0
- data/spec/normal_spec.rb +27 -0
- data/spec/orthogonal_spec.rb +27 -0
- data/spec/permutation_spec.rb +33 -0
- data/spec/plus_spec.rb +45 -0
- data/spec/rank_spec.rb +24 -0
- data/spec/real_spec.rb +45 -0
- data/spec/rect_spec.rb +7 -0
- data/spec/rectangular_spec.rb +7 -0
- data/spec/regular_spec.rb +34 -0
- data/spec/round_spec.rb +24 -0
- data/spec/row_size_spec.rb +12 -0
- data/spec/row_spec.rb +37 -0
- data/spec/row_vector_spec.rb +33 -0
- data/spec/row_vectors_spec.rb +25 -0
- data/spec/rows_spec.rb +40 -0
- data/spec/scalar/Fail_spec.rb +5 -0
- data/spec/scalar/Raise_spec.rb +5 -0
- data/spec/scalar/divide_spec.rb +5 -0
- data/spec/scalar/exponent_spec.rb +5 -0
- data/spec/scalar/included_spec.rb +5 -0
- data/spec/scalar/initialize_spec.rb +5 -0
- data/spec/scalar/minus_spec.rb +5 -0
- data/spec/scalar/multiply_spec.rb +5 -0
- data/spec/scalar/plus_spec.rb +5 -0
- data/spec/scalar_spec.rb +66 -0
- data/spec/shared/collect.rb +29 -0
- data/spec/shared/conjugate.rb +21 -0
- data/spec/shared/determinant.rb +53 -0
- data/spec/shared/equal_value.rb +35 -0
- data/spec/shared/identity.rb +21 -0
- data/spec/shared/imaginary.rb +22 -0
- data/spec/shared/inverse.rb +42 -0
- data/spec/shared/rectangular.rb +20 -0
- data/spec/shared/trace.rb +14 -0
- data/spec/shared/transpose.rb +23 -0
- data/spec/singular_spec.rb +34 -0
- data/spec/spec_helper.rb +35 -0
- data/spec/square_spec.rb +29 -0
- data/spec/symmetric_spec.rb +30 -0
- data/spec/t_spec.rb +5 -0
- data/spec/to_a_spec.rb +5 -0
- data/spec/to_s_spec.rb +5 -0
- data/spec/tr_spec.rb +6 -0
- data/spec/trace_spec.rb +6 -0
- data/spec/transpose_spec.rb +5 -0
- data/spec/unit_spec.rb +5 -0
- data/spec/unitary_spec.rb +29 -0
- data/spec/upper_triangular_spec.rb +26 -0
- data/spec/vector/cross_product_spec.rb +25 -0
- data/spec/vector/each2_spec.rb +52 -0
- data/spec/vector/eql_spec.rb +17 -0
- data/spec/vector/inner_product_spec.rb +25 -0
- data/spec/vector/normalize_spec.rb +19 -0
- data/spec/zero_spec.rb +55 -0
- metadata +341 -0
data/spec/plus_spec.rb
ADDED
@@ -0,0 +1,45 @@
|
|
1
|
+
require File.expand_path('../fixtures/classes', __FILE__)
|
2
|
+
require 'matrix'
|
3
|
+
|
4
|
+
describe "Matrix#+" do
|
5
|
+
before :each do
|
6
|
+
@a = Matrix[ [1,2], [3,4] ]
|
7
|
+
@b = Matrix[ [4,5], [6,7] ]
|
8
|
+
end
|
9
|
+
|
10
|
+
it "returns the result of adding the corresponding elements of self and other" do
|
11
|
+
(@a + @b).should == Matrix[ [5,7], [9,11] ]
|
12
|
+
end
|
13
|
+
|
14
|
+
it "returns an instance of Matrix" do
|
15
|
+
(@a + @b).should be_kind_of(Matrix)
|
16
|
+
end
|
17
|
+
|
18
|
+
it "raises a Matrix::ErrDimensionMismatch if the matrices are different sizes" do
|
19
|
+
lambda { @a + Matrix[ [1] ] }.should raise_error(Matrix::ErrDimensionMismatch)
|
20
|
+
end
|
21
|
+
|
22
|
+
it "raises a ExceptionForMatrix::ErrOperationNotDefined if other is a Numeric Type" do
|
23
|
+
lambda { @a + 2 }.should raise_error(ExceptionForMatrix::ErrOperationNotDefined)
|
24
|
+
lambda { @a + 1.2 }.should raise_error(ExceptionForMatrix::ErrOperationNotDefined)
|
25
|
+
lambda { @a + bignum_value }.should raise_error(ExceptionForMatrix::ErrOperationNotDefined)
|
26
|
+
end
|
27
|
+
|
28
|
+
ruby_bug "redmine:2365", "1.8.7" do
|
29
|
+
it "raises a TypeError if other is of wrong type" do
|
30
|
+
lambda { @a + nil }.should raise_error(TypeError)
|
31
|
+
lambda { @a + "a" }.should raise_error(TypeError)
|
32
|
+
lambda { @a + [ [1, 2] ] }.should raise_error(TypeError)
|
33
|
+
lambda { @a + Object.new }.should raise_error(TypeError)
|
34
|
+
end
|
35
|
+
end
|
36
|
+
|
37
|
+
ruby_bug "redmine #5307", "1.9.3" do
|
38
|
+
describe "for a subclass of Matrix" do
|
39
|
+
it "returns an instance of that subclass" do
|
40
|
+
m = MatrixSub.ins
|
41
|
+
(m+m).should be_an_instance_of(MatrixSub)
|
42
|
+
end
|
43
|
+
end
|
44
|
+
end
|
45
|
+
end
|
data/spec/rank_spec.rb
ADDED
@@ -0,0 +1,24 @@
|
|
1
|
+
require 'matrix'
|
2
|
+
|
3
|
+
describe "Matrix#rank" do
|
4
|
+
it "returns the rank of the Matrix" do
|
5
|
+
Matrix[ [7,6], [3,9] ].rank.should == 2
|
6
|
+
end
|
7
|
+
|
8
|
+
# The examples below seems to indicate a regression in MRI 1.9.1, as it works
|
9
|
+
# in 1.8.7. On 1.9.1 it just hangs.
|
10
|
+
ruby_bug "#1020", "1.9.1.500" do
|
11
|
+
it "doesn't loop forever" do
|
12
|
+
Matrix[ [1,2,3], [4,5,6], [7,8,9] ].rank.should == 2
|
13
|
+
Matrix[ [1, 2, 0, 3], [1, -2, 3, 0], [0, 0, 4, 8], [2, 4, 0, 6] ].rank.
|
14
|
+
should == 3
|
15
|
+
end
|
16
|
+
end
|
17
|
+
|
18
|
+
ruby_bug "revision 24969", "1.8.7" do
|
19
|
+
it "works for some easy rectangular matrices" do
|
20
|
+
Matrix[[0,0],[0,0],[1,0]].rank.should == 1
|
21
|
+
Matrix[[0,1],[0,0],[1,0]].rank.should == 2
|
22
|
+
end
|
23
|
+
end
|
24
|
+
end
|
data/spec/real_spec.rb
ADDED
@@ -0,0 +1,45 @@
|
|
1
|
+
require File.expand_path('../fixtures/classes', __FILE__)
|
2
|
+
require 'matrix'
|
3
|
+
|
4
|
+
ruby_version_is "1.9" do
|
5
|
+
describe "Matrix#real?" do
|
6
|
+
it "returns true for matrices with all real entries" do
|
7
|
+
Matrix[ [1, 2], [3, 4] ].real?.should be_true
|
8
|
+
Matrix[ [1.9, 2], [3, 4] ].real?.should be_true
|
9
|
+
end
|
10
|
+
|
11
|
+
it "returns true for empty matrices" do
|
12
|
+
Matrix.empty.real?.should be_true
|
13
|
+
end
|
14
|
+
|
15
|
+
it "returns false if one element is a Complex" do
|
16
|
+
Matrix[ [Complex(1,1), 2], [3, 4] ].real?.should be_false
|
17
|
+
end
|
18
|
+
|
19
|
+
conflicts_with :CMath do
|
20
|
+
it "returns false if one element is a Complex whose imaginary part is 0" do
|
21
|
+
Matrix[ [Complex(1,0), 2], [3, 4] ].real?.should be_false
|
22
|
+
end
|
23
|
+
end
|
24
|
+
end
|
25
|
+
|
26
|
+
describe "Matrix#real" do
|
27
|
+
it "returns a matrix with the real part of the elements of the receiver" do
|
28
|
+
Matrix[ [1, 2], [3, 4] ].real.should == Matrix[ [1, 2], [3, 4] ]
|
29
|
+
Matrix[ [1.9, Complex(1,1)], [Complex(-0.42, 0), 4] ].real.should == Matrix[ [1.9, 1], [-0.42, 4] ]
|
30
|
+
end
|
31
|
+
|
32
|
+
it "returns empty matrices on the same size if empty" do
|
33
|
+
Matrix.empty(0, 3).real.should == Matrix.empty(0, 3)
|
34
|
+
Matrix.empty(3, 0).real.should == Matrix.empty(3, 0)
|
35
|
+
end
|
36
|
+
|
37
|
+
ruby_bug "redmine #5307", "1.9.3" do
|
38
|
+
describe "for a subclass of Matrix" do
|
39
|
+
it "returns an instance of that subclass" do
|
40
|
+
MatrixSub.ins.real.should be_an_instance_of(MatrixSub)
|
41
|
+
end
|
42
|
+
end
|
43
|
+
end
|
44
|
+
end
|
45
|
+
end
|
data/spec/rect_spec.rb
ADDED
@@ -0,0 +1,34 @@
|
|
1
|
+
require 'matrix'
|
2
|
+
|
3
|
+
describe "Matrix#regular?" do
|
4
|
+
|
5
|
+
ruby_bug "#1020", "1.8.7" do
|
6
|
+
it "returns false for singular matrices" do
|
7
|
+
m = Matrix[ [1,2,3], [3,4,3], [0,0,0] ]
|
8
|
+
m.regular?.should be_false
|
9
|
+
|
10
|
+
m = Matrix[ [1,2,9], [3,4,9], [1,2,9] ]
|
11
|
+
m.regular?.should be_false
|
12
|
+
end
|
13
|
+
|
14
|
+
it "returns true if the Matrix is regular" do
|
15
|
+
Matrix[ [0,1], [1,0] ].regular?.should be_true
|
16
|
+
end
|
17
|
+
end
|
18
|
+
|
19
|
+
ruby_bug "", "1.8.7" do
|
20
|
+
it "returns true for an empty 0x0 matrix" do
|
21
|
+
Matrix.empty(0,0).regular?.should be_true
|
22
|
+
end
|
23
|
+
|
24
|
+
it "raises an error for rectangular matrices" do
|
25
|
+
lambda {
|
26
|
+
Matrix[[1], [2], [3]].regular?
|
27
|
+
}.should raise_error(Matrix::ErrDimensionMismatch)
|
28
|
+
|
29
|
+
lambda {
|
30
|
+
Matrix.empty(3,0).regular?
|
31
|
+
}.should raise_error(Matrix::ErrDimensionMismatch)
|
32
|
+
end
|
33
|
+
end
|
34
|
+
end
|
data/spec/round_spec.rb
ADDED
@@ -0,0 +1,24 @@
|
|
1
|
+
require File.expand_path('../fixtures/classes', __FILE__)
|
2
|
+
require 'matrix'
|
3
|
+
|
4
|
+
ruby_version_is "1.9.3" do
|
5
|
+
describe "Matrix#round" do
|
6
|
+
it "returns a matrix with all entries rounded" do
|
7
|
+
Matrix[ [1, 2.34], [5.67, 8] ].round.should == Matrix[ [1, 2], [6, 8] ]
|
8
|
+
Matrix[ [1, 2.34], [5.67, 8] ].round(1).should == Matrix[ [1, 2.3], [5.7, 8] ]
|
9
|
+
end
|
10
|
+
|
11
|
+
it "returns empty matrices on the same size if empty" do
|
12
|
+
Matrix.empty(0, 3).round.should == Matrix.empty(0, 3)
|
13
|
+
Matrix.empty(3, 0).round(42).should == Matrix.empty(3, 0)
|
14
|
+
end
|
15
|
+
|
16
|
+
ruby_bug "redmine #5307", "1.9.3" do
|
17
|
+
describe "for a subclass of Matrix" do
|
18
|
+
it "returns an instance of that subclass" do
|
19
|
+
MatrixSub.ins.round.should be_an_instance_of(MatrixSub)
|
20
|
+
end
|
21
|
+
end
|
22
|
+
end
|
23
|
+
end
|
24
|
+
end
|
@@ -0,0 +1,12 @@
|
|
1
|
+
require 'matrix'
|
2
|
+
|
3
|
+
describe "Matrix#row_size" do
|
4
|
+
it "returns the number rows" do
|
5
|
+
Matrix[ [1,2], [3, 4], [5, 6] ].row_size.should == 3
|
6
|
+
end
|
7
|
+
|
8
|
+
it "returns the number rows even for some empty matrices" do
|
9
|
+
Matrix[ [], [], [] ].row_size.should == 3
|
10
|
+
end
|
11
|
+
|
12
|
+
end
|
data/spec/row_spec.rb
ADDED
@@ -0,0 +1,37 @@
|
|
1
|
+
require 'matrix'
|
2
|
+
|
3
|
+
describe "Matrix#row" do
|
4
|
+
before :all do
|
5
|
+
@m = Matrix[ [1, 2], [2, 3], [3, 4] ]
|
6
|
+
end
|
7
|
+
|
8
|
+
it "returns a Vector when called without a block" do
|
9
|
+
@m.row(0).should == Vector[1,2]
|
10
|
+
end
|
11
|
+
|
12
|
+
it "yields the elements of the row when called with a block" do
|
13
|
+
a = []
|
14
|
+
@m.row(0) {|x| a << x}
|
15
|
+
a.should == [1,2]
|
16
|
+
end
|
17
|
+
|
18
|
+
it "counts backwards for negative argument" do
|
19
|
+
@m.row(-1).should == Vector[3, 4]
|
20
|
+
end
|
21
|
+
|
22
|
+
ruby_bug "redmine:1532", "1.8.7" do
|
23
|
+
it "returns self when called with a block" do
|
24
|
+
@m.row(0) { |x| x }.should equal(@m)
|
25
|
+
end
|
26
|
+
|
27
|
+
it "returns nil when out of bounds" do
|
28
|
+
@m.row(3).should == nil
|
29
|
+
@m.row(-4).should == nil
|
30
|
+
end
|
31
|
+
|
32
|
+
it "never yields when out of bounds" do
|
33
|
+
lambda { @m.row(3){ raise } }.should_not raise_error
|
34
|
+
lambda { @m.row(-4){ raise } }.should_not raise_error
|
35
|
+
end
|
36
|
+
end
|
37
|
+
end
|
@@ -0,0 +1,33 @@
|
|
1
|
+
require File.expand_path('../fixtures/classes', __FILE__)
|
2
|
+
require 'matrix'
|
3
|
+
|
4
|
+
describe "Matrix.row_vector" do
|
5
|
+
|
6
|
+
it "returns a Matrix" do
|
7
|
+
Matrix.row_vector([]).should be_an_instance_of(Matrix)
|
8
|
+
end
|
9
|
+
|
10
|
+
it "returns a single-row Matrix with the specified values" do
|
11
|
+
Matrix.row_vector([1,2]).should == Matrix[ [1,2] ]
|
12
|
+
end
|
13
|
+
|
14
|
+
it "returns a 1x0 matrix when called with an empty Array" do
|
15
|
+
Matrix.row_vector([]).should == Matrix[ [] ]
|
16
|
+
end
|
17
|
+
|
18
|
+
ruby_version_is ""..."1.9" do
|
19
|
+
describe "for a subclass of Matrix" do
|
20
|
+
it "returns an instance of Matrix" do
|
21
|
+
MatrixSub.row_vector([1]).should be_an_instance_of(Matrix)
|
22
|
+
end
|
23
|
+
end
|
24
|
+
end
|
25
|
+
|
26
|
+
ruby_version_is "1.9" do
|
27
|
+
describe "for a subclass of Matrix" do
|
28
|
+
it "returns an instance of that subclass" do
|
29
|
+
MatrixSub.row_vector([1]).should be_an_instance_of(MatrixSub)
|
30
|
+
end
|
31
|
+
end
|
32
|
+
end
|
33
|
+
end
|
@@ -0,0 +1,25 @@
|
|
1
|
+
require 'matrix'
|
2
|
+
|
3
|
+
describe "Matrix#row_vectors" do
|
4
|
+
|
5
|
+
before(:each) do
|
6
|
+
@vectors = Matrix[ [1,2], [3,4] ].row_vectors
|
7
|
+
end
|
8
|
+
|
9
|
+
it "returns an Array" do
|
10
|
+
Matrix[ [1,2], [3,4] ].row_vectors.should be_an_instance_of(Array)
|
11
|
+
end
|
12
|
+
|
13
|
+
it "returns an Array of Vectors" do
|
14
|
+
@vectors.all? {|v| v.should be_an_instance_of(Vector)}
|
15
|
+
end
|
16
|
+
|
17
|
+
it "returns each row as a Vector" do
|
18
|
+
@vectors.should == [Vector[1,2], Vector[3,4]]
|
19
|
+
end
|
20
|
+
|
21
|
+
it "returns an empty Array for empty matrices" do
|
22
|
+
Matrix[].row_vectors.should == []
|
23
|
+
Matrix[ [] ].row_vectors.should == [ Vector[] ]
|
24
|
+
end
|
25
|
+
end
|
data/spec/rows_spec.rb
ADDED
@@ -0,0 +1,40 @@
|
|
1
|
+
require File.expand_path('../fixtures/classes', __FILE__)
|
2
|
+
require 'matrix'
|
3
|
+
|
4
|
+
describe "Matrix.rows" do
|
5
|
+
before :each do
|
6
|
+
@a = [1, 2]
|
7
|
+
@b = [3, 4]
|
8
|
+
@m = Matrix.rows([@a, @b])
|
9
|
+
end
|
10
|
+
|
11
|
+
it "returns a Matrix" do
|
12
|
+
@m.should be_kind_of(Matrix)
|
13
|
+
end
|
14
|
+
|
15
|
+
it "creates a matrix from argument rows" do
|
16
|
+
@m.row(0).to_a.should == @a
|
17
|
+
@m.row(1).to_a.should == @b
|
18
|
+
end
|
19
|
+
|
20
|
+
it "copies the original rows by default" do
|
21
|
+
@a << 3
|
22
|
+
@b << 6
|
23
|
+
@m.row(0).should_not equal(@a)
|
24
|
+
@m.row(1).should_not equal(@b)
|
25
|
+
end
|
26
|
+
|
27
|
+
it "references the original rows if copy is false" do
|
28
|
+
@m_ref = Matrix.rows([@a, @b], false)
|
29
|
+
@a << 3
|
30
|
+
@b << 6
|
31
|
+
@m_ref.row(0).to_a.should == @a
|
32
|
+
@m_ref.row(1).to_a.should == @b
|
33
|
+
end
|
34
|
+
|
35
|
+
describe "for a subclass of Matrix" do
|
36
|
+
it "returns an instance of that subclass" do
|
37
|
+
MatrixSub.rows([[0, 1], [0, 1]]).should be_an_instance_of(MatrixSub)
|
38
|
+
end
|
39
|
+
end
|
40
|
+
end
|
data/spec/scalar_spec.rb
ADDED
@@ -0,0 +1,66 @@
|
|
1
|
+
require 'matrix'
|
2
|
+
|
3
|
+
describe "Matrix.scalar" do
|
4
|
+
|
5
|
+
before(:each) do
|
6
|
+
@side = 3
|
7
|
+
@value = 8
|
8
|
+
@a = Matrix.scalar(@side, @value)
|
9
|
+
end
|
10
|
+
|
11
|
+
it "returns a Matrix" do
|
12
|
+
@a.should be_kind_of(Matrix)
|
13
|
+
end
|
14
|
+
|
15
|
+
it "returns a n x n matrix" do
|
16
|
+
@a.row_size.should == @side
|
17
|
+
@a.column_size.should == @side
|
18
|
+
end
|
19
|
+
|
20
|
+
it "initializes diagonal to value" do
|
21
|
+
(0...@a.row_size).each do |i|
|
22
|
+
@a[i, i].should == @value
|
23
|
+
end
|
24
|
+
end
|
25
|
+
|
26
|
+
it "initializes all non-diagonal values to 0" do
|
27
|
+
(0...@a.row_size).each do |i|
|
28
|
+
(0...@a.column_size).each do |j|
|
29
|
+
if i != j
|
30
|
+
@a[i, j].should == 0
|
31
|
+
end
|
32
|
+
end
|
33
|
+
end
|
34
|
+
end
|
35
|
+
|
36
|
+
before(:each) do
|
37
|
+
@side = 3
|
38
|
+
@value = 8
|
39
|
+
@a = Matrix.scalar(@side, @value)
|
40
|
+
end
|
41
|
+
|
42
|
+
it "returns a Matrix" do
|
43
|
+
@a.should be_kind_of(Matrix)
|
44
|
+
end
|
45
|
+
|
46
|
+
it "returns a square matrix, where the first argument specifies the side of the square" do
|
47
|
+
@a.row_size.should == @side
|
48
|
+
@a.column_size.should == @side
|
49
|
+
end
|
50
|
+
|
51
|
+
it "puts the second argument in all diagonal values" do
|
52
|
+
(0...@a.row_size).each do |i|
|
53
|
+
@a[i, i].should == @value
|
54
|
+
end
|
55
|
+
end
|
56
|
+
|
57
|
+
it "fills all values not on the main diagonal with 0" do
|
58
|
+
(0...@a.row_size).each do |i|
|
59
|
+
(0...@a.column_size).each do |j|
|
60
|
+
if i != j
|
61
|
+
@a[i, j].should == 0
|
62
|
+
end
|
63
|
+
end
|
64
|
+
end
|
65
|
+
end
|
66
|
+
end
|