function R = qr_givens_full (A) %QR_GIVENS_FULL Givens-rotation QR factorization, for full matrices. % Example: % R = qr_givens_full (A) % See also: cs_demo % CXSparse, Copyright (c) 2006-2022, Timothy A. Davis. All Rights Reserved. % SPDX-License-Identifier: LGPL-2.1+ [m n] = size (A) ; for i = 2:m for k = 1:min(i-1,n) A ([k i],k:n) = givens2 (A(k,k), A(i,k)) * A ([k i],k:n) ; A (i,k) = 0 ; end end R = A ;