fix test

Author: jClugstor

test/interface/sparsediff_tests.jl

@@ -85,19 +85,15 @@ for f in [f_oop, f_ip]
     end
 end
 
+# test for https://github.com/SciML/OrdinaryDiffEq.jl/issues/2653#issuecomment-2778430025
+
 function sparse_f!(du, u, p, t)
     du[1] = u[1] + u[2]
     du[2] = u[3]^2
     return du[3] = u[1]^2
 end
 
-backend_allow = AutoSparse(
-    AutoForwardDiff();
-    sparsity_detector = TracerSparsityDetector(),
-    coloring_algorithm = GreedyColoringAlgorithm(; allow_denser = true)
-)
-
-backend_no_allow = AutoSparse(
+backend = AutoSparse(
     AutoForwardDiff();
     sparsity_detector = TracerSparsityDetector(),
     coloring_algorithm = GreedyColoringAlgorithm()
@@ -107,33 +103,19 @@ u = ones(3)
 du = zero(u)
 p = t = nothing
 
-prep_allow = DI.prepare_jacobian(
-    sparse_f!, du, backend_allow, u, DI.Constant(p), DI.Constant(t))
-prep_no_allow = DI.prepare_jacobian(
-    sparse_f!, du, backend_no_allow, u, DI.Constant(p), DI.Constant(t))
+prep = DI.prepare_jacobian(
+    sparse_f!, du, backend, u, DI.Constant(p), DI.Constant(t))
 # this is what the user may typically provide to the ODE problem
 
-function inplace_jac_allow!(J, u, p, t)
+function inplace_jac!(J, u, p, t)
     return DI.jacobian!(
-        sparse_f!, zeros(3), J, prep_allow, backend_allow, u, DI.Constant(p), DI.Constant(t))
+        sparse_f!, zeros(3), J, prep, backend, u, DI.Constant(p), DI.Constant(t))
 end
 
-function inplace_jac_no_allow!(J, u, p, t)
-    return DI.jacobian!(
-        sparse_f!, zeros(3), J, prep_no_allow, backend_no_allow, u, DI.Constant(p), DI.Constant(t))
-end
-
-jac_prototype = similar(sparsity_pattern(prep_allow), eltype(u))
-
-ode_f_allow = ODEFunction(
-    sparse_f!, jac = inplace_jac_allow!, jac_prototype = jac_prototype)
-prob_allow = ODEProblem(ode_f_allow, [1, 1, 1], (0.0, 1.0))
-
-ode_f_no_allow = ODEFunction(
-    sparse_f!, jac = inplace_jac_no_allow!, jac_prototype = jac_prototype)
-prob_no_allow = ODEProblem(ode_f_no_allow, [1, 1, 1], (0.0, 1.0))
+jac_prototype = similar(sparsity_pattern(prep), eltype(u))
 
-sol = solve(prob_allow, Rodas5())
+ode_f = ODEFunction(sparse_f!, jac = inplace_jac!, jac_prototype = jac_prototype)
+prob = ODEProblem(ode_f, [1, 1, 1], (0.0, 1.0))
 
-@test_throws DimensionMismatch sol=solve(prob_no_allow, Rodas5())
+@test_no_warn sol = solve(prob, Rodas5())