Add tests for Ax_mul_B!(Matrix, D, Matrix)

JuliaLang · Apr 11, 2017 · 90c5cc4 · 90c5cc4
1 parent 07ffd2b
commit 90c5cc4
Showing 1 changed file with 5 additions and 0 deletions.
diff --git a/test/linalg/diagonal.jl b/test/linalg/diagonal.jl
@@ -143,6 +143,11 @@ srand(1)
         @test (r = full(D) * v   ; A_mul_B!(vv, D, v)  ≈ r ≈ vv)
         @test (r = full(D)' * v  ; Ac_mul_B!(vv, D, v) ≈ r ≈ vv)
         @test (r = full(D).' * v ; At_mul_B!(vv, D, v) ≈ r ≈ vv)
+
+        UU = similar(U)
+        @test (r = full(D) * U   ; A_mul_B!(UU, D, U) ≈ r ≈ UU)
+        @test (r = full(D)' * U  ; Ac_mul_B!(UU, D, U) ≈ r ≈ UU)
+        @test (r = full(D).' * U ; At_mul_B!(UU, D, U) ≈ r ≈ UU)
     end
     @testset "triu/tril" begin
         @test istriu(D)