What about `gamma_inc_inv`'s derivative ?

[lrnv]

I am hitting the following issue :

ERROR: MethodError: no method matching __gamma_inc_inv(::ForwardDiff.Dual{...}, ::ForwardDiff.Dual{...}, ::Bool)

Is there somewhere an implementation for such a derivative ? I guess the routines for the inverse of the gamma incomplete function are based on numerical root finding, so it that even possible ?