Example program for zggqrf

This example solves the general Gauss–Markov linear model problem
minx y2   subject to   d=Ax+By  
where
A = 0.96-0.81i -0.03+0.96i -0.91+2.06i -0.98+1.98i -1.20+0.19i -0.66+0.42i 0.62-0.46i 1.01+0.02i 0.63-0.17i 1.08-0.28i 0.20-0.12i -0.07+1.23i ,  
B = 0.5-1.0i 0.0i+0.0 0.0i+0.0 0.0i+0.0 0.0i+0.0 1.0-2.0i 0.0i+0.0 0.0i+0.0 0.0i+0.0 0.0i+0.0 2.0-3.0i 0.0i+0.0 0.0i+0.0 0.0i+0.0 0.0i+0.0 5.0-4.0i  
and
d = 6.00-0.40i -5.27+0.90i 2.72-2.13i -1.30-2.80i .  
The solution is obtained by first computing a generalized QR factorization of the matrix pair A,B. The example illustrates the general solution process, although the above data corresponds to a simple weighted least squares problem.
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.