Fernando Duarte

#### stacktrace for `ForwardDiff.hessian(testf,p)` MethodError: no method matching Float64(::ForwardDiff.Dual{ForwardDiff.Tag{typeof(testf), Float64}, Float64, 2}) Closest candidates are: (::Type{T})(::Real, ::RoundingMode) where T

Thank you for your quick response, much appreciated! I think the code you link to will take me a while to understand. I will give it a shot, not sure...

Hi David! Thanks for your reply. I am working with @smitch151 on the same project. You are right that we were passing a scalar as the input to the `(Turing)Dirichlet`...

I realized the issue is now quite different. I will open a new issue and maybe @smitch151 can close this one?

+1 on fixing this issue, just run into it today

Thank you!

Maybe related: ```julia Symbolics.scalarize( (z[3:4]') * A) # error ``` ``` KeyError: key 1 not found Stacktrace: [1] getindex(h::Dict{Sym, Any}, key::Int64) @ Base ./dict.jl:482 [2] (::Symbolics.var"#93#96"{Dict{Sym, Any}})(::Tuple{Int64, Int64}) @ Symbolics...