Adding CUTEst.jl to the optimization benchmarks

.buildkite/launch_benchmarks.yml.signature

@@ -1,2 +1,3 @@
-Salted__�c?�ު��DwL�J�/���)��|p��5+c��\~o
�.��	RR�hAZR��
-$'���� 6���!�m�^����S�j�Y�����
+Salted__�c?
+�ު��DwL�J�/���)��|p��5+c��\~o
�.��	RR�hAZR��
+$'���� 6���!�m�^����S�j�Y�����

.buildkite/run_benchmark.yml

@@ -19,8 +19,8 @@ steps:
       - JuliaCI/julia#v1:
           version: "1.10"
       - staticfloat/sandbox:
-          rootfs_url: "https://github.com/thazhemadam/openmodelica-rootfs-image/releases/download/v1.23.0/rootfs-openmodelica-v1.23.0.amd64.tar.gz"
-          rootfs_treehash: "82970243dc4f188e599a976abc20951f4aba2912"
+          rootfs_url: "https://github.com/thazhemadam/ScimlBenchmarks-rootfs-image/releases/download/rootfs-image/rootfs-script.amd64.tar.gz"
+          rootfs_treehash: "943bde824554f0afad5a939341f05e3bdb393051"
           uid: 1000
           gid: 1000
           workspaces:

benchmarks/AstroChem/Project.toml

@@ -10,5 +10,22 @@ ODEInterfaceDiffEq = "09606e27-ecf5-54fc-bb29-004bd9f985bf"
 OrdinaryDiffEq = "1dea7af3-3e70-54e6-95c3-0bf5283fa5ed"
 OrdinaryDiffEqCore = "bbf590c4-e513-4bbe-9b18-05decba2e5d8"
 Plots = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
+RecursiveFactorization = "f2c3362d-daeb-58d1-803e-2bc74f2840b4"
 Sundials = "c3572dad-4567-51f8-b174-8c6c989267f4"
 Symbolics = "0c5d862f-8b57-4792-8d23-62f2024744c7"
+
+[compat]
+BenchmarkTools = "1"
+Catalyst = "14"
+DiffEqDevTools = "2"
+LSODA = "0.7"
+LinearSolve = "3"
+NaNMath = "1"
+ODEInterface = "0.5"
+ODEInterfaceDiffEq = "3"
+OrdinaryDiffEq = "6.89"
+OrdinaryDiffEqCore = "1.9"
+Plots = "1"
+RecursiveFactorization = "0.2.23"
+Sundials = "4"
+Symbolics = "6"

benchmarks/AstroChem/astrochem.jmd

@@ -10,6 +10,7 @@ using Plots
 using Symbolics
 using DiffEqDevTools
 using Sundials, ODEInterface, ODEInterfaceDiffEq, LSODA
+using RecursiveFactorization
 ```
 
 ## Without Temperature Dynamics
@@ -383,8 +384,6 @@ refsol = solve(oprob, Rodas5P(), abstol=1e-14, reltol=1e-14)
 ```julia
 refsol = solve(oprob, Rodas5P(), abstol=1e-13, reltol=1e-13)
 
-solve(oprob, RadauIIA9(), abstol=1e-9, reltol=1e-9)
-
 # Run Benchmark
 
 abstols = 1.0 ./ 10.0 .^ (9:10)
@@ -393,17 +392,15 @@ reltols = 1.0 ./ 10.0 .^ (9:10)
 setups = [
           Dict(:alg=>FBDF()),
           Dict(:alg=>QNDF()),
-          #Dict(:alg=>Rodas4P()),
           Dict(:alg=>CVODE_BDF()),
           #Dict(:alg=>ddebdf()),
-          #Dict(:alg=>Rodas4()),
           Dict(:alg=>Rodas5P()),
           Dict(:alg=>KenCarp4()),
           Dict(:alg=>KenCarp47()),
-          Dict(:alg=>RadauIIA9()),
+          #Dict(:alg=>RadauIIA9()),
           #Dict(:alg=>rodas()),
           #Dict(:alg=>radau()),
-          #Dict(:alg=>lsoda()),
+          Dict(:alg=>lsoda()),
           #Dict(:alg=>ImplicitEulerExtrapolation(min_order = 5, init_order = 3,threading = OrdinaryDiffEqCore.PolyesterThreads())),
           #Dict(:alg=>ImplicitEulerExtrapolation(min_order = 5, init_order = 3,threading = false)),
           #Dict(:alg=>ImplicitEulerBarycentricExtrapolation(min_order = 5, threading = OrdinaryDiffEqCore.PolyesterThreads())),