Fix SpeedMapping 0.4 compatibility

SpeedMapping 0.4 introduced breaking API changes:
- `converged::Bool` -> `status::Symbol` in return type
- `tol` -> `abstol`
- `Lp` -> `pnorm`
- `store_info` -> `store_trace`
- `σ_min` -> `initial_learning_rate`
- `orders::Vector{Int}` -> `orders::Tuple`

This commit adds version detection to support both 0.3 and 0.4+ APIs.

🤖 Generated with [Claude Code](https://claude.com/claude-code)

Co-Authored-By: Claude Opus 4.5 <noreply@anthropic.com>

Author: ChrisRackauckas

ext/NonlinearSolveSpeedMappingExt.jl

@@ -1,11 +1,15 @@
 module NonlinearSolveSpeedMappingExt
 
-using SpeedMapping: speedmapping
+using SpeedMapping: SpeedMapping, speedmapping
 
 using NonlinearSolveBase: NonlinearSolveBase
 using NonlinearSolve: NonlinearSolve, SpeedMappingJL
 using SciMLBase: SciMLBase, NonlinearProblem, ReturnCode
 
+# Check SpeedMapping version to determine API
+const SPEEDMAPPING_VERSION = pkgversion(SpeedMapping)
+const SPEEDMAPPING_NEW_API = SPEEDMAPPING_VERSION >= v"0.4"
+
 function SciMLBase.__solve(
         prob::NonlinearProblem, alg::SpeedMappingJL, args...;
         abstol = nothing, maxiters = 1000, alias = SciMLBase.NonlinearAliasSpecifier(alias_u0 = false),
@@ -28,18 +32,37 @@ function SciMLBase.__solve(
 
     time_limit = ifelse(maxtime === nothing, 1000, maxtime)
 
-    sol = speedmapping(
-        u; m!, tol, Lp = Inf, maps_limit = maxiters, alg.orders,
-        alg.check_obj, store_info = store_trace isa Val{true}, alg.σ_min, alg.stabilize,
-        time_limit
-    )
+    sol = if SPEEDMAPPING_NEW_API
+        # SpeedMapping 0.4+ API
+        speedmapping(
+            u; m!, abstol = tol, pnorm = Inf, maps_limit = maxiters,
+            orders = Tuple(alg.orders),
+            initial_learning_rate = alg.σ_min == 0.0 ? 1.0 : alg.σ_min,
+            store_trace = store_trace isa Val{true},
+            time_limit
+        )
+    else
+        # SpeedMapping 0.3 API
+        speedmapping(
+            u; m!, tol, Lp = Inf, maps_limit = maxiters, alg.orders,
+            alg.check_obj, store_info = store_trace isa Val{true}, alg.σ_min, alg.stabilize,
+            time_limit
+        )
+    end
     res = prob.u0 isa Number ? first(sol.minimizer) : sol.minimizer
     resid = NonlinearSolveBase.Utils.evaluate_f(prob, res)
 
+    # SpeedMapping 0.3 returns `converged::Bool`, 0.4+ returns `status::Symbol`
+    converged = if hasproperty(sol, :converged)
+        sol.converged
+    else
+        sol.status == :first_order
+    end
+
     return SciMLBase.build_solution(
         prob, alg, res, resid;
         original = sol, stats = SciMLBase.NLStats(sol.maps, 0, 0, 0, sol.maps),
-        retcode = ifelse(sol.converged, ReturnCode.Success, ReturnCode.Failure)
+        retcode = ifelse(converged, ReturnCode.Success, ReturnCode.Failure)
     )
 end