* and \: support const SparseMatrixCSC

Author: alyst

Project.toml

@@ -13,6 +13,7 @@ Preferences = "21216c6a-2e73-6563-6e65-726566657250"
 Printf = "de0858da-6303-5e67-8744-51eddeeeb8d7"
 Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 SpecialFunctions = "276daf66-3868-5448-9aa4-cd146d93841b"
+SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
 
 [weakdeps]
 StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"

src/ForwardDiff.jl

@@ -5,7 +5,9 @@ using DiffResults: DiffResult, MutableDiffResult
 using Preferences
 using Random
 using LinearAlgebra
+using SparseArrays
 using Base: require_one_based_indexing
+
 import Printf
 import NaNMath
 import SpecialFunctions

src/dual.jl

@@ -824,17 +824,20 @@ end
 # implement ldiv!() of dual vector by a matrix
 LinearAlgebra.ldiv!(y::StridedVector{T},
                     m::Union{LowerTriangular{<:LinearAlgebra.BlasFloat},
-                             UpperTriangular{<:LinearAlgebra.BlasFloat}},
+                             UpperTriangular{<:LinearAlgebra.BlasFloat},
+                             SparseMatrixCSC{<:LinearAlgebra.BlasFloat}},
                     x::StridedVector{T}) where T <: Dual =
     (ldiv!(reinterpret(reshape, valtype(T), y)', m, reinterpret(reshape, valtype(T), x)'); y)
 
 Base.:\(m::Union{LowerTriangular{<:LinearAlgebra.BlasFloat},
-                 UpperTriangular{<:LinearAlgebra.BlasFloat}},
+                 UpperTriangular{<:LinearAlgebra.BlasFloat},
+                 SparseMatrixCSC{<:LinearAlgebra.BlasFloat}},
         x::StridedVector{<:Dual}) = ldiv!(similar(x), m, x)
 
 for MT in (StridedMatrix{<:LinearAlgebra.BlasFloat},
            LowerTriangular{<:LinearAlgebra.BlasFloat},
            UpperTriangular{<:LinearAlgebra.BlasFloat},
+           SparseMatrixCSC{<:LinearAlgebra.BlasFloat},
            )
 @eval begin