/
zhgeqz_example.r
executable file
·38 lines (32 loc) · 1.82 KB
/
zhgeqz_example.r
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
ZHGEQZ Example Program Results
Matrix A after balancing
1 2 3 4
1 ( 1.0000, 3.0000) ( 1.0000, 4.0000) ( 0.1000, 0.5000) ( 0.1000, 0.6000)
2 ( 2.0000, 2.0000) ( 4.0000, 3.0000) ( 0.8000, 0.4000) ( 1.6000, 0.5000)
3 ( 0.3000, 0.1000) ( 0.9000, 0.2000) ( 0.2700, 0.0300) ( 0.8100, 0.0400)
4 ( 0.4000, 0.0000) ( 1.6000, 0.1000) ( 0.6400, 0.0200) ( 2.5600, 0.0300)
Matrix B after balancing
1 2 3 4
1 ( 1.0000, 0.0000) ( 2.0000, 1.0000) ( 0.3000, 0.2000) ( 0.4000, 0.3000)
2 ( 1.0000, 1.0000) ( 4.0000, 2.0000) ( 0.9000, 0.3000) ( 1.6000, 0.4000)
3 ( 0.1000, 0.2000) ( 0.8000, 0.3000) ( 0.2700, 0.0400) ( 0.6400, 0.0500)
4 ( 0.1000, 0.3000) ( 1.6000, 0.4000) ( 0.8100, 0.0500) ( 2.5600, 0.0600)
Matrix A in Hessenberg form
1 2 3 4
1 ( -2.868, -1.595) ( -0.809, -0.328) ( -4.900, -0.987) ( -0.048, 1.163)
2 ( -2.672, 0.595) ( -0.790, 0.049) ( -4.955, -0.163) ( -0.439, -0.574)
3 ( 0.000, 0.000) ( -0.098, -0.011) ( -1.168, -0.137) ( -1.756, -0.205)
4 ( 0.000, 0.000) ( 0.000, 0.000) ( 0.087, 0.004) ( 0.032, 0.001)
Matrix B is triangular
1 2 3 4
1 ( -1.775, 0.000) ( -0.721, 0.043) ( -5.021, 1.190) ( -0.145, 0.726)
2 ( 0.000, 0.000) ( -0.218, 0.035) ( -2.541, -0.146) ( -0.823, -0.418)
3 ( 0.000, 0.000) ( 0.000, 0.000) ( -1.396, -0.163) ( -1.747, -0.204)
4 ( 0.000, 0.000) ( 0.000, 0.000) ( 0.000, 0.000) ( -0.100, -0.004)
Minimal required LWORK = 4
Actual value of LWORK = 24
Generalized eigenvalues
1 ( -0.635, 1.653)
2 ( 0.493, 0.910)
3 ( 0.458, -0.843)
4 ( 0.674, -0.050)