Add LinearAlgebra.lu etc. for ExtendableSparse

Author: j-fu

src/extendable.jl

@@ -270,25 +270,22 @@ if USE_GPL_LIBS
 
 for (Tv) in (:Float64,:ComplexF64)
     @eval begin
-"""
-$(TYPEDSIGNATURES)
-
-[`\\`](@ref) for ExtendableSparse for $($Tv)
-"""
         function LinearAlgebra.:\(ext::ExtendableSparseMatrix{$Tv,Ti}, B::AbstractVecOrMat{$Tv}) where Ti
             flush!(ext)
             ext.cscmatrix\B
         end
     end
 end
-
+    
 end # USE_GPL_LIBS
 
 
 """
-$(TYPEDSIGNATURES)
+    A\b
 
-[`\\`](@ref) for ExtendableSparse for generic floating point. This calls Sparspak.jl.
+[`\\`](@ref) for ExtendableSparse. It calls the LU factorization form Sparspak.jl, unless GPL components
+are allowed  in the Julia sysimage and the floating point type of the matrix is Float64 or Complex64.
+In that case, Julias standard `\` is called, which is realized via UMFPACK.
 """
 function LinearAlgebra.:\(ext::ExtendableSparseMatrix{Tv,Ti}, b::AbstractVector{Tv}) where {Tv,Ti}
     flush!(ext)

src/factorizations.jl

@@ -80,7 +80,29 @@ If `nothing` is passed as first parameter, [`factorize`](@ref) is called.
 """
 LinearAlgebra.lu!(lufact::AbstractFactorization, A::ExtendableSparseMatrix)=factorize!(lufact,A)
 
-LinearAlgebra.lu(A::ExtendableSparseMatrix)=factorize!(LUFactorization(),A)
+
+
+"""
+```
+lu(matrix)
+```
+
+Create LU factorization. It calls the LU factorization form Sparspak.jl, unless GPL components
+are allowed  in the Julia sysimage and the floating point type of the matrix is Float64 or Complex64.
+In that case, Julias standard `lu` is called, which is realized via UMFPACK.
+"""
+LinearAlgebra.lu(A::ExtendableSparseMatrix{Tv,Ti}) where {Tv, Ti} =factorize!(SparspakLU{Tv,Ti}(),A)
+
+if USE_GPL_LIBS
+
+for (Tv) in (:Float64,:ComplexF64)
+    @eval begin
+        LinearAlgebra.lu(A::ExtendableSparseMatrix{$Tv,Ti}) where Ti =factorize!(LUFactorization{$Tv,Ti}(),A)
+    end
+end
+
+end # USE_GPL_LIBS
+
 
 
 """

src/ilu0.jl

@@ -71,7 +71,7 @@ function update!(precon::ILU0Preconditioner{Tv,Ti}) where {Tv,Ti}
 end
 
 
-function  LinearAlgebra.ldiv!(u::AbstractArray{T,1}, precon::ILU0Preconditioner, v::AbstractArray{T,1}) where T
+function  LinearAlgebra.ldiv!(u::AbstractArray{T,1}, precon::ILU0Preconditioner{Tv,Ti}, v::AbstractArray{T,1}) where {T,Tv,Ti}
     cscmatrix=precon.A.cscmatrix
     colptr=cscmatrix.colptr
     rowval=cscmatrix.rowval
@@ -81,15 +81,15 @@ function  LinearAlgebra.ldiv!(u::AbstractArray{T,1}, precon::ILU0Preconditioner,
     idiag=precon.idiag
 
     @inbounds for j=1:n
-        x=zero(T)
+        x=zero(Tv)
         @inbounds for k=colptr[j]:idiag[j]-1
             x+=nzval[k]*u[rowval[k]]
         end
         u[j]=xdiag[j]*(v[j]-x)
     end
 
     @inbounds for j=n:-1:1
-        x=zero(T)
+        x=zero(Tv)
         @inbounds for k=idiag[j]+1:colptr[j+1]-1
             x+=u[rowval[k]]*nzval[k]
         end