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

JuliaLang · Apr 11, 2017 · c5e00f0 · c5e00f0
1 parent cf5ad80
commit c5e00f0
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)