Fix format and add compat

Author: tansongchen

Project.toml

@@ -65,6 +65,7 @@ SciMLBase = "2.37.0"
 Setfield = "1.1.1"
 StaticArrays = "1.9"
 StaticArraysCore = "1.4.2"
+TaylorDiff = "0.2.5"
 Test = "1.10"
 Tracker = "0.2.33"
 Zygote = "0.6.69"
@@ -88,9 +89,10 @@ ReTestItems = "817f1d60-ba6b-4fd5-9520-3cf149f6a823"
 Reexport = "189a3867-3050-52da-a836-e630ba90ab69"
 ReverseDiff = "37e2e3b7-166d-5795-8a7a-e32c996b4267"
 StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
+TaylorDiff = "b36ab563-344f-407b-a36a-4f200bebf99c"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 Tracker = "9f7883ad-71c0-57eb-9f7f-b5c9e6d3789c"
 Zygote = "e88e6eb3-aa80-5325-afca-941959d7151f"
 
 [targets]
-test = ["AllocCheck", "Aqua", "DiffEqBase", "ExplicitImports", "FiniteDiff", "ForwardDiff", "Hwloc", "InteractiveUtils", "LinearAlgebra",  "NonlinearProblemLibrary", "Pkg", "PolyesterForwardDiff", "Random", "ReTestItems", "Reexport", "ReverseDiff",  "StaticArrays", "Test", "Tracker", "Zygote"]
+test = ["AllocCheck", "Aqua", "DiffEqBase", "ExplicitImports", "FiniteDiff", "ForwardDiff", "Hwloc", "InteractiveUtils", "LinearAlgebra",  "NonlinearProblemLibrary", "Pkg", "PolyesterForwardDiff", "Random", "ReTestItems", "Reexport", "ReverseDiff",  "StaticArrays", "TaylorDiff", "Test", "Tracker", "Zygote"]

ext/SimpleNonlinearSolveTaylorDiffExt.jl

@@ -1,13 +1,15 @@
 module SimpleNonlinearSolveTaylorDiffExt
 using SimpleNonlinearSolve
-using SimpleNonlinearSolve: ImmutableNonlinearProblem, ReturnCode, build_solution, check_termination, init_termination_cache
+using SimpleNonlinearSolve: ImmutableNonlinearProblem, ReturnCode, build_solution,
+                            check_termination, init_termination_cache
 using SimpleNonlinearSolve: __maybe_unaliased, _get_fx, __fixed_parameter_function
 using MaybeInplace: @bb
 using SciMLBase: isinplace
 
 import TaylorDiff
 
-@inline function __get_higher_order_derivatives(::SimpleHouseholder{N}, prob, f, x, fx) where N
+@inline function __get_higher_order_derivatives(
+        ::SimpleHouseholder{N}, prob, f, x, fx) where {N}
     vN = Val(N)
     l = map(one, x)
     t = TaylorDiff.make_seed(x, l, vN)
@@ -28,7 +30,7 @@ end
 
 function SciMLBase.__solve(prob::ImmutableNonlinearProblem, alg::SimpleHouseholder{N},
         args...; abstol = nothing, reltol = nothing, maxiters = 1000,
-        termination_condition = nothing, alias_u0 = false, kwargs...) where N
+        termination_condition = nothing, alias_u0 = false, kwargs...) where {N}
     x = __maybe_unaliased(prob.u0, alias_u0)
     length(x) == 1 ||
         throw(ArgumentError("SimpleHouseholder only supports scalar problems"))

test/core/rootfind_tests.jl

@@ -1,6 +1,7 @@
 @testsetup module RootfindingTesting
 using Reexport
-@reexport using AllocCheck, StaticArrays, Random, LinearAlgebra, ForwardDiff, DiffEqBase
+@reexport using AllocCheck, StaticArrays, Random, LinearAlgebra, ForwardDiff, DiffEqBase,
+                TaylorDiff
 import PolyesterForwardDiff
 
 quadratic_f(u, p) = u .* u .- p
@@ -92,17 +93,17 @@ end
 end
 
 @testitem "SimpleHouseholder" setup=[RootfindingTesting] tags=[:core] begin
-    using TaylorDiff
     @testset "AutoDiff: TaylorDiff.jl" for order in (2, 3, 4)
-        @testset "[OOP] u0: $(nameof(typeof(u0)))" for u0 in (
-            [1.0], @SVector[1.0], 1.0)
-            sol = benchmark_nlsolve_oop(quadratic_f, u0; solver = SimpleHouseholder{order}())
+        @testset "[OOP] u0: $(nameof(typeof(u0)))" for u0 in ([1.0], @SVector[1.0], 1.0)
+            sol = benchmark_nlsolve_oop(
+                quadratic_f, u0; solver = SimpleHouseholder{order}())
             @test SciMLBase.successful_retcode(sol)
             @test all(abs.(sol.u .* sol.u .- 2) .< 1e-9)
         end
 
         @testset "[IIP] u0: $(nameof(typeof(u0)))" for u0 in ([1.0],)
-            sol = benchmark_nlsolve_iip(quadratic_f!, u0; solver = SimpleHouseholder{order}())
+            sol = benchmark_nlsolve_iip(
+                quadratic_f!, u0; solver = SimpleHouseholder{order}())
             @test SciMLBase.successful_retcode(sol)
             @test all(abs.(sol.u .* sol.u .- 2) .< 1e-9)
         end
@@ -112,7 +113,8 @@ end
         u0 in (1.0, [1.0], @SVector[1.0])
 
         probN = NonlinearProblem(quadratic_f, u0, 2.0)
-        @test all(solve(probN, SimpleHouseholder{2}(); termination_condition).u .≈ sqrt(2.0))
+        @test all(solve(probN, SimpleHouseholder{2}(); termination_condition).u .≈
+                  sqrt(2.0))
     end
 end