fix: :wrappers testing

Author: avik-pal

ext/NonlinearSolveMINPACKExt.jl

@@ -28,7 +28,8 @@ function SciMLBase.__solve(
         original = MINPACK.fsolve(
             f!, u0, m; tol, show_trace, tracing, method, iterations = maxiters)
     else
-        _jac! = NonlinearSolve.__construct_extension_jac(prob, alg, u0, resid; alg.autodiff)
+        autodiff = alg.autodiff === missing ? nothing : alg.autodiff
+        _jac! = NonlinearSolve.__construct_extension_jac(prob, alg, u0, resid; autodiff)
         jac! = @closure (J, u) -> (_jac!(J, u); Cint(0))
         original = MINPACK.fsolve(
             f!, jac!, u0, m; tol, show_trace, tracing, method, iterations = maxiters)

ext/NonlinearSolveNLsolveExt.jl

@@ -18,7 +18,8 @@ function SciMLBase.__solve(
     if prob.f.jac === nothing && alg.autodiff isa Symbol
         df = OnceDifferentiable(f!, u0, resid; alg.autodiff)
     else
-        jac! = NonlinearSolve.__construct_extension_jac(prob, alg, u0, resid; alg.autodiff)
+        autodiff = alg.autodiff isa Symbol ? nothing : alg.autodiff
+        jac! = NonlinearSolve.__construct_extension_jac(prob, alg, u0, resid; autodiff)
         if prob.f.jac_prototype === nothing
             J = similar(
                 u0, promote_type(eltype(u0), eltype(resid)), length(u0), length(resid))

ext/NonlinearSolveSIAMFANLEquationsExt.jl

@@ -91,10 +91,11 @@ function SciMLBase.__solve(prob::NonlinearProblem, alg::SIAMFANLEquationsJL, arg
                         f, u, m, zeros(T, N, 2 * m + 4); atol, rtol, maxit = maxiters, beta)
                 end
             else
+                autodiff = alg.autodiff === missing ? nothing : alg.autodiff
                 FPS = prob.f.jac_prototype !== nothing ? zero(prob.f.jac_prototype) :
                       __zeros_like(u, N, N)
                 jac = NonlinearSolve.__construct_extension_jac(
-                    prob, alg, u, resid; alg.autodiff)
+                    prob, alg, u, resid; autodiff)
                 AJ! = @closure (J, u, x) -> jac(J, x)
                 if method == :newton
                     sol = nsol(f, u, FS, FPS, AJ!; sham = 1, atol,

src/algorithms/pseudo_transient.jl

@@ -17,7 +17,7 @@ This implementation specifically uses "switched evolution relaxation"
 """
 function PseudoTransient(;
         concrete_jac = nothing, linsolve = nothing, linesearch = nothing,
-        precs = DEFAULT_PRECS,  alpha_initial = 1e-3, autodiff = nothing,
+        precs = DEFAULT_PRECS, alpha_initial = 1e-3, autodiff = nothing,
         jvp_autodiff = nothing, vjp_autodiff = nothing)
     descent = DampedNewtonDescent(; linsolve, precs, initial_damping = alpha_initial,
         damping_fn = SwitchedEvolutionRelaxation())