-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathlinalg.180433.out
108 lines (107 loc) · 4.01 KB
/
linalg.180433.out
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
A =
[[ 2.+2.j 2.+2.j 0.+0.j 0.+0.j 0.+0.j]
[ 2.+2.j 2.-2.j -4.-4.j 0.+0.j 0.+0.j]
[ 0.+0.j 2.-2.j 2.-2.j 8.+4.j 0.+0.j]
[ 0.+0.j 0.+0.j 2.+1.j 8.+2.j 4.-2.j]
[ 0.+0.j 0.+0.j 0.+0.j 4.-2.j 4.+1.j]]
Ainv =
[[ 0.1753-0.5632j 0.0747+0.3132j -0.6379-0.4885j 0.7701+0.2414j -0.8046+0.3448j]
[ 0.0747+0.3132j -0.0747-0.3132j 0.6379+0.4885j -0.7701-0.2414j 0.8046-0.3448j]
[ 0.2443-0.319j -0.2443+0.319j -0.0747-0.5632j 0.2644+0.5057j -0.5747-0.2299j]
[-0.0302+0.0963j 0.0302-0.0963j 0.0661+0.1264j -0.0287-0.1782j 0.1494+0.1264j]
[-0.0431-0.1006j 0.0431+0.1006j -0.1437-0.0575j 0.1494+0.1264j 0.023 -0.0575j]]
b =
[[ 1.+0.j 0.+0.j 0.-0.j -0.+0.j 0.-0.j]
[ 0.+0.j 1.+0.j 0.+0.j 0.+0.j 0.-0.j]
[ 0.+0.j 0.+0.j 1.-0.j 0.+0.j 0.+0.j]
[ 0.+0.j -0.+0.j -0.-0.j 1.+0.j -0.-0.j]
[ 0.-0.j -0.+0.j 0.-0.j 0.-0.j 1.+0.j]]
numpy result:
[[ 0.1753-0.5632j 0.0747+0.3132j -0.6379-0.4885j 0.7701+0.2414j -0.8046+0.3448j]
[ 0.0747+0.3132j -0.0747-0.3132j 0.6379+0.4885j -0.7701-0.2414j 0.8046-0.3448j]
[ 0.2443-0.319j -0.2443+0.319j -0.0747-0.5632j 0.2644+0.5057j -0.5747-0.2299j]
[-0.0302+0.0963j 0.0302-0.0963j 0.0661+0.1264j -0.0287-0.1782j 0.1494+0.1264j]
[-0.0431-0.1006j 0.0431+0.1006j -0.1437-0.0575j 0.1494+0.1264j 0.023 -0.0575j]]
tr = (18+1j)
A =
[[ 2.+2.j 2.+2.j 0.+0.j 0.+0.j 0.+0.j]
[ 2.+2.j 2.-2.j -4.-4.j 0.+0.j 0.+0.j]
[ 0.+0.j 2.-2.j 2.-2.j 8.+4.j 0.+0.j]
[ 0.+0.j 0.+0.j 2.+1.j 8.+2.j 4.-2.j]
[ 0.+0.j 0.+0.j 0.+0.j 4.-2.j 4.+1.j]]
A* =
[[ 2.-2.j 2.-2.j 0.-0.j 0.-0.j 0.-0.j]
[ 2.-2.j 2.+2.j 2.+2.j 0.-0.j 0.-0.j]
[ 0.-0.j -4.+4.j 2.+2.j 2.-1.j 0.-0.j]
[ 0.-0.j 0.-0.j 8.-4.j 8.-2.j 4.+2.j]
[ 0.-0.j 0.-0.j 0.-0.j 4.+2.j 4.-1.j]]
[[ True True True True True]
[ True True True True True]
[ True True True True True]
[ True True True True True]
[ True True True True True]]
numpy time = 0.0001239776611328125
my time = 0.300947904586792
[[ True True True True True]
[ True True True True True]
[ True True True True True]
[ True True True True True]
[ True True True True True]]
numpy time = 4.315376281738281e-05
my time = 0.001973867416381836
[[ True True True True True]
[ True True True True True]
[ True True True True True]
[ True True True True True]
[ True True True True True]]
numpy time = 2.9087066650390625e-05
my time = 0.0018906593322753906
[[ True True True True True]
[ True True True True True]
[ True True True True True]
[ True True True True True]
[ True True True True True]]
numpy time = 2.288818359375e-05
my time = 0.0019621849060058594
[[ True True True True True]
[ True True True True True]
[ True True True True True]
[ True True True True True]
[ True True True True True]]
numpy time = 2.288818359375e-05
my time = 0.0018897056579589844
[[ True True True True True]
[ True True True True True]
[ True True True True True]
[ True True True True True]
[ True True True True True]]
numpy time = 2.193450927734375e-05
my time = 0.0018849372863769531
[[ True True True True True]
[ True True True True True]
[ True True True True True]
[ True True True True True]
[ True True True True True]]
numpy time = 2.2411346435546875e-05
my time = 0.002045869827270508
[[ True True True True True]
[ True True True True True]
[ True True True True True]
[ True True True True True]
[ True True True True True]]
numpy time = 2.1696090698242188e-05
my time = 0.0018770694732666016
[[ True True True True True]
[ True True True True True]
[ True True True True True]
[ True True True True True]
[ True True True True True]]
numpy time = 2.193450927734375e-05
my time = 0.0021581649780273438
[[ True True True True True]
[ True True True True True]
[ True True True True True]
[ True True True True True]
[ True True True True True]]
numpy time = 2.2172927856445312e-05
my time = 0.0020952224731445312