Add shallow types and KLU support, initial UMFPACK support as well (#90)

* enable use of Plain Ol' Arrays using shallow pack

* port the KLU experiments, expand pack interface slightly

* klu working

* klu working

* KLU can now use GBVector, add a copyto!, real and imag functions extended to collections

* followup to prev commit

* broaden broadcasting definitions, prep for internal pack

* add shallow types, rm gbpointer

* finish rm gbpoiinter

* enable refactor, pack and unpack improvements + packbroadcast

* add klu to testing set up. flesh out shallow types (note: need to add repack interface)

* fix #88

* pull in Umfpack work, untested except by hand

* rm extraneous kwargs

Author: Will Kimmerer

Project.toml

@@ -22,7 +22,9 @@ SuiteSparse = "4607b0f0-06f3-5cda-b6b1-a6196a1729e9"
 ChainRulesCore = "1"
 MacroTools = "0.5"
 Preferences = "1"
-SpecialFunctions = "2"
 SSGraphBLAS_jll = "7.2"
+SpecialFunctions = "2"
 StorageOrders = "0.2"
 julia = "1.6"
+KLU = "0.4"
+

lib/libinternal.jl

@@ -23,6 +23,7 @@ using SpecialFunctions: lgamma, gamma, erf, erfc
 using Base.Broadcast
 using Serialization
 using StorageOrders
+using KLU
 
 export ColMajor, RowMajor, storageorder #reexports from StorageOrders
 include("abstracts.jl")
@@ -31,7 +32,7 @@ include("libutils.jl")
 include("../lib/LibGraphBLAS_gen.jl")
 using .LibGraphBLAS
 
-include("operators/libgbops.jl")
+include("operators/libgbops.jl") 
 
 include("gbtypes.jl")
 include("types.jl")
@@ -61,7 +62,12 @@ include("vector.jl")
 include("matrix.jl")
 include("abstractgbarray.jl")
 include("random.jl")
-# 
+# Miscellaneous Operations
+include("print.jl")
+include("pack.jl")
+include("unpack.jl")
+include("options.jl")
+# Core operations (mul, elementwise, etc)
 include("operations/operationutils.jl")
 include("operations/transpose.jl")
 include("operations/mul.jl")
@@ -74,10 +80,6 @@ include("operations/concat.jl")
 include("operations/resize.jl")
 include("operations/sort.jl")
 # 
-include("print.jl")
-include("pack.jl")
-include("unpack.jl")
-include("options.jl")
 
 include("operations/broadcasts.jl")
 include("chainrules/chainruleutils.jl")
@@ -95,6 +97,9 @@ include("misc.jl")
 include("mmread.jl")
 # include("iterator.jl")
 include("oriented.jl")
+include("solvers/klu.jl")
+include("solvers/umfpack.jl")
+include("shallowtypes.jl")
 
 export SparseArrayCompat
 export LibGraphBLAS