Remaining tests for SDiagonal

Author: mschauer

test/SDiagonal.jl

@@ -4,25 +4,42 @@
         @test SDiagonal{1,Float64}((1,)).diag === SVector{1,Float64}((1,))
 
         @test SDiagonal{4,Float64}((1, 1.0, 1, 1)).diag.data === (1.0, 1.0, 1.0, 1.0)
+        @test SDiagonal{4}((1, 1.0, 1, 1)).diag.data === (1.0, 1.0, 1.0, 1.0)
+        @test SDiagonal((1, 1.0, 1, 1)).diag.data === (1.0, 1.0, 1.0, 1.0)
 
         # Bad input
         @test_throws Exception SMatrix{1,Int}()
         @test_throws Exception SMatrix{2,Int}((1,))
 
         # From SMatrix
         @test SDiagonal(SMatrix{2,2,Int}((1,2,3,4))).diag.data === (1,4)
+        
+        @test SDiagonal{1,Int}(SDiagonal{1,Float64}((1,))).diag[1] === 1
 
     end
 
     @testset "Methods" begin
 
         @test StaticArrays.scalem(@SMatrix([1 1 1;1 1 1; 1 1 1]), @SVector [1,2,3]) === @SArray [1 2 3; 1 2 3; 1 2 3]
+        @test StaticArrays.scalem(@SVector([1,2,3]),@SMatrix [1 1 1;1 1 1; 1 1 1])' === @SArray [1 2 3; 1 2 3; 1 2 3]
 
         m = SDiagonal(@SVector [11, 12, 13, 14])
-        m2 = diagm([11, 12, 13, 14])
 
-        b = @SVector [2,-1,2,1]
-        b2 = Vector(b)
+    
+        
+        @test diag(m) === m.diag
+        
+     
+        m2 = diagm([11, 12, 13, 14])
+     
+        @test logdet(m) == logdet(m2)
+        @test logdet(im*m) ≈ logdet(im*m2)
+        @test det(m) == det(m2)
+        @test trace(m) == trace(m2)
+        @test logm(m) == logm(m2)
+        @test expm(m) == expm(m2)
+        @test sqrtm(m) == sqrtm(m2)
+     
 
         @test isimmutable(m) == true
 
@@ -55,21 +72,30 @@
 
         @test_throws Exception m[1] = 1
 
+        b = @SVector [2,-1,2,1]
+        b2 = Vector(b)
+  
+        
         @test m*b ==  @SVector [22,-12,26,14]
         @test (b'*m)' ==  @SVector [22,-12,26,14]
 
         @test m\b == m2\b
+        
+        @test b'/m == b'/m2
+        @test_throws Exception b/m
         @test m*m == m2*m
 
         @test ishermitian(m) == ishermitian(m2)
+        @test ishermitian(m/2)
+        
         @test isposdef(m) == isposdef(m2)
         @test issymmetric(m) == issymmetric(m2)
 
-        @test m' == m
+        @test (2*m/2)' == m
         @test 2m == m + m
-        @test 0m == m - m
+        @test m*0 == m - m
 
-        @test m\m == eye(SDiagonal{4,Float64})
+        @test m/m == m\m == eye(SDiagonal{4,Float64})