Sparsity detection missing values due to temporary storage? about sparsitydetection.jl HOT 6 CLOSED

julia> const Z = Int[1,2]
2-element Array{Int64,1}:
 1
 2
# This works! --
julia> sparsity!([1,2],[2,3]) do du,u
                 Z[:] .= u[:]
                du[1] = Z[1]
              end
Explored path: SparsityDetection.Path(Bool[], 1)
2×2 SparseMatrixCSC{Bool,Int64} with 1 stored entry:
  [1, 1]  =  1

# This doesn't --
julia> sparsity!([1,2],[2,3]) do du,u
                 Z[:] .= u[:]
                du[:] .= Z[:]
              end
Explored path: SparsityDetection.Path(Bool[], 1)
2×2 SparseMatrixCSC{Bool,Int64} with 0 stored entries

const N = 12

const My = Tridiagonal([1.0 for i in 1:N-1],[-2.0 for i in 1:N],[1.0 for i in 1:N-1])
My[1,2] = 2.0
My[end,end-1] = 2.0

function f(du,u,p,t)
  mul!(du,My,u)
end

using SparsityDetection
input = rand(N,N)
output = similar(input)
sparsity_pattern = sparsity!(f,output,input,nothing,0.0)
# empty