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 = [[4.13467 0.0126659 -1.9657 3.09283 0.826586]
[0.0126659 2.33374 -0.224443 -2.06036 0.0580411]
[-1.9657 -0.224443 1.34721 -1.1825 -1.06629]
[3.09283 -2.06036 -1.1825 5.31376 1.73426]
[0.826586 0.0580411 -1.06629 1.73426 5.01824]]
norm(A*V-V*diag(d)) = 0
Real non-symmetric matrix
A = [[0.611692 -0.534532 0.103149 0.14747 -0.166006]
[-0.0799272 1.51576 -0.86184 -1.17376 0.445168]
[-0.101654 -0.116042 -0.420697 0.114017 0.361577]
[1.02917 1.18937 0.297823 0.522963 -0.741776]
[-1.37943 -0.105928 -0.286478 -0.556702 0.485322]]
norm(A*V-V*diag(d)) = 0
Complex hermitian matrix
A = [[1.85815+0i -1.27532-1.01874i 1.1571-0.793181i 0.00889625+0.683968i -0.493022-0.532008i]
[-1.27532+1.01874i 4.45757+0i 0.955997+0.88651i 0.816075-1.57469i -0.102042+1.6868i]
[1.1571+0.793181i 0.955997-0.88651i 6.50557+0i 1.34837+0.148368i -1.02759-1.46678i]
[0.00889625-0.683968i 0.816075+1.57469i 1.34837-0.148368i 2.66655+0i -2.56785+0.166685i]
[-0.493022+0.532008i -0.102042-1.6868i -1.02759+1.46678i -2.56785-0.166685i 3.72102+0i]]
norm(A*V-V*diag(d)) = 0
Complex non-hermitian matrix
A = [[0.219214-0.141661i 0.176139-0.276432i 0.83024-1.29979i -0.291467+0.780305i 0.204412+0.926516i]
[0.100501+1.05596i 1.01785-0.561079i -0.248716-0.630004i -0.675195-0.038722i 1.29624-0.103014i]
[0.126522-0.348387i -0.23416+0.548687i -0.183664+1.09633i -0.761494+1.61663i -0.093321+0.28428i]
[-0.2209+0.166731i 0.319975+0.714349i -0.147774-0.213195i -0.294119-0.799369i -0.806121+0.623665i]
[-0.540764+0.476032i 0.045842-1.01049i -0.169277-0.432241i 0.267765-1.08489i -0.869934-0.703924i]]
norm(A*V-V*diag(d)) = 0
|
