similar doesn't preserve the sparsity pattern

JuliaDiff · Oct 6, 2023 · b71300d · b71300d
1 parent 04810f5
commit b71300d
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 1 deletion.
diff --git a/Project.toml b/Project.toml
@@ -1,7 +1,7 @@
 name = "SparseDiffTools"
 uuid = "47a9eef4-7e08-11e9-0b38-333d64bd3804"
 authors = ["Pankaj Mishra <[email protected]>", "Chris Rackauckas <[email protected]>"]
-version = "2.7.0"
+version = "2.8.0"
 
 [deps]
 ADTypes = "47edcb42-4c32-4615-8424-f2b9edc5f35b"

diff --git a/src/highlevel/common.jl b/src/highlevel/common.jl
@@ -260,6 +260,7 @@ function init_jacobian(c::AbstractMaybeSparseJacobianCache)
 end
 init_jacobian(::Nothing, ::Type{T}, fx, x) where {T} = similar(fx, T, length(fx), length(x))
 init_jacobian(J, ::Type{T}, _, _) where {T} = similar(J, T, size(J, 1), size(J, 2))
+init_jacobian(J::SparseMatrixCSC, ::Type{T}, _, _) where {T} = T.(J)
 
 __maybe_copy_x(_, x) = x
 __maybe_copy_x(_, ::Nothing) = nothing