add sparse matrix-based functions

Author: alyst

Project.toml

@@ -6,10 +6,11 @@ version = "0.1.1"
 
 [deps]
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
+SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
 Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
 
 [compat]
-julia="1"
+julia = "1"
 
 [extras]
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"

src/DiffTests.jl

@@ -1,6 +1,7 @@
 module DiffTests
 
 using LinearAlgebra: det, norm, dot, tr, Diagonal, LowerTriangular, UpperTriangular
+using SparseArrays: sparse
 using Statistics: mean
 
 #=
@@ -272,15 +273,23 @@ diag_lmul(x::VecOrMat) = diag_matrix(x) * x
 dense_lmul(x::VecOrMat) = test_matrix(x) * x
 utriag_lmul(x::VecOrMat) =  UpperTriangular(test_matrix(x)) * x
 ltriag_lmul(x::VecOrMat) =  LowerTriangular(test_matrix(x)) * x
+sparse_lmul(x::VecOrMat) = sparse(test_matrix(x)) * x
+sp_utriag_lmul(x::VecOrMat) = UpperTriangular(sparse(test_matrix(x))) * x
+sp_ltriag_lmul(x::VecOrMat) = LowerTriangular(sparse(test_matrix(x))) * x
 
 # left division by a constant matrix
 diag_ldiv(x::VecOrMat) = diag_matrix(x) \ x
 dense_ldiv(x::VecOrMat) = test_matrix(x) \ x
 utriag_ldiv(x::VecOrMat) =  UpperTriangular(test_matrix(x)) \ x
 ltriag_ldiv(x::VecOrMat) =  LowerTriangular(test_matrix(x)) \ x
+sparse_ldiv(x::VecOrMat) = sparse(test_matrix(x)) \ x
+sp_utriag_ldiv(x::VecOrMat) = UpperTriangular(sparse(test_matrix(x))) \ x
+sp_ltriag_ldiv(x::VecOrMat) = LowerTriangular(sparse(test_matrix(x))) \ x
 
 const VECTOR_TO_VECTOR_FUNCS = (diag_lmul, dense_lmul, utriag_lmul, ltriag_lmul,
-                                diag_ldiv, utriag_ldiv, ltriag_ldiv)
+                                sparse_lmul, sp_utriag_lmul, sp_ltriag_lmul,
+                                diag_ldiv, utriag_ldiv, ltriag_ldiv,
+                                sparse_ldiv, sp_utriag_ldiv, sp_ltriag_ldiv,)
 
 ######################
 # f(x::Array)::Array #
@@ -309,6 +318,8 @@ const ARRAY_TO_ARRAY_FUNCS = (-, chebyquad, brown_almost_linear, trigonometric,
 
 const MATRIX_TO_MATRIX_FUNCS = (inv,
                                 diag_lmul, dense_lmul, utriag_lmul, ltriag_lmul,
-                                diag_ldiv, utriag_ldiv, ltriag_ldiv)
+                                sparse_lmul, sp_utriag_lmul, sp_ltriag_lmul,
+                                diag_ldiv, utriag_ldiv, ltriag_ldiv,
+                                sparse_ldiv, sp_utriag_ldiv, sp_ltriag_ldiv,)
 
 end # module