Add LinearAlgebra.:\ for Symmetric{ExtendableSparse} and Hermitian{ExtendableSparse}

Author: j-fu

Project.toml

@@ -1,7 +1,7 @@
 name = "ExtendableSparse"
 uuid = "95c220a8-a1cf-11e9-0c77-dbfce5f500b3"
 authors = ["Juergen Fuhrmann <juergen.fuhrmann@wias-berlin.de>"]
-version = "0.3.7"
+version = "0.3.9"
 
 [deps]
 DocStringExtensions = "ffbed154-4ef7-542d-bbb7-c09d3a79fcae"

src/extendable.jl

@@ -178,6 +178,7 @@ function Base.show(io::IO,::MIME"text/plain",ext::ExtendableSparseMatrix)
 end
 
 
+
 """
 $(SIGNATURES)
 
@@ -283,6 +284,32 @@ function LinearAlgebra.:\(ext::ExtendableSparseMatrix,X::Union{AbstractArray{T,1
 end
 
 
+"""
+$(SIGNATURES)
+
+[`\\`](@ref) for Symmetric{ExtenableSparse}
+"""
+function LinearAlgebra.:\(symm_ext::Symmetric{Tm, ExtendableSparseMatrix{Tm, Ti}}, B::Union{AbstractArray{T,1}, AbstractArray{T,2}} where T) where{Tm,Ti}
+    flush!(symm_ext.data)
+    symm_csc=Symmetric(symm_ext.data.cscmatrix,Symbol(symm_ext.uplo))
+    symm_csc\B
+end
+
+
+"""
+$(SIGNATURES)
+
+[`\\`](@ref) for Hermitian{ExtenableSparse}
+"""
+function LinearAlgebra.:\(symm_ext::Hermitian{Tm, ExtendableSparseMatrix{Tm, Ti}}, B::Union{AbstractArray{T,1}, AbstractArray{T,2}} where T) where{Tm,Ti}
+    flush!(symm_ext.data)
+    symm_csc=Hermitian(symm_ext.data.cscmatrix,Symbol(symm_ext.uplo))
+    symm_csc\B
+end
+
+
+
+
 """
 $(SIGNATURES)
 

test/runtests.jl

@@ -177,3 +177,47 @@ end
     @test   all(isapprox.(test_precon(ParallelJacobiPreconditioner,20,20,20), (true, 2.0406032775945658e-5), rtol=1.0e-5))
 end
 
+
+
+function test_symmetric(n,uplo)
+    A=ExtendableSparseMatrix(n,n)
+    sprand_sdd!(A)
+    b=rand(n)
+    flush!(A)
+    SA=Symmetric(A,uplo)
+    Scsc=Symmetric(A.cscmatrix,uplo)
+    SA\b≈ Scsc\b
+end
+
+
+
+@testset "symmetric" begin
+    @test test_symmetric(3,:U)
+    @test test_symmetric(3,:L)
+    @test test_symmetric(30,:U)
+    @test test_symmetric(30,:L)
+    @test test_symmetric(300,:U)
+    @test test_symmetric(300,:L)
+end
+
+
+
+function test_hermitian(n,uplo)
+    A=ExtendableSparseMatrix{ComplexF64,Int64}(n,n)
+    sprand_sdd!(A)
+    flush!(A)
+    A.cscmatrix.nzval.=(1.0+0.01im)*A.cscmatrix.nzval
+    b=rand(n)
+    HA=Hermitian(A,uplo)
+    Hcsc=Hermitian(A.cscmatrix,uplo)
+    HA\b≈ Hcsc\b
+end
+
+@testset "hermitian" begin
+    @test test_hermitian(3,:U)
+    @test test_hermitian(3,:L)
+    @test test_hermitian(30,:U)
+    @test test_hermitian(30,:L)
+    @test test_hermitian(300,:U)
+    @test test_hermitian(300,:L)
+end