The Perl Toolchain Summit 2025 Needs You: You can help 🙏 Learn more

PJACKLAM / Math-Matrix-0.94 / t / kron.t 

            

              
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

              
#!perl
use strict;
use warnings;
use Math::Matrix;
use Test::More tests => 16;
note('two 2-by-2 matrices');
{
    my $x = Math::Matrix -> new([[ 1, -2],
                                 [-1,  0]]);
    my $y = Math::Matrix -> new([[ 4, -3],
                                 [ 2,  3]]);
    my $z = $x -> kron($y);
    is(ref($z), 'Math::Matrix', '$z is a Math::Matrix');
    is_deeply([ @$z ], [[ 4, -3, -8,  6],
                        [ 2,  3, -4, -6],
                        [-4,  3,  0,  0],
                        [-2, -3,  0,  0]],
              '$z has the right values');
    # Verify that modifying $y does not modify $x or $y.
    my ($nrowz, $ncolz) = $z -> size();
    for my $i (0 .. $nrowz - 1) {
        for my $j (0 .. $ncolz - 1) {
            $z -> [$i][$j] += 100;
        }
    }
    is_deeply([ @$x ], [[ 1, -2],
                        [-1,  0]], '$x is unmodified');
    is_deeply([ @$y ], [[ 4, -3],
                        [ 2,  3]], '$y is unmodified');
}
note('2-by-2 matrix and empty matrix');
{
    my $x = Math::Matrix -> new([[ 1, -2],
                                 [-1,  0]]);
    my $y = Math::Matrix -> new([]);
    my $z = $x -> kron($y);
    is(ref($z), 'Math::Matrix', '$z is a Math::Matrix');
    is_deeply([ @$z ], [], '$z has the right values');
    # Verify that $x or $y are not modified.
    is_deeply([ @$x ], [[ 1, -2],
                        [-1,  0]], '$x is unmodified');
    is_deeply([ @$y ], [], '$y is unmodified');
}
note('empty matrix and 2-by-2 matrix');
{
    my $x = Math::Matrix -> new([]);
    my $y = Math::Matrix -> new([[ 4, -3],
                                 [ 2,  3]]);
    my $z = $x -> kron($y);
    is(ref($z), 'Math::Matrix', '$z is a Math::Matrix');
    is_deeply([ @$z ], [], '$z has the right values');
    # Verify that $x or $y are not modified.
    is_deeply([ @$x ], [], '$x is unmodified');
    is_deeply([ @$y ], [[ 4, -3],
                        [ 2,  3]], '$y is unmodified');
}
note('two empty matrices');
{
    my $x = Math::Matrix -> new([]);
    my $y = Math::Matrix -> new([]);
    my $z = $x -> kron($y);
    is(ref($z), 'Math::Matrix', '$z is a Math::Matrix');
    is_deeply([ @$z ], [], '$z has the right values');
    # Verify that $x or $y are not modified.
    is_deeply([ @$x ], [], '$x is unmodified');
    is_deeply([ @$y ], [], '$y is unmodified');
}