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
|
================================
Test of eigenvalue routines
================================
Real symmetric matrix
A = [[6.75107 -4.11024 7.57353 -6.9009 0.784665]
[-4.11024 3.8051 -6.47556 4.54544 -2.45177]
[7.57353 -6.47556 13.9033 -8.85583 3.75844]
[-6.9009 4.54544 -8.85583 9.84391 -1.34372]
[0.784665 -2.45177 3.75844 -1.34372 3.26664]]
norm(A*V-V*diag(d)) = 0
Real non-symmetric matrix
A = [[-0.227649 1.07507 1.12342 0.494639 -0.179719]
[-1.83729 -0.396638 1.46072 0.721932 0.549174]
[0.0511841 0.51065 -0.269572 0.246282 -0.808968]
[1.44132 1.10323 -1.00539 0.672038 -2.29735]
[-1.48058 0.910165 1.87195 -1.00149 -0.453797]]
norm(A*V-V*diag(d)) = 0
Complex hermitian matrix
A = [[5.66035+0i -0.828828+0.624351i -0.795847+0.598926i 1.01174-1.23943i 0.243402-0.0850845i]
[-0.828828-0.624351i 4.69172+0i 0.169611-1.02905i -1.92101+1.73425i 1.5394-0.718265i]
[-0.795847-0.598926i 0.169611+1.02905i 1.48641+0i -0.549066-0.315334i 0.543691+0.138787i]
[1.01174+1.23943i -1.92101-1.73425i -0.549066+0.315334i 2.5618+0i -2.13434-0.0421936i]
[0.243402+0.0850845i 1.5394+0.718265i 0.543691-0.138787i -2.13434+0.0421936i 4.30329+0i]]
norm(A*V-V*diag(d)) = 0
Complex non-hermitian matrix
A = [[0.61333+0.506973i 0.350946-0.839171i -0.431011-0.173918i 0.43428-0.328699i 0.289013+0.366511i]
[-0.295419+0.675185i 0.530681-0.233433i -0.666422-0.180822i -0.497205+0.369685i 0.840226+0.358469i]
[-0.475898+0.120633i -0.376251+0.957435i 0.135151+0.431443i 1.30784+0.644844i -0.553286+0.0173776i]
[-0.85725-0.0647025i 0.959656+0.279574i 0.735259+0.373431i -0.64014-0.467499i 0.166865+0.139432i]
[0.564002+0.646938i -0.165729-2.07819i -0.43385-0.288087i 0.248209-0.110927i -0.304968+1.67774i]]
norm(A*V-V*diag(d)) = 0
|
