#example from comment#7 in https://savannah.gnu.org/bugs/?51950
A=sparse(5,6);
A(3,1)=0.8;
A(2,2)=1.4;
A(1,6)=-0.5;
[Q,R,E]=qr(A);
P=eye(size(E,1))(:,E);
tmp = qr(R);
R
Q
Q*R*P'
Q*R*P' - A
norm(Q*R*P' -A,1)
#output
#{
R =
Compressed Column Sparse (rows = 5, cols = 6, nnz = 3 [10%])

  (1, 1) ->  0.80000
  (2, 2) ->  1.4000
  (3, 3) -> -0.50000

Q =

   0   0   1   0   0
   0   1   0   0   0
   1   0   0   0   0
   0   0   0   0   1
   0   0   0   1   0

ans =

   0.00000   0.00000   0.00000   0.00000   0.00000  -0.50000
   0.00000   1.40000   0.00000   0.00000   0.00000   0.00000
   0.80000   0.00000   0.00000   0.00000   0.00000   0.00000
   0.00000   0.00000   0.00000   0.00000   0.00000   0.00000
   0.00000   0.00000   0.00000   0.00000   0.00000   0.00000

ans =

   0   0   0   0   0   0
   0   0   0   0   0   0
   0   0   0   0   0   0
   0   0   0   0   0   0
   0   0   0   0   0   0

ans = 0
#}


#example from original submission in https://savannah.gnu.org/bugs/?51950
for iter = 1:100 
  A = sprandn(50,100,0.1);
  [Q,R,E] = qr(A);
  tmp = qr(R);
  P = eye(size(E,1))(:,E);
  norm(Q*R*P' -A, 1)
endfor
#output
#{
ans = 0
ans =    1.4271e-14
ans =    1.1788e-14
ans =    2.0439e-14
ans =    1.0154e-14
ans =    1.9055e-14
ans =    1.4896e-14
ans =    1.2781e-14
ans =    1.2125e-14
ans =    1.4393e-14
ans =    1.0944e-14
ans =    1.0554e-14
ans =    4.0611e-15
ans =    2.1929e-14
ans =    1.4191e-14
ans =    1.0631e-14
ans =    1.0893e-14
ans =    1.6126e-14
ans =    1.0326e-14
ans =    1.3363e-14
ans =    1.0637e-14
ans =    1.1998e-14
ans =    1.6871e-14
ans =    1.0460e-14
ans =    1.0773e-14
ans =    1.1534e-14
ans =    1.1489e-14
ans =    1.2182e-14
ans =    1.3023e-14
ans =    1.2662e-14
ans =    1.2145e-14
ans =    1.2032e-14
ans =    1.6777e-14
ans =    9.7472e-15
ans =    8.8055e-15
ans =    1.3443e-14
ans =    1.1380e-14
ans =    9.2924e-15
ans =    1.1688e-14
ans =    1.1191e-14
ans =    6.8161e-15
ans =    1.0483e-14
ans =    1.4870e-14
ans =    1.2284e-14
ans =    9.6079e-15
ans =    1.1169e-14
ans =    8.6018e-15
ans =    1.0529e-14
ans =    1.6749e-14
ans =    8.5819e-15
ans =    1.2312e-14
ans =    8.7115e-15
ans =    1.3746e-14
ans =    1.2676e-14
ans =    1.4547e-14
ans =    1.1168e-14
ans =    1.2189e-14
ans =    9.2493e-15
ans =    1.6768e-14
ans =    1.5188e-14
ans =    1.2573e-14
ans =    1.3793e-14
ans =    1.1316e-14
ans =    1.4129e-14
ans =    9.7887e-15
ans =    1.1227e-14
ans =    1.1644e-14
ans =    1.0137e-14
ans =    1.1683e-14
ans =    1.3745e-14
ans =    8.7235e-15
ans =    1.4010e-14
ans =    8.2370e-15
ans =    1.3020e-14
ans =    1.4285e-14
ans =    1.1485e-14
ans = 0
ans =    1.4934e-14
ans =    1.4685e-14
ans =    1.7177e-14
ans =    1.1796e-14
ans =    1.1626e-14
ans =    1.0720e-14
ans =    1.6754e-14
ans =    1.2484e-14
ans =    1.7783e-14
ans = 0
ans =    1.2210e-14
ans =    1.5358e-14
ans =    9.3984e-15
ans =    1.2565e-14
ans =    1.1641e-14
ans =    1.3474e-14
ans =    1.1512e-14
ans =    1.0719e-14
ans =    1.0236e-14
ans =    1.6947e-14
ans =    1.2800e-14
ans =    9.9990e-15
ans =    1.0704e-14
ans =    5.9517e-15
#}