Updates

Author: ParamThakkar123

benchmarks/DynamicalODE/Henon-Heiles_energy_conservation_benchmark.jmd

@@ -19,22 +19,20 @@ In order to separate these two, we will use `iip` for the in-place names and `oo
 using OrdinaryDiffEq, Plots, DiffEqCallbacks
 using SciMLBenchmarks
 using TaylorIntegration, LinearAlgebra, StaticArrays
-
-# Set plotting formats
 gr(fmt=:png)
 default(fmt=:png)
 
 T(p) = 1//2 * norm(p)^2
-V(q) = 1//2 * (q[1]^2 + q[2]^2 + 2*q[1]^2 * q[2] - 2//3 * q[2]^3)
-H(p, q, params) = T(p) + V(q)
+V(q) = 1//2 * (q[1]^2 + q[2]^2 + 2q[1]^2 * q[2]- 2//3 * q[2]^3)
+H(p,q, params) = T(p) + V(q)
 
-function iip_dq(dq, p, q, params, t)
+function iip_dq(dq,p,q,params,t)
     dq[1] = p[1]
     dq[2] = p[2]
 end
 
-function iip_dp(dp, p, q, params, t)
-    dp[1] = -q[1] * (1 + 2 * q[2])
+function iip_dp(dp,p,q,params,t)
+    dp[1] = -q[1] * (1 + 2q[2])
     dp[2] = -q[2] - (q[1]^2 - q[2]^2)
 end
 
@@ -46,7 +44,7 @@ function oop_dq(p, q, params, t)
 end
 
 function oop_dp(p, q, params, t)
-    dp1 = -q[1] * (1 + 2 * q[2])
+    dp1 = -q[1] * (1 + 2q[2])
     dp2 = -q[2] - (q[1]^2 - q[2]^2)
     @SVector [dp1, dp2]
 end
@@ -58,33 +56,19 @@ function hamilton(du, u, p, t)
     dq, q = @views u[3:4], du[3:4]
     dp, p = @views u[1:2], du[1:2]
 
-    dp[1] = -q[1] * (1 + 2 * q[2])
+    dp[1] = -q[1] * (1 + 2q[2])
     dp[2] = -q[2] - (q[1]^2 - q[2]^2)
     dq .= p
 
     return nothing
 end
 
-function manifold_jacobian!(J, u)
-    p, q = u[1:2], u[3:4]
-    
-    J[1, 1] = p[1]
-    J[1, 2] = p[2]
-    
-    J[1, 3] = q[1] + 4 * q[1] * q[2]
-    J[1, 4] = q[2] + 2 * q[1]^2 - 2 * q[2]^2
-    
-    J[2:4, :] .= 0
-    
-    nothing
-end
-
 function g(resid, u, p)
     resid[1] = H([u[1], u[2]], [u[3], u[4]], nothing) - E
     resid[2:4] .= 0
 end
 
-const cb = ManifoldProjection(g; nlsolve = (f,u0,p) -> (res = nlsolve(f, u0; ftol=1e-13); res.zero), manifold_jacobian=manifold_jacobian!)
+const cb = ManifoldProjection(g, nlopts=Dict(:ftol => 1e-13))
 const E = H(iip_p0, iip_q0, nothing)
 ```
 

benchmarks/DynamicalODE/Manifest.toml

@@ -1,6 +1,6 @@
 # This file is machine-generated - editing it directly is not advised
 
-julia_version = "1.10.5"
+julia_version = "1.10.4"
 manifest_format = "2.0"
 project_hash = "268db6928ccffd4d6f79409a36ddb0a44c29d12b"
 
@@ -385,15 +385,10 @@ version = "6.158.3"
     Unitful = "1986cc42-f94f-5a68-af5c-568840ba703d"
 
 [[deps.DiffEqCallbacks]]
-deps = ["DataStructures", "DiffEqBase", "ForwardDiff", "Functors", "LinearAlgebra", "Markdown", "NonlinearSolve", "Parameters", "RecipesBase", "RecursiveArrayTools", "SciMLBase", "StaticArraysCore"]
-git-tree-sha1 = "19dbd44d18bbfdfcf5e56c99cea9b0ed23df350a"
+deps = ["ConcreteStructs", "DataStructures", "DiffEqBase", "DifferentiationInterface", "Functors", "LinearAlgebra", "Markdown", "RecipesBase", "RecursiveArrayTools", "SciMLBase", "StaticArraysCore"]
+git-tree-sha1 = "7f700fa4fb6e55f4672f8218ef228107245a2e9d"
 uuid = "459566f4-90b8-5000-8ac3-15dfb0a30def"
-version = "3.9.1"
-
-    [deps.DiffEqCallbacks.weakdeps]
-    OrdinaryDiffEq = "1dea7af3-3e70-54e6-95c3-0bf5283fa5ed"
-    OrdinaryDiffEqCore = "bbf590c4-e513-4bbe-9b18-05decba2e5d8"
-    Sundials = "c3572dad-4567-51f8-b174-8c6c989267f4"
+version = "4.0.0"
 
 [[deps.DiffEqNoiseProcess]]
 deps = ["DiffEqBase", "Distributions", "GPUArraysCore", "LinearAlgebra", "Markdown", "Optim", "PoissonRandom", "QuadGK", "Random", "Random123", "RandomNumbers", "RecipesBase", "RecursiveArrayTools", "ResettableStacks", "SciMLBase", "StaticArraysCore", "Statistics"]
@@ -469,9 +464,9 @@ version = "0.6.16"
 
 [[deps.DispatchDoctor]]
 deps = ["MacroTools", "Preferences"]
-git-tree-sha1 = "810cae95ca4db1c563be2d7d38fa9b7289d302db"
+git-tree-sha1 = "453df2ce2aef7de59c69a56d31dcd2ec3384dd77"
 uuid = "8d63f2c5-f18a-4cf2-ba9d-b3f60fc568c8"
-version = "0.4.16"
+version = "0.4.17"
 weakdeps = ["ChainRulesCore", "EnzymeCore"]
 
     [deps.DispatchDoctor.extensions]
@@ -779,9 +774,9 @@ version = "0.73.8+0"
 
 [[deps.GenericLinearAlgebra]]
 deps = ["LinearAlgebra", "Printf", "Random", "libblastrampoline_jll"]
-git-tree-sha1 = "f47136cac29a9b7a8a88dbce1195394978091edb"
+git-tree-sha1 = "c4f9c87b74aedf20920034bd4db81d0bffc527d2"
 uuid = "14197337-ba66-59df-a3e3-ca00e7dcff7a"
-version = "0.3.13"
+version = "0.3.14"
 
 [[deps.GenericSchur]]
 deps = ["LinearAlgebra", "Printf"]
@@ -2720,7 +2715,7 @@ version = "0.15.2+0"
 [[deps.libblastrampoline_jll]]
 deps = ["Artifacts", "Libdl"]
 uuid = "8e850b90-86db-534c-a0d3-1478176c7d93"
-version = "5.11.0+0"
+version = "5.8.0+1"
 
 [[deps.libdecor_jll]]
 deps = ["Artifacts", "Dbus_jll", "JLLWrappers", "Libdl", "Libglvnd_jll", "Pango_jll", "Wayland_jll", "xkbcommon_jll"]

benchmarks/DynamicalODE/Quadrupole_boson_Hamiltonian_energy_conservation_benchmark.jmd

@@ -29,23 +29,6 @@ function iip_dq(dq, p, q, params, t)
     dq[2] = A * p[2]
 end
 
-function manifold_jacobian!(J, u)
-    p₁, p₂ = u[1], u[2]
-    q₁, q₂ = u[3], u[4]
-    
-    q₁_sq = q₁^2
-    q₂_sq = q₂^2
-    q_sum_sq = q₁_sq + q₂_sq
-    
-    J[1,1] = A * p₁
-    J[1,2] = A * p₂
-    J[1,3] = A * q₁ + (B/√2) * (3q₂_sq - 3q₁_sq) + D * q₁ * q_sum_sq
-    J[1,4] = A * q₂ + (B/√2) * 6q₁q₂ + D * q₂ * q_sum_sq
-    J[2:4,:] .= 0
-    
-    nothing
-end
-
 function iip_dp(dp, p, q, params, t)
     dp[1] = -A * q[1] - 3 * B / sqrt(2) * (q[2]^2 - q[1]^2) - D * q[1] * (q[1]^2 + q[2]^2)
     dp[2] = -q[2] * (A + 3 * sqrt(2) * B * q[1] + D * (q[1]^2 + q[2]^2))
@@ -84,7 +67,7 @@ function g(resid, u, p)
 end
 
 const E = H(iip_p0, iip_q0, nothing)
-const cb = ManifoldProjection(g;  nlsolve = (f,u0) -> (res = nlsolve(f, u0; ftol=1e-13); res.zero),  manifold_jacobian=manifold_jacobian!)
+const cb = ManifoldProjection(g, nlopts=Dict(:ftol=>1e-13));
 ```
 
 For the comparison we will use the following function