Added linear algbra methods

Author: j-fu

src/extendable.jl

@@ -324,27 +324,51 @@ end
 """
 $(TYPEDSIGNATURES)
 
-Drop in replacement for LU factorization.
-
+Delegating LU factorization.
 """
 function LinearAlgebra.lu(E::ExtendableSparseMatrix)
     @inbounds flush!(E)
     return LinearAlgebra.lu(E.cscmatrix)
 end
 
+"""
+$(SIGNATURES)
 
+Delegating Matrix multiplication
+"""
+function  LinearAlgebra.mul!(r::AbstractArray{T,1} where T, E::ExtendableSparse.ExtendableSparseMatrix, x::AbstractArray{T,1} where T)
+    @inbounds flush!(E)
+    return LinearAlgebra.mul!(r,E.cscmatrix,x)
+end
 
 
 """
 $(SIGNATURES)
 
-Flush and delegate to cscmatrix.
+Delegating Matrix ldiv
 """
-function LinearAlgebra.mul!(r::AbstractArray{Tv,1},
-                            E::ExtendableSparse.ExtendableSparseMatrix{Tv,Ti},
-                            x::AbstractArray{Tv,1}) where{Tv,Ti<:Integer}
+function  LinearAlgebra.ldiv!(r::AbstractArray{T,1} where T, E::ExtendableSparse.ExtendableSparseMatrix, x::AbstractArray{T,1} where T)
+    @inbounds flush!(E)
+    return LinearAlgebra.ldiv!(r,E.cscmatrix,x)
+end
+
+"""
+$(SIGNATURES)
+
+Delegating Matrix multiplication
+"""
+function  LinearAlgebra.mul!(r::AbstractArray{T,2} where T, E::ExtendableSparse.ExtendableSparseMatrix, x::AbstractArray{T,2} where T)
     @inbounds flush!(E)
     return LinearAlgebra.mul!(r,E.cscmatrix,x)
 end
 
 
+"""
+$(SIGNATURES)
+
+Delegating Matrix ldiv
+"""
+function  LinearAlgebra.ldiv!(r::AbstractArray{T,2} where T, E::ExtendableSparse.ExtendableSparseMatrix, x::AbstractArray{T,2} where T)
+    @inbounds flush!(E)
+    return LinearAlgebra.ldiv!(r,E.cscmatrix,x)
+end

test/ExtendableSparseTest.jl

@@ -78,4 +78,24 @@ function constructors()
 
     true
 end
+
+
+function optest(n)
+    println("optest:")
+    A=ExtendableSparseMatrix(n,n)
+    nrow=min(div(n,4),5)
+    @time for i=1:n
+        A[i,i]=2*n
+        for k=1:nrow
+            j=rand((1:n))
+            if i!=j
+                A[i,j]=-rand()
+            end
+        end
+    end
+    b=rand(n)
+    @time x=A\b
+    Ax=A*x
+    b≈ Ax
+end
 end

test/runtests.jl

@@ -29,5 +29,7 @@ include("ExtendableSparseTest.jl")
         @test ExtendableSparseTest.check(m=m,n=n,nnz=nnz,nsplice=nsplice)
     end
 
-    
+    @test ExtendableSparseTest.optest(10)
+    @test ExtendableSparseTest.optest(100)
+    @test ExtendableSparseTest.optest(1000)
 end