The Perl and Raku Conference 2025: Greenville, South Carolina - June 27-29 Learn more

PJACKLAM / Math-Matrix-0.94 / t / sub.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
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127

              
#!perl
use strict;
use warnings;
use Math::Matrix;
use Test::More tests => 20;
note("sub() with matrices");
{
    my $x = Math::Matrix -> new([[  5, -2,  4 ],
                                 [ -1,  3,  2 ]]);
    my $y = Math::Matrix -> new([[  4,  3,  2 ],
                                 [ -1, -2,  6 ]]);
    my $z = $x -> sub($y);
    is(ref($z), 'Math::Matrix', '$z is a Math::Matrix');
    is_deeply([ @$z ], [[ 1, -5,  2],
                        [ 0,  5, -4]],
              '$z has the right values');
    # Verify that modifying $z does not modify $x or $y.
    my ($nrowy, $ncoly) = $y -> size();
    for my $i (0 .. $nrowy - 1) {
        for my $j (0 .. $ncoly - 1) {
            $z -> [$i][$j] += 100;
        }
    }
    is_deeply([ @$x ], [[  5, -2,  4 ],
                        [ -1,  3,  2 ]], '$x is unmodified');
    is_deeply([ @$y ], [[  4,  3,  2 ],
                        [ -1, -2,  6 ]], '$y is unmodified');
}
note("sub() with matrix and scalar");
{
    my $x = Math::Matrix -> new([[  5, -2,  4 ],
                                 [ -1,  3,  2 ]]);
    my $y = Math::Matrix -> new([[ 3 ]]);
    my $z = $x -> sub($y);
    is(ref($z), 'Math::Matrix', '$z is a Math::Matrix');
    is_deeply([ @$z ], [[  2, -5,  1 ],
                        [ -4,  0, -1 ]],
              '$z has the right values');
    # Verify that modifying $z does not modify $x or $y.
    my ($nrowy, $ncoly) = $y -> size();
    for my $i (0 .. $nrowy - 1) {
        for my $j (0 .. $ncoly - 1) {
            $z -> [$i][$j] += 100;
        }
    }
    is_deeply([ @$x ], [[  5, -2,  4 ],
                        [ -1,  3,  2 ]], '$x is unmodified');
    is_deeply([ @$y ], [[ 3 ]], '$y is unmodified');
}
note("sub() with scalar and matrix");
{
    my $x = Math::Matrix -> new([[ 3 ]]);
    my $y = Math::Matrix -> new([[  5, -2,  4 ],
                                 [ -1,  3,  2 ]]);
    my $z = $x -> sub($y);
    is(ref($z), 'Math::Matrix', '$z is a Math::Matrix');
    is_deeply([ @$z ], [[ -2,  5, -1 ],
                        [  4,  0,  1 ]],
              '$z has the right values');
    # Verify that modifying $z does not modify $x or $y.
    my ($nrowy, $ncoly) = $y -> size();
    for my $i (0 .. $nrowy - 1) {
        for my $j (0 .. $ncoly - 1) {
            $z -> [$i][$j] += 100;
        }
    }
    is_deeply([ @$x ], [[ 3 ]], '$x is unmodified');
    is_deeply([ @$y ], [[  5, -2,  4 ],
                        [ -1,  3,  2 ]], '$y is unmodified');
}
note("sub() with empty matrices");
{
    my $x = Math::Matrix -> new([]);
    my $y = Math::Matrix -> new([]);
    my $z = $x -> sub($y);
    is(ref($z), 'Math::Matrix', '$z is a Math::Matrix');
    is_deeply([ @$z ], [], '$z has the right values');
}
note("overloading");
{
    my $x = Math::Matrix -> new([[3]]);
    my $y = $x - 4;
    is(ref($y), 'Math::Matrix', '$y is a Math::Matrix');
    is_deeply([ @$y ], [[-1]],
              '$y has the right values');
}
{
    my $x = Math::Matrix -> new([[3]]);
    my $y = 4 - $x;
    is(ref($y), 'Math::Matrix', '$y is a Math::Matrix');
    is_deeply([ @$y ], [[1]],
              '$y has the right values');
}
{
    my $x = Math::Matrix -> new([[3]]);
    $x -= 4;
    is(ref($x), 'Math::Matrix', '$x is a Math::Matrix');
    is_deeply([ @$x ], [[-1]],
              '$x has the right values');
}