Detection with complex fails

[ChrisRackauckas]

using SparsityDetection
function f(du, u, p, t)
    Q = reinterpret(Complex{Float64}, u)
    dQ = reinterpret(Complex{Float64}, du)
    dQ .= Q
end

input = rand(Float64, 4)
output = similar(input)
sparsity_pattern = sparsity!((dx, x) -> f(dx, x, [], 0), output, input)
jac = Float64.(sparse(sparsity_pattern))
colors = matrix_colors(jac)

[atombear]

sorry it took me a few hours to get to this!

[ChrisRackauckas]

@shashi might this have a quick solution?

[shashi]

So each complex element comes from 2 inputs?

I think the problem is reinterpret. It would work if you just used 1.0 + 2.0im.

[atombear]

the problem is almost certainly reinterpret... but i want to use reinterpret!

…

On Thu, Apr 2, 2020 at 9:13 PM Shashi Gowda ***@***.***> wrote: So each complex element comes from 2 inputs? I think the problem is reinterpret. It would work if you just used 1.0 + 2.0im. — You are receiving this because you commented. Reply to this email directly, view it on GitHub <#21 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AEABWVGMGM2DGVL75ZTQP3LRKVO7JANCNFSM4LF6XW6Q> .