add a missing method for transpose(A::Hermitian{<:Real}) (#22280)

JuliaLang · Jun 8, 2017 · 2d8c0f9 · 2d8c0f9
1 parent 0061f2a
commit 2d8c0f9
Show file tree

Hide file tree

Showing 2 changed files with 12 additions and 5 deletions.
diff --git a/base/linalg/symmetric.jl b/base/linalg/symmetric.jl
@@ -210,6 +210,7 @@ issymmetric(A::Hermitian{<:Real}) = true
 issymmetric(A::Hermitian{<:Complex}) = isreal(A)
 issymmetric(A::Symmetric) = true
 transpose(A::Symmetric) = A
+transpose(A::Hermitian{<:Real}) = A
 ctranspose(A::Symmetric{<:Real}) = A
 function ctranspose(A::Symmetric)
     AC = ctranspose(A.data)

diff --git a/test/linalg/symmetric.jl b/test/linalg/symmetric.jl
@@ -93,14 +93,20 @@ let n=10
             @test ishermitian(Symmetric(b + b'))
         end
 
-        #transpose, ctranspose
+        # transpose, ctranspose
+        S = Symmetric(asym)
+        H = Hermitian(asym)
         if eltya <: Real
-            @test transpose(Symmetric(asym)) == asym
+            @test  transpose(S) === S == asym
+            @test ctranspose(S) === S == asym
+            @test  transpose(H) === H == asym
+            @test ctranspose(H) === H == asym
         else
-            @test transpose(Hermitian(asym)) == transpose(asym)
+            @test  transpose(S) === S
+            @test ctranspose(S) ==  Symmetric(conj(asym))
+            @test  transpose(H) ==  Hermitian(transpose(asym))
+            @test ctranspose(H) === H == asym
         end
-        @test ctranspose(Symmetric(asym)) == Symmetric(conj(asym))
-        @test ctranspose(Hermitian(asym)) == asym
 
         #tril/triu
         for di in -n:n