add test for Quaternions

Author: araujoms

test/generic.jl

@@ -123,6 +123,31 @@ end
     @test_throws DimensionMismatch axpy!(α, x, Vector(1:3), y, Vector(1:5))
 end
 
+@testset "generic syrk & herk" begin
+    for T ∈ (BigFloat, Complex{BigFloat}, Quaternion{Float64})
+        α = randn(T)
+        a = randn(T, 3, 4)
+        csmall = similar(a, 3, 3)
+        csmall_fallback = similar(a, 3, 3)
+        cbig = similar(a, 4, 4)
+        cbig_fallback = similar(a, 4, 4)
+        mul!(csmall, a, a', real(α), false)
+        LinearAlgebra._generic_matmatmul!(csmall_fallback, a, a', real(α), false)
+        @test ishermitian(csmall)
+        @test csmall ≈ csmall_fallback
+        mul!(cbig, a', a, real(α), false)
+        LinearAlgebra._generic_matmatmul!(cbig_fallback, a', a, real(α), false)
+        @test ishermitian(cbig)
+        @test cbig ≈ cbig_fallback
+        mul!(csmall, a, transpose(a), α, false)
+        LinearAlgebra._generic_matmatmul!(csmall_fallback, a, transpose(a), α, false)
+        @test csmall ≈ csmall_fallback
+        mul!(cbig, transpose(a), a, α, false)
+        LinearAlgebra._generic_matmatmul!(cbig_fallback, transpose(a), a, α, false)
+        @test cbig ≈ cbig_fallback
+    end
+end
+
 @test !issymmetric(fill(1,5,3))
 @test !ishermitian(fill(1,5,3))
 @test (x = fill(1,3); cross(x,x) == zeros(3))

test/matmul.jl

@@ -537,22 +537,6 @@ end
     @test_throws DimensionMismatch LinearAlgebra.herk_wrapper!(A5x5, 'N', A6x5)
 end
 
-@testset "generic syrk & herk" begin
-    for T ∈ (BigFloat, Complex{BigFloat})
-        a = randn(T, 3, 4)
-        csmall = similar(a, 3, 3)
-        cbig = similar(a, 4, 4)
-        LinearAlgebra._generic_matmatmul!(csmall, a, a', true, false)
-        @test csmall ≈ a * a'
-        LinearAlgebra._generic_matmatmul!(csmall, a, transpose(a), true, false)
-        @test csmall ≈ a * transpose(a)
-        LinearAlgebra._generic_matmatmul!(cbig, a', a, true, false)
-        @test cbig ≈ a' * a
-        LinearAlgebra._generic_matmatmul!(cbig, transpose(a), a, true, false)
-        @test cbig ≈ transpose(a) * a
-    end
-end
-
 @testset "matmul for types w/o sizeof (issue #1282)" begin
     AA = fill(complex(1, 1), 10, 10)
     for A in (copy(AA), view(AA, 1:10, 1:10))