@@ -188,13 +188,6 @@ rectangularQ(Q::LinearAlgebra.AbstractQ) = Matrix(Q)
188
188
@test mul! (c, b, q' ) ≈ b* q'
189
189
@test_throws DimensionMismatch mul! (Matrix {eltya} (I, n+ 1 , n), q, b)
190
190
191
- b = similar (a[:,1 ]); rand! (b)
192
- c = similar (a[:,1 ])
193
- d = similar (a[:,1 ])
194
- @test mul! (c, q, b) ≈ q* b
195
- @test mul! (c, q' , b) ≈ q' * b
196
- @test_throws DimensionMismatch mul! (Vector {eltya} (undef, n+ 1 ), q, b)
197
-
198
191
qra = qr (a[:,1 : n1], NoPivot ())
199
192
q, r = qra. Q, qra. R
200
193
@test rmul! (copy (squareQ (q)' ), q) ≈ Matrix (I, n, n)
@@ -212,6 +205,13 @@ rectangularQ(Q::LinearAlgebra.AbstractQ) = Matrix(Q)
212
205
@test mul! (c, b, q) ≈ b* q
213
206
@test mul! (c, b, q' ) ≈ b* q'
214
207
@test_throws DimensionMismatch mul! (Matrix {eltya} (I, n+ 1 , n), q, b)
208
+
209
+ b = similar (a[:,1 ]); rand! (b)
210
+ c = similar (a[:,1 ])
211
+ d = similar (a[:,1 ])
212
+ @test mul! (c, q, b) ≈ q* b
213
+ @test mul! (c, q' , b) ≈ q' * b
214
+ @test_throws DimensionMismatch mul! (Vector {eltya} (undef, n+ 1 ), q, b)
215
215
end
216
216
end
217
217
end
0 commit comments