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
|
use Test::More tests => 5;
use File::Spec;
use lib File::Spec->catfile("..","lib");
use Math::MatrixReal;
do 'funcs.pl';
$matrix = Math::MatrixReal->new_from_string(<<"MATRIX");
[ 1 0 0 0 1 ]
[ 0 2 0 0 2 ]
[ 0 0 3 0 0 ]
[ 0 0 0 4 0 ]
[ 0 0 0 0 1 ]
MATRIX
ok( $matrix->decompose_LR->rank_LR == 5, 'matrix has rank 5' );
########################
$matrix = Math::MatrixReal->new_from_string(<<"MATRIX");
[ 1 0 0 0 0 ]
[ 0 3 0 0 0 ]
[ 0 0 4 0 0 ]
[ 1 0 0 1 0 ]
[ 1 1 1 1 0 ]
MATRIX
ok( $matrix->decompose_LR->rank_LR == 4, 'matrix has rank 4' );
$matrix = Math::MatrixReal->new_from_string(<<"MATRIX");
[ 1 ]
MATRIX
ok($matrix->decompose_LR->rank_LR == 1, 'matrix has rank 1' );
$matrix->zero;
ok($matrix->decompose_LR->rank_LR == 0, 'zero matrix has rank 0' );
$matrix = Math::MatrixReal->new_from_string(<<"MATRIX");
[ 1 0 0 0 0e-3 ]
[ 0 0 0 0 0.0 ]
[ 0 0 0 0 0.0 ]
[ 0 0 0 0 0.0 ]
[ 0 0 0 0 0.0 ]
MATRIX
ok($matrix->decompose_LR->rank_LR == 1, 'matrix has rank 1' );
|
