Automatic chain differential for system of equations? about forwarddiff.jl HOT 3 CLOSED

julia> using ForwardDiff

julia> f1(s) = (s/3)^2
f1 (generic function with 1 method)

julia> f2(s, l) = s^2 + exp(3*l/sqrt(4))
f2 (generic function with 1 method)

julia> function f(x)
                  return f1(x[1]) + f2(x[1], x[2])
              end
f (generic function with 1 method)

julia> g = gradient_func(f, Partials{2,Float64}, mutates=false)
gradf (generic function with 1 method)

julia> g([0.1, 56.2])
2-element Array{Float64,1}:
 0.222222
 6.12514e36

julia> g = forwarddiff_gradient(f, Float64)

julia> g = forwarddiff_gradient(f, Float64, fadtype=:typed, n=2)