Add tests for AbstractSteadyStateProblem default dispatches

This commit adds comprehensive tests for the new AbstractSteadyStateProblem
dispatches, verifying that:
- Both in-place and out-of-place SteadyStateProblems work with default algorithm
- The init and solve interfaces work correctly
- StaticArrays are supported
- Problem conversion from SteadyStateProblem to NonlinearProblem works
- Various problem configurations are handled correctly

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

Co-Authored-By: Claude <[email protected]>

Author: ChrisRackauckas

src/default.jl

@@ -38,31 +38,15 @@ function SciMLBase.__solve(prob::NonlinearProblem, ::Nothing, args...; kwargs...
 end
 
 function SciMLBase.__init(prob::SciMLBase.AbstractSteadyStateProblem, ::Nothing, args...; kwargs...)
-    must_use_jacobian = Val(SciMLBase.has_jac(prob.f))
-    return SciMLBase.__init(
-        prob,
-        FastShortcutNonlinearPolyalg(
-            eltype(prob.u0); must_use_jacobian, u0_len = length(prob.u0)
-        ),
-        args...;
-        kwargs...
-    )
+    # Convert SteadyStateProblem to NonlinearProblem and use its default
+    nlprob = SciMLBase.NonlinearProblem(prob)
+    return SciMLBase.__init(nlprob, nothing, args...; kwargs...)
 end
 
 function SciMLBase.__solve(prob::SciMLBase.AbstractSteadyStateProblem, ::Nothing, args...; kwargs...)
-    must_use_jacobian = Val(SciMLBase.has_jac(prob.f))
-    prefer_simplenonlinearsolve = Val(prob.u0 isa StaticArray)
-    return SciMLBase.__solve(
-        prob,
-        FastShortcutNonlinearPolyalg(
-            eltype(prob.u0);
-            must_use_jacobian,
-            prefer_simplenonlinearsolve,
-            u0_len = length(prob.u0)
-        ),
-        args...;
-        kwargs...
-    )
+    # Convert SteadyStateProblem to NonlinearProblem and use its default
+    nlprob = SciMLBase.NonlinearProblem(prob)
+    return SciMLBase.__solve(nlprob, nothing, args...; kwargs...)
 end
 
 function SciMLBase.__init(prob::NonlinearLeastSquaresProblem, ::Nothing, args...; kwargs...)

test/default_alg_tests.jl

@@ -0,0 +1,109 @@
+@testitem "Default Algorithm for AbstractSteadyStateProblem" tags=[:core] begin
+    using SciMLBase, StaticArrays
+
+    # Test with in-place function
+    function f_iip(du, u, p, t)
+        du[1] = 2 - 2u[1]
+        du[2] = u[1] - 4u[2]
+    end
+    
+    u0 = zeros(2)
+    prob_iip = SteadyStateProblem(f_iip, u0)
+    
+    @testset "In-place SteadyStateProblem" begin
+        # Test with default algorithm (nothing)
+        sol = solve(prob_iip)
+        @test SciMLBase.successful_retcode(sol.retcode)
+        @test maximum(abs, sol.resid) < 1e-6
+        
+        # Test with explicit nothing
+        sol = solve(prob_iip, nothing)
+        @test SciMLBase.successful_retcode(sol.retcode)
+        @test maximum(abs, sol.resid) < 1e-6
+        
+        # Test init interface
+        cache = init(prob_iip)
+        sol = solve!(cache)
+        @test SciMLBase.successful_retcode(sol.retcode)
+        @test maximum(abs, sol.resid) < 1e-6
+        
+        # Test init with nothing
+        cache = init(prob_iip, nothing)
+        sol = solve!(cache)
+        @test SciMLBase.successful_retcode(sol.retcode)
+        @test maximum(abs, sol.resid) < 1e-6
+    end
+    
+    # Test with out-of-place function
+    f_oop(u, p, t) = [2 - 2u[1], u[1] - 4u[2]]
+    u0 = zeros(2)
+    prob_oop = SteadyStateProblem(f_oop, u0)
+    
+    @testset "Out-of-place SteadyStateProblem" begin
+        # Test with default algorithm (nothing)
+        sol = solve(prob_oop)
+        @test SciMLBase.successful_retcode(sol.retcode)
+        @test maximum(abs, sol.resid) < 1e-6
+        
+        # Test with explicit nothing
+        sol = solve(prob_oop, nothing)
+        @test SciMLBase.successful_retcode(sol.retcode)
+        @test maximum(abs, sol.resid) < 1e-6
+        
+        # Test init interface
+        cache = init(prob_oop)
+        sol = solve!(cache)
+        @test SciMLBase.successful_retcode(sol.retcode)
+        @test maximum(abs, sol.resid) < 1e-6
+        
+        # Test init with nothing
+        cache = init(prob_oop, nothing)
+        sol = solve!(cache)
+        @test SciMLBase.successful_retcode(sol.retcode)
+        @test maximum(abs, sol.resid) < 1e-6
+    end
+    
+    # Test that SteadyStateProblem conversion works
+    @testset "Problem conversion" begin
+        # Create equivalent NonlinearProblem
+        function f_nl(u, p)
+            [2 - 2u[1], u[1] - 4u[2]]
+        end
+        
+        prob_nl = NonlinearProblem(f_nl, u0)
+        
+        # Convert SteadyStateProblem to NonlinearProblem
+        prob_converted = NonlinearProblem(prob_oop)
+        
+        # Both should solve to the same solution
+        sol_nl = solve(prob_nl)
+        sol_converted = solve(prob_converted)
+        
+        @test sol_nl.u ≈ sol_converted.u atol=1e-10
+    end
+    
+    # Test with StaticArrays
+    @testset "StaticArrays support" begin
+        f_static(u, p, t) = @SVector [2 - 2u[1], u[1] - 4u[2]]
+        u0_static = @SVector [0.0, 0.0]
+        prob_static = SteadyStateProblem(f_static, u0_static)
+        
+        sol = solve(prob_static)
+        @test SciMLBase.successful_retcode(sol.retcode)
+        @test maximum(abs, sol.resid) < 1e-6
+    end
+    
+    # Test that solve works with various problem types
+    @testset "Mixed problem types" begin
+        # Regular arrays
+        prob1 = SteadyStateProblem(f_oop, [0.5, 0.5])
+        sol1 = solve(prob1)
+        @test SciMLBase.successful_retcode(sol1.retcode)
+        
+        # With parameters
+        f_param(u, p, t) = [p[1] - 2u[1], u[1] - 4u[2]]
+        prob2 = SteadyStateProblem(f_param, [0.5, 0.5], [2.0])
+        sol2 = solve(prob2)
+        @test SciMLBase.successful_retcode(sol2.retcode)
+    end
+end