NAG Toolbox: nag_lapack_zgeqpf (f08bs)

Purpose

Syntax

Description

References

Parameters

Compulsory Input Parameters

1:     $\mathrm{a}\left(\mathit{lda},:\right)$ – complex array

The first dimension of the array a must be at least $\max \left(1,\mathbf{m}\right)$.

The second dimension of the array a must be at least $\max \left(1,\mathbf{n}\right)$.

The $m$ by $n$ matrix $A$.

2:     $\mathrm{jpvt}\left(:\right)$ – int64int32nag_int array

The dimension of the array jpvt must be at least $\max \left(1,\mathbf{n}\right)$

If $\mathbf{jpvt}\left(i\right)\ne 0$, then the $i$ th column of $A$ is moved to the beginning of $AP$ before the decomposition is computed and is fixed in place during the computation. Otherwise, the $i$ th column of $A$ is a free column (i.e., one which may be interchanged during the computation with any other free column).

Optional Input Parameters

1:     $\mathrm{m}$ – int64int32nag_int scalar

Default: the first dimension of the array a.

$m$, the number of rows of the matrix $A$.

Constraint: $\mathbf{m}\ge 0$.

2:     $\mathrm{n}$ – int64int32nag_int scalar

Default: the second dimension of the array a.

$n$, the number of columns of the matrix $A$.

Constraint: $\mathbf{n}\ge 0$.

Output Parameters

1:     $\mathrm{a}\left(\mathit{lda},:\right)$ – complex array

The first dimension of the array a will be $\max \left(1,\mathbf{m}\right)$.

The second dimension of the array a will be $\max \left(1,\mathbf{n}\right)$.

If $m\ge n$, the elements below the diagonal store details of the unitary matrix $Q$ and the upper triangle stores the corresponding elements of the $n$ by $n$ upper triangular matrix $R$.

If $m<n$, the strictly lower triangular part stores details of the unitary matrix $Q$ and the remaining elements store the corresponding elements of the $m$ by $n$ upper trapezoidal matrix $R$.

The diagonal elements of $R$ are real.

2:     $\mathrm{jpvt}\left(:\right)$ – int64int32nag_int array

The dimension of the array jpvt will be $\max \left(1,\mathbf{n}\right)$

Details of the permutation matrix $P$. More precisely, if $\mathbf{jpvt}\left(i\right)=k$, then the $k$th column of $A$ is moved to become the $i$ th column of $AP$; in other words, the columns of $AP$ are the columns of $A$ in the order $\mathbf{jpvt}\left(1\right),\mathbf{jpvt}\left(2\right),\dots ,\mathbf{jpvt}\left(n\right)$.

3:     $\mathrm{tau}\left(\min \left(\mathbf{m},\mathbf{n}\right)\right)$ – complex array

Further details of the unitary matrix $Q$.

4:     $\mathrm{info}$ – int64int32nag_int scalar

$\mathbf{info}=0$ unless the function detects an error (see Error Indicators and Warnings).

Error Indicators and Warnings

   $\mathbf{info}=-i$

If $\mathbf{info}=-i$, parameter $i$ had an illegal value on entry. The parameters are numbered as follows:

1: m, 2: n, 3: a, 4: lda, 5: jpvt, 6: tau, 7: work, 8: rwork, 9: info.

It is possible that info refers to a parameter that is omitted from the MATLAB interface. This usually indicates that an error in one of the other input parameters has caused an incorrect value to be inferred.

Accuracy

Further Comments

[a, info] = f08at(a(:,1:m), tau);

[a, info] = f08at(a, tau);

[c, info] = f08au('Left','Conjugate Transpose', a(:,min(m,n)), tau, c);

Example

Open in the MATLAB editor: f08bs_example

function f08bs_example


fprintf('f08bs example results\n\n');

a = [ 0.47 - 0.34i, -0.40 + 0.54i,  0.60 + 0.01i,  0.80 - 1.02i;
     -0.32 - 0.23i, -0.05 + 0.20i, -0.26 - 0.44i, -0.43 + 0.17i;
      0.35 - 0.60i, -0.52 - 0.34i,  0.87 - 0.11i, -0.34 - 0.09i;
      0.89 + 0.71i, -0.45 - 0.45i, -0.02 - 0.57i,  1.14 - 0.78i;
     -0.19 + 0.06i,  0.11 - 0.85i,  1.44 + 0.80i,  0.07 + 1.14i];
b = [-0.85 - 1.63i,  2.49 + 4.01i;
     -2.16 + 3.52i, -0.14 + 7.98i;
      4.57 - 5.71i,  8.36 - 0.28i;
      6.38 - 7.40i, -3.55 + 1.29i;
      8.41 + 9.39i, -6.72 + 5.03i];
[m,n] = size(a);
jpvt = zeros(n,1,'int64');

% Compute the QR factorization of a
[a, jpvt, tau, info] = f08bs( ...
			      a, jpvt);

% Choose tol to reflect the relative accuracy of the input data
tol = 0.01;

% Determine which columns of R to use
k = find(abs(diag(a)) <= tol*abs(a(1,1)));
if numel(k) == 0
  k = numel(diag(a));
else
  k = k(1)-1;
end

% Compute c = (q^H)*b,
[c, info] = f08au( ...
		   'Left', 'Conjugate Transpose', a, tau, b);

% Compute least-squares solution by backsubstitution in r*b = c
c(1:k, :) = inv(triu(a(1:k,1:k)))*c(1:k,:);
c(k+1:4, :) = 0;

% Unscramble the least-squares solution stored in c
x = zeros(4, 2);
for i=1:4
  x(jpvt(i), :) = c(i, :);
end

fprintf('\nLeast-squares solution\n');
disp(x);

f08bs example results


Least-squares solution
   0.0000 + 0.0000i   0.0000 + 0.0000i
   2.6925 + 8.0446i  -2.0563 - 2.9759i
   2.7602 + 2.5455i   1.0588 + 1.4635i
   2.7383 + 0.5123i  -1.4150 + 0.2982i