[julia-users] least squares residuals from QR ldiv?

[Simon Kornblith]

Is there a way to get the least squares residuals for qrfact(X)\y without 
having to compute X*beta, as is given by gels? I believe this should be 
||Q2'y||; Q2'y appears to be computed but explicitly zeroed in A_ldiv_B!.

Simon