trixi-framework · ranocha · Feb 7, 2025 · Feb 7, 2025 · Feb 7, 2025 · Feb 7, 2025
diff --git a/.github/workflows/Downgrade.yml b/.github/workflows/Downgrade.yml
@@ -72,7 +72,7 @@ jobs:
       - uses: julia-actions/cache@v2
       - uses: julia-actions/julia-downgrade-compat@v1
         with:
-          skip: LinearAlgebra,Printf,SparseArrays,UUIDs,DiffEqBase,DelimitedFiles,Test,Downloads,Random
+          skip: LinearAlgebra,Printf,SparseArrays,UUIDs,DelimitedFiles,Test,Downloads,Random
           projects: ., test
       - uses: julia-actions/julia-buildpkg@v1
         env:

diff --git a/Project.toml b/Project.toml
@@ -9,7 +9,6 @@ CodeTracking = "da1fd8a2-8d9e-5ec2-8556-3022fb5608a2"
 ConstructionBase = "187b0558-2788-49d3-abe0-74a17ed4e7c9"
 DataStructures = "864edb3b-99cc-5e75-8d2d-829cb0a9cfe8"
 DelimitedFiles = "8bb1440f-4735-579b-a4ab-409b98df4dab"
-DiffEqBase = "2b5f629d-d688-5b77-993f-72d75c75574e"
 DiffEqCallbacks = "459566f4-90b8-5000-8ac3-15dfb0a30def"
 Downloads = "f43a241f-c20a-4ad4-852c-f6b1247861c6"
 EllipsisNotation = "da5c29d0-fa7d-589e-88eb-ea29b0a81949"
@@ -69,8 +68,7 @@ ConstructionBase = "1.5"
 Convex = "0.16"
 DataStructures = "0.18.15"
 DelimitedFiles = "1"
-DiffEqBase = "6 - 6.143"
-DiffEqCallbacks = "2.35"
+DiffEqCallbacks = "2.35, 4"
 Downloads = "1.6"
 ECOS = "1.1.2"
 EllipsisNotation = "1.0"

diff --git a/examples/tree_1d_dgsem/elixir_advection_diffusion.jl b/examples/tree_1d_dgsem/elixir_advection_diffusion.jl
@@ -1,4 +1,4 @@
-using OrdinaryDiffEq
+using OrdinaryDiffEq, ADTypes
 using Trixi
 
 ###############################################################################
@@ -87,7 +87,8 @@ callbacks = CallbackSet(summary_callback, analysis_callback, alive_callback, sav
 # OrdinaryDiffEq's `solve` method evolves the solution in time and executes the passed callbacks
 time_int_tol = 1.0e-10
 time_abs_tol = 1.0e-10
-sol = solve(ode, KenCarp4(autodiff = false), abstol = time_abs_tol, reltol = time_int_tol,
+sol = solve(ode, KenCarp4(autodiff = AutoFiniteDiff()),
+            abstol = time_abs_tol, reltol = time_int_tol,
             save_everystep = false, callback = callbacks)
 
 # Print the timer summary

diff --git a/examples/tree_1d_dgsem/elixir_advection_diffusion_restart.jl b/examples/tree_1d_dgsem/elixir_advection_diffusion_restart.jl
@@ -1,4 +1,4 @@
-using OrdinaryDiffEq
+using OrdinaryDiffEq, ADTypes
 using Trixi
 
 ###############################################################################
@@ -22,7 +22,8 @@ callbacks = CallbackSet(summary_callback, analysis_callback, alive_callback)
 ###############################################################################
 # run the simulation
 
-sol = solve(ode, KenCarp4(autodiff = false), abstol = time_abs_tol, reltol = time_int_tol,
+sol = solve(ode, KenCarp4(autodiff = AutoFiniteDiff()),
+            abstol = time_abs_tol, reltol = time_int_tol,
             save_everystep = false, callback = callbacks)
 
 # Print the timer summary

diff --git a/examples/tree_3d_dgsem/elixir_lbm_taylor_green_vortex.jl b/examples/tree_3d_dgsem/elixir_lbm_taylor_green_vortex.jl
@@ -73,5 +73,5 @@ callbacks = CallbackSet(summary_callback,
 sol = solve(ode, CarpenterKennedy2N54(williamson_condition = false),
             dt = 1.0, # solve needs some value here but it will be overwritten by the stepsize_callback
             save_everystep = false, callback = callbacks,
-            save_start = false, alias_u0 = true);
+            save_start = false, alias = ODEAliasSpecifier(alias_u0 = true));
 summary_callback() # print the timer summary
diff --git a/test/Project.toml b/test/Project.toml
@@ -1,4 +1,5 @@
 [deps]
+ADTypes = "47edcb42-4c32-4615-8424-f2b9edc5f35b"
 Aqua = "4c88cf16-eb10-579e-8560-4a9242c79595"
 CairoMakie = "13f3f980-e62b-5c42-98c6-ff1f3baf88f0"
 Convex = "f65535da-76fb-5f13-bab9-19810c17039a"
@@ -21,6 +22,7 @@ StableRNGs = "860ef19b-820b-49d6-a774-d7a799459cd3"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 
 [compat]
+ADTypes = "0.2, 1.12"
 Aqua = "0.8"
 CairoMakie = "0.12, 0.13"
 Convex = "0.16"

diff --git a/test/test_mpi.jl b/test/test_mpi.jl
@@ -34,13 +34,13 @@ end # MPI
 
     t = 0.5
     let u = 1.0
-        @test ode_norm(u, t) ≈ OrdinaryDiffEq.ODE_DEFAULT_NORM(u, t)
+        @test ode_norm(u, t) ≈ OrdinaryDiffEq.DiffEqBase.ODE_DEFAULT_NORM(u, t)
     end
     let u = [1.0, -2.0]
-        @test ode_norm(u, t) ≈ OrdinaryDiffEq.ODE_DEFAULT_NORM(u, t)
+        @test ode_norm(u, t) ≈ OrdinaryDiffEq.DiffEqBase.ODE_DEFAULT_NORM(u, t)
     end
     let u = [SVector(1.0, -2.0), SVector(0.5, -0.1)]
-        @test ode_norm(u, t) ≈ OrdinaryDiffEq.ODE_DEFAULT_NORM(u, t)
+        @test ode_norm(u, t) ≈ OrdinaryDiffEq.DiffEqBase.ODE_DEFAULT_NORM(u, t)
     end
 end # MPI supporting functionality
 

diff --git a/test/test_parabolic_1d.jl b/test/test_parabolic_1d.jl
@@ -60,8 +60,8 @@ end
     # Create a CallbackSet to collect all callbacks such that they can be passed to the ODE solver
     callbacks = CallbackSet(summary_callback, analysis_callback, alive_callback,
                             amr_callback)
-    sol = solve(ode, KenCarp4(autodiff = false), abstol = time_abs_tol,
-                reltol = time_int_tol,
+    sol = solve(ode, KenCarp4(autodiff = AutoFiniteDiff()),
+                abstol = time_abs_tol, reltol = time_int_tol,
                 save_everystep = false, callback = callbacks)
     l2_error, linf_error = analysis_callback(sol)
     @test l2_error ≈ [6.4878111416468355e-6]