This example finds the basic solutions for the linear least squares problems
where
and
are the columns of the matrix
,
A
factorization is performed on the first
rows of
using
dgeqrt after which the first
rows of
are updated by applying
using
dgemqrt. The remaining row is added by performing a
update using
dtpqrt;
is updated by applying the new
using
dtpmqrt; the solution is finally obtained by triangular solve using
from the updated
.