It seems like this should work. All that is needed is: ```jl Base.big(x::Dual{T}) where {T} = Dual{T}(big(x.value), big(x.partials)) Base.big(p::Partials) = Partials(big.(p.values)) ``` (Mentioned [on discourse](https://discourse.julialang.org/t/optim-jl-cant-optimize-function-which-includes-bigfloat-variables/65095).)