Example program for zgbrfs
This example solves the system of equations
A
X
=
B
using iterative refinement and to compute the forward and backward error bounds, where
A
=
-
1.65
+
2.26
i
-
2.05
-
0.85
i
0.97
-
2.84
i
0.00
+
0.00
i
0.00
+
6.30
i
-
1.48
-
1.75
i
-
3.99
+
4.01
i
0.59
-
0.48
i
0.00
+
0.00
i
-
0.77
+
2.83
i
-
1.06
+
1.94
i
3.33
-
1.04
i
0.00
+
0.00
i
0.00
+
0.00
i
4.48
-
1.09
i
-
0.46
-
1.72
i
and
B
=
-
1.06
+
21.50
i
12.85
+
0
2.84
i
-
22.72
-
53.90
i
-
70.22
+
21.57
i
28.24
-
38.60
i
-
20.73
-
0
1.23
i
-
34.56
+
16.73
i
26.01
+
31.97
i
.
Here
A
is nonsymmetric and is treated as a band matrix, which must first be factorized by
zgbtrf
.
f90 source code
Input data file
Expected results