Incompatibility with CuArrays and ForwardDiff

[charleskawczynski]

MWE:

using CUDA
using NonlinearSolve
f(u,p) = u .* u .- 2
u0 = CUDA.CuArray([1.0, 1.0])
probN = NonlinearProblem{false}(f, u0)
solver = solve(probN, NewtonRaphson(), tol = 1e-9)