>
> Assume lamnbert_W is defined through this repetitive procedure, will
> DualNumbers.jl work for this function? I did not succeed:
>
> julia> lambert_W(Dual(e, 1.0))
> ERROR: no method lambert_W(Dual{Float64},)
>
As Jason did with your previous function, if you remove "::Real" from the
function signature, it will work fine:
julia> lambert_W(Dual(e, 1.0))
1.0 + 0.18393972058572117du
If you find performance lacking without the type, you would parametrize the
function (see
http://docs.julialang.org/en/release-0.2/manual/methods/#parametric-methods
).
Cheers!
Kevin
