Put trimmable code into package

See TrimTest.jl comment for rationale.

Author: RomeoV

test/trim/Project.toml

@@ -1,3 +1,6 @@
+name = "TrimTest"
+uuid = "7e54ada7-ece5-4046-aa01-512d530850d8"
+
 [deps]
 ADTypes = "47edcb42-4c32-4615-8424-f2b9edc5f35b"
 DiffEqBase = "2b5f629d-d688-5b77-993f-72d75c75574e"

test/trim/clean_optimization.jl

@@ -5,8 +5,6 @@ using LinearAlgebra
 using StaticArrays
 using LinearSolve
 const LS = LinearSolve
-using SciMLBase: successful_retcode
-using JET
 
 function f(u, p)
     L, U = cholesky(p.Σ)
@@ -32,6 +30,3 @@ function minimize(x)
     sol = solve(prob, alg)
     return sol
 end
-
-@test successful_retcode(minimize(1.0).retcode)
-@test_opt minimize(1.0)

test/trim/main.jl

@@ -1,11 +1,7 @@
 using TrimTest
 
 function (@main)(argv::Vector{String})::Cint
-    λ = try
-        parse(Float64, argv[2])
-    catch
-        parse(Float64, argv[1])
-    end
+    λ = parse(Float64, argv[2])
     sol = TrimTest.minimize(λ)
     println(Core.stdout, sum(sol.u))
     return 0

test/trim/src/TrimTest.jl

@@ -1,3 +1,26 @@
 module TrimTest
+#=
+Currently, trimming only works if the target code is in a package. I.e., trying to trim
+```julia
+include("trimmable_optimization.jl")
+function (@main)(argv::Vector{String})::Cint
+    minimize(1.0)
+    return 0
+end
+```
+or even
+```julia
+mod MyMod
+    include("trimmable_optimization.jl")
+end
+function (@main)(argv::Vector{String})::Cint
+    MyMod.minimize(1.0)
+    return 0
+end
+```
+segfaults `juliac`. Looking at the segfault stacktrace it seems the culprit is
+`const cache = init(...)`. Either way, we circumvent the segfault by putting
+this below code into a package definition.
+=#
 include("../trimmable_optimization.jl")
 end

test/trim/trimmable_optimization.jl

@@ -1,13 +1,12 @@
 using NonlinearSolveFirstOrder
-using ADTypes: AutoForwardDiff
 using DiffEqBase
+using ADTypes: AutoForwardDiff
 using ForwardDiff
 using LinearAlgebra
 using StaticArrays
 using LinearSolve
+import SciMLBase
 const LS = LinearSolve
-using SciMLBase: successful_retcode
-using JET
 
 function f(u, p)
     L, U = cholesky(p.Σ)
@@ -24,7 +23,6 @@ struct MyParams{T, M}
     λ::T
     Σ::M
 end
-DiffEqBase.anyeltypedual(::MyParams) = Any
 
 const autodiff = AutoForwardDiff(; chunksize = 1)
 const alg = TrustRegion(; autodiff, linsolve = LS.CholeskyFactorization())
@@ -37,6 +35,3 @@ function minimize(x)
     solve!(cache)
     return cache
 end
-
-@test successful_retcode(minimize(1.0).retcode)
-@test_opt minimize(1.0)