Given a QR factorization of a real or complex matrix A = Q*R, Q unitary and R upper trapezoidal, return the QR factorization of A + u*v', where u and v are column vectors (rank-1 update).
If the matrix Q is not square, the matrix A is updated by Q*Q'*u*v' instead of u*v'.