make sure sparse matrix is constructed with eltype(u)

Author: jClugstor

lib/OrdinaryDiffEqDifferentiation/src/derivative_utils.jl

@@ -767,7 +767,7 @@ function build_J_W(alg, u, uprev, p, t, dt, f::F, jac_config, ::Type{uEltypeNoUn
                     !isnothing(jac_config) ?
                     convert.(
                         eltype(u), SparseMatrixColorings.sparsity_pattern(jac_config[1])) :
-                    spzeros(length(u), length(u))
+                    spzeros(eltype(u), length(u), length(u))
                 elseif eltype(f.sparsity) == Bool
                     convert.(eltype(u), f.sparsity)
                 else
@@ -798,7 +798,7 @@ function build_J_W(alg, u, uprev, p, t, dt, f::F, jac_config, ::Type{uEltypeNoUn
 
                 if isnothing(f.sparsity)
                     !isnothing(jac_config) ? convert.(eltype(u), SparseMatrixColorings.sparsity_pattern(jac_config[1])) :
-                    spzeros(length(u), length(u))
+                    spzeros(eltype(u), length(u), length(u))
                 elseif eltype(f.sparsity) == Bool
                     convert.(eltype(u), f.sparsity)
                 else