diff --git a/src/implementations/lq.jl b/src/implementations/lq.jl
index 07ae7136..bef88d4f 100644
--- a/src/implementations/lq.jl
+++ b/src/implementations/lq.jl
@@ -164,9 +164,10 @@ function lq_householder!(
     end
 
     if computeL
-        L̃ = lowertriangular!(view(A, axes(L)...))
-        positive && gaugefix!(lq_householder!, L̃, Q, diagview(A))
-        copyto!(L, L̃)
+        # we need to first copy then gaugefix - avoiding aliasing between L and Ld for broadcast
+        Ld = diagview(A)
+        copyto!(L, lowertriangular!(view(A, axes(L)...)))
+        positive && gaugefix!(lq_householder!, L, Q, Ld)
     else
         gaugefix!(lq_householder!, nothing, Q, diagview(A))
     end