using DiffEqOperators
using SparsityDetection, SparseArrays
nknots = 10
h = 1.0/(nknots+1)
const Δ = CenteredDifference(2, 2, h, nknots)
u0 = rand(nknots, nknots)
function st!(du,u,p,t)
du .= 0.0
du[1,:] .= Δ*Dirichlet0BC(Float64)*view(u,1,:)
return nothing
end
outp = similar(u0)
# st!(outp,u0, nothing, 0.0)
sparsity_pattern = jacobian_sparsity(st!,outp,u0, nothing, 0.0)