The reason is because this is incorrect during differentiation: you aren’t using Float64 values all of the time. Instead you’re looking for A = Array{eltype(u),2}(undef,3,3).
2 Likes
The reason is because this is incorrect during differentiation: you aren’t using Float64 values all of the time. Instead you’re looking for A = Array{eltype(u),2}(undef,3,3).