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
|
(conj mx1) =
1.0 2.0 3.0
4.0 5.0 6.0
7.0 8.0 9.0
10.0 11.0 12.0
OK
(transpose mx1) =
1.0 4.0 7.0 10.0
2.0 5.0 8.0 11.0
3.0 6.0 9.0 12.0
OK
(herm mx1) =
1.0 4.0 7.0 10.0
2.0 5.0 8.0 11.0
3.0 6.0 9.0 12.0
OK
mx4 =
-2.0+1.0i 1.0+0.0i
3.5-3.5i 1.5+1.0i
-10.0+0.0i 10.0+1.0i
OK
(conj mx4) =
-2.0-1.0i 1.0+-0.0i
3.5+3.5i 1.5-1.0i
-10.0+-0.0i 10.0-1.0i
OK
(transpose mx4) =
-2.0+1.0i 3.5-3.5i -10.0+0.0i
1.0+0.0i 1.5+1.0i 10.0+1.0i
OK
(herm mx4) =
-2.0-1.0i 3.5+3.5i -10.0+-0.0i
1.0+-0.0i 1.5-1.0i 10.0-1.0i
OK
(is-instance? mx8 (:diagonal-matrix <complex>)) = #t OK
(conj mx5) =
(diag 1.0 -1.5 2.0 0.0 5.0) OK
(transpose mx5) =
(diag 1.0 -1.5 2.0 0.0 5.0) OK
(herm mx5) =
(diag 1.0 -1.5 2.0 0.0 5.0) OK
(conj mx8) =
(diag 0.0-2.0i -10.5+3.5i 2.0-2.0i 1.0+-0.0i) OK
(transpose mx8) =
(diag 0.0+2.0i -10.5-3.5i 2.0+2.0i 1.0+0.0i) OK
(herm mx8) =
(diag 0.0-2.0i -10.5+3.5i 2.0-2.0i 1.0+-0.0i) OK
(diag-of mx1) =
(diag 1.0 5.0 9.0) OK
(diag-of mx9) =
(diag 1.0 5.0 9.0) OK
(generate-matrix 3 4 (lambda (((i1 <integer>) (i2 <integer>)) <real> pure) (if (= i1 i2) 1.0 0.0))) =
1.0 0.0 0.0 0.0
0.0 1.0 0.0 0.0
0.0 0.0 1.0 0.0
OK
(generate-diag-matrix 4 integer->real) =
(diag 0.0 1.0 2.0 3.0) OK
(matrix-map r-square mx1) =
1.0 4.0 9.0
16.0 25.0 36.0
49.0 64.0 81.0
100.0 121.0 144.0
OK
(matrix-map-w-ind (lambda (((i1 <integer>) (i2 <integer>) (r-elem <real>)) <real> pure) (if (= i1 i2) r-elem 1.0)) mx1) =
1.0 1.0 1.0
1.0 5.0 1.0
1.0 1.0 9.0
1.0 1.0 1.0
OK
(diag-matrix-map-w-ind (lambda (((i <integer>) (r-elem <real>)) <real> pure) (+ r-elem (integer->real i))) mx5) =
(diag 1.0 -0.5 4.0 3.0 9.0) OK
|
