Example program for zgelsd
This example solves the linear least squares problem
min
x
b
-
A
x
2
for the solution,
x
, of minimum norm, where
A
=
0.47
-
0.34
i
-
0.32
-
0.23
i
0.35
-
0.60
i
0.89
+
0.71
i
-
0.19
+
0.06
i
-
0.40
+
0.54
i
-
0.05
+
0.20
i
-
0.52
-
0.34
i
-
0.45
-
0.45
i
0.11
-
0.85
i
0.60
+
0.01
i
-
0.26
-
0.44
i
0.87
-
0.11
i
-
0.02
-
0.57
i
1.44
+
0.80
i
0.80
-
1.02
i
-
0.43
+
0.17
i
-
0.34
-
0.09
i
1.14
-
0.78
i
0.07
+
1.14
i
and
b
=
2.15
-
0.20
i
-
2.24
+
1.82
i
4.45
-
4.28
i
5.70
-
6.25
i
.
A tolerance of
0.01
is used to determine the effective rank of
A
.
Note that the block size (NB) of
64
assumed in this example is not realistic for such a small problem, but should be suitable for large problems.
f90 source code
Input data file
Expected results