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
|
mx1 =
(diag 1.0 1.0 1.0 1.0 1.0) OK
(diagonal-matrix-ref mx5 0) = 3.0 OK
(matrix-ref mx5 1 1) = 4.0 OK
(matrix-ref mx5 0 1) = 0.0 OK
mx5 =
(diag 3.0 4.0 2.0 2.0 2.0) OK
(number-of-rows mx2) = 4 OK
(number-of-columns mx2) = 4 OK
(diag-matrix-numel mx2) = 4 OK
(diag->normal-matrix mx2) =
2.0 0.0 0.0 0.0
0.0 -1.0 0.0 0.0
0.0 0.0 2.5 0.0
0.0 0.0 0.0 2.0
OK
(is-instance? mx6 <bvm-matrix>) = #t OK
(->real-diag-matrix mx2) =
(diag 2.0 -1.0 2.5 2.0) OK
(->single-diag-matrix mx2) =
(diag 2.0 -1.0 2.5 2.0) OK
(->complex-diag-matrix mx2) =
(diag 2.0+0.0i -1.0+0.0i 2.5+0.0i 2.0+0.0i) OK
(->single-complex-diag-matrix mx2) =
(diag 2.0+0.0i -1.0+0.0i 2.5+0.0i 2.0+0.0i) OK
(+ mx2 mx3) =
(diag 3.0 1.0 5.0 -1.0) OK
(+ mx2 mx7) =
3.0 2.0 3.0 1.0
4.0 4.0 6.0 -1.0
7.0 8.0 11.5 10.0
-1.0 -2.0 -3.0 -2.0
OK
(+ mx7 mx2) =
3.0 2.0 3.0 1.0
4.0 4.0 6.0 -1.0
7.0 8.0 11.5 10.0
-1.0 -2.0 -3.0 -2.0
OK
(- mx2 mx3) =
(diag 1.0 -3.0 0.0 5.0) OK
(- mx2 mx7) =
1.0 -2.0 -3.0 -1.0
-4.0 -6.0 -6.0 1.0
-7.0 -8.0 -6.5 -10.0
1.0 2.0 3.0 6.0
OK
(- mx7 mx2) =
-1.0 2.0 3.0 1.0
4.0 6.0 6.0 -1.0
7.0 8.0 6.5 10.0
-1.0 -2.0 -3.0 -6.0
OK
(- mx2) =
(diag -2.0 1.0 -2.5 -2.0) OK
(* mx2 mx3) =
(diag 2.0 -2.0 6.25 -6.0) OK
(* mx2 mx4) =
2.0 4.0 6.0
-4.0 -5.0 -6.0
17.5 20.0 22.5
-2.0 -4.0 -6.0
OK
(* (transpose mx4) mx2) =
2.0 -4.0 17.5 -2.0
4.0 -5.0 20.0 -4.0
6.0 -6.0 22.5 -6.0
OK
(* 5.0 mx2) =
(diag 10.0 -5.0 12.5 10.0) OK
(* mx2 5.0) =
(diag 10.0 -5.0 12.5 10.0) OK
(= mx2 mx2) = #t OK
(= mx2 mx3) = #f OK
(= (diag->normal-matrix mx2) mx2) = #t OK
(= (diag->normal-matrix mx3) mx2) = #f OK
(conj mx2) =
(diag 2.0 -1.0 2.5 2.0) OK
(transpose mx2) =
(diag 2.0 -1.0 2.5 2.0) OK
(herm mx2) =
(diag 2.0 -1.0 2.5 2.0) OK
(diag-of mx4) =
(diag 1.0 5.0 9.0) OK
(bvm-generate-diag-matrix 3 (lambda (((i <integer>)) <real> pure) (- (integer->real i)))) =
(diag -0.0 -1.0 -2.0) OK
(matrix-map r-square mx2) =
(diag 4.0 1.0 6.25 4.0) OK
(diag-matrix-map-w-ind (lambda (((i <integer>) (r-elem <real>)) <real> pure) (+ (integer->real i) r-elem)) mx2) =
(diag 2.0 0.0 4.5 5.0) OK
(+ mx2 (constant-real-bvm-diag-matrix (quote (1.0 2.0))))
exception add-op-dm-dm:dim-mismatch OK
(- mx2 (constant-real-bvm-diag-matrix (quote (1.0 2.0))))
exception add-op-dm-dm:dim-mismatch OK
(+ mx2 mx4)
exception add-op-dm-mx:dim-mismatch OK
(+ mx4 mx2)
exception add-op-dm-mx:dim-mismatch OK
(- mx2 mx4)
exception add-op-dm-mx:dim-mismatch OK
(- mx4 mx2)
exception add-op-mx-dm:dim-mismatch OK
(* mx2 (constant-real-bvm-diag-matrix (quote (1.0 2.0))))
exception mult-op-dm-dm:dim-mismatch OK
(* mx2 (constant-real-bvm-matrix (quote ((1.0 2.0) (3.0 4.0)))))
exception mult-op-dm-mx:dim-mismatch OK
(* (constant-real-bvm-matrix (quote ((1.0 2.0) (3.0 4.0)))) mx2)
exception mult-op-mx-dm:dim-mismatch OK
|
