Update test_blockdiagonal.jl

 - Stop using `Symmtric` in test. Instead use dense matrices.
 - Adding tests for matrix functions.

Author: mipals

test/test_blockdiagonal.jl

@@ -2,13 +2,13 @@ using BlockDiagonalMatrices
 using LinearAlgebra
 using SparseArrays
 using StaticArrays
+using Test
 
 function get_random_sized_psdm(n)
     b1 = rand(n,n);
     Q = qr(b1).Q
     B1 = Q*Diagonal((rand(n) .+ 0.5 ))*Q'
-    B1 = Symmetric(B1)
-    return B1
+    return (B1 + B1')/2
 end
 
 #### Testing equal block sizes
@@ -39,6 +39,18 @@ SB = sparse(B)
 @test B[1,1] == SB[1,1]
 @test B[end,end] == SB[end,end]
 
+B_dense = Matrix(B)
+b1 = size(B.blocks[1],1)
+bn = size(B.blocks[end],1) - 1
+for func in (:log, :sqrt, :sin, :tan, :cos, :sinh, :tanh)
+    @eval begin
+        func_dense = ($func)(B_dense)
+        func_block = ($func)(B)
+        @test func_dense[1:b1,1:b1] ≈ func_block.blocks[1]
+        @test func_dense[end-bn:end,end-bn:end] ≈ func_block.blocks[end]
+    end
+end
+
 
 # Make 3x3 blocks - But with StaticArrays for speedup
 C = BlockDiagonal([SMatrix{3,3,eltype(block)}(block) for block in B.blocks])
@@ -75,6 +87,17 @@ Xr = randn(size(Br,1),3)
 @test A*Xr ≈ Br*Xr
 @test A\Xr ≈ Br\Xr
 
+A_dense = Matrix(Br)
+r1 = size(Br.blocks[1],1)
+rn = size(Br.blocks[end],1) - 1
+for func in (:log, :sqrt, :sin, :tan, :cos, :sinh, :tanh)
+    @eval begin
+        func_dense = ($func)(A_dense)
+        func_block = ($func)(Br)
+        @test func_dense[1:r1,1:r1] ≈ func_block.blocks[1]
+        @test func_dense[end-rn:end,end-rn:end] ≈ func_block.blocks[end]
+    end
+end
 
 # Testing logs, determinants, and variants
 @test logdet(A) ≈ logdet(Br)