# Optimization program which previously worked is no longer working

**URL:** https://discourse.julialang.org/t/optimization-program-which-previously-worked-is-no-longer-working/118900
**Category:** General Usage
**Tags:** question, diffeq, optimization, differentialequation
**Created:** [September 1, 2024, 6:19am UTC](https://discourse.julialang.org/t/optimization-program-which-previously-worked-is-no-longer-working/118900 "2024-09-01T06:19:31Z")
**Posts on this page:** 13
**Page:** 1

<div class="post-metadata">

### Author: ![tighem](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/tighem/32/210486_2.png) [@tighem](https://discourse.julialang.org/u/tighem)
#### Post date: [September 1, 2024, 6:19am UTC](https://discourse.julialang.org/t/optimization-program-which-previously-worked-is-no-longer-working/118900/1 "2024-09-01T06:19:32Z")

</div>

I am running julia version 1.10.5, about a month ago I had written and successfully run a piece of optimization code which uses LBFGS to optimize a system of 218 differential equations dependent on 12 parameters. Returning to this piece of code I am no longer able to run it (however it still works when run on smaller inputs). Since then my package environment _has_ changed so perhaps this is the issue? I am unsure. IWhat are the possible causes and solutions to this issue? Now when I attempt to run the optimization code on the same input that worked a month ago I get the following error:

`[545371] signal (15): Terminated  
in expression starting at /home/tighemcasey/StickySquaresProject/Optimize3x3/BestDeltasRun.jl:110  
jl\_gc\_szclass at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/julia\_internal.h:442 [inlined]  
jl\_gc\_alloc\_ at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/julia\_internal.h:473 [inlined]  
jl\_gc\_alloc at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gc.c:3583  
ijl\_new\_bits at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/datatype.c:970  
ijl\_get\_nth\_field at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/datatype.c:1598  
do\_apply at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/builtins.c:707

- at ./operators.jl:587 [inlined]  
f1! at /home/tighemcasey/StickySquaresProject/Optimize3x3/Equations.jl:177  
Void at /home/tighemcasey/.julia/packages/SciMLBase/HkaWE/src/utils.jl:482 [inlined]  
CallWrapper at /home/tighemcasey/.julia/packages/FunctionWrappers/Q5cBx/src/FunctionWrappers.jl:65  
unknown function (ip: 0x7f1e5ed96404)  
macro expansion at /home/tighemcasey/.julia/packages/FunctionWrappers/Q5cBx/src/FunctionWrappers.jl:137 [inlined]  
do\_ccall at /home/tighemcasey/.julia/packages/FunctionWrappers/Q5cBx/src/FunctionWrappers.jl:125 [inlined]  
FunctionWrapper at /home/tighemcasey/.julia/packages/FunctionWrappers/Q5cBx/src/FunctionWrappers.jl:144 [inlined]  
\_call at /home/tighemcasey/.julia/packages/FunctionWrappersWrappers/9XR0m/src/FunctionWrappersWrappers.jl:12 [inlined]  
FunctionWrappersWrapper at /home/tighemcasey/.julia/packages/FunctionWrappersWrappers/9XR0m/src/FunctionWrappersWrappers.jl:10 [inlined]  
ODEFunction at /home/tighemcasey/.julia/packages/SciMLBase/HkaWE/src/scimlfunctions.jl:2335 [inlined]  
perform\_step! at /home/tighemcasey/.julia/packages/OrdinaryDiffEqTsit5/DHYtz/src/tsit\_perform\_step.jl:204  
perform\_step! at /home/tighemcasey/.julia/packages/OrdinaryDiffEqTsit5/DHYtz/src/tsit\_perform\_step.jl:181 [inlined]  
solve! at /home/tighemcasey/.julia/packages/OrdinaryDiffEqCore/4A2vD/src/solve.jl:551  
#\_\_solve#61 at /home/tighemcasey/.julia/packages/OrdinaryDiffEqCore/4A2vD/src/solve.jl:7 [inlined]  
\_\_solve at /home/tighemcasey/.julia/packages/OrdinaryDiffEqCore/4A2vD/src/solve.jl:1 [inlined]  
#solve\_call#44 at /home/tighemcasey/.julia/packages/DiffEqBase/sCsah/src/solve.jl:612 [inlined]  
solve\_call at /home/tighemcasey/.julia/packages/DiffEqBase/sCsah/src/solve.jl:569 [inlined]  
#solve\_up#53 at /home/tighemcasey/.julia/packages/DiffEqBase/sCsah/src/solve.jl:1080 [inlined]  
solve\_up at /home/tighemcasey/.julia/packages/DiffEqBase/sCsah/src/solve.jl:1066 [inlined]  
#solve#51 at /home/tighemcasey/.julia/packages/DiffEqBase/sCsah/src/solve.jl:1003 [inlined]  
solve at /home/tighemcasey/.julia/packages/DiffEqBase/sCsah/src/solve.jl:993  
unknown function (ip: 0x7f1e5edad93f)  
\_jl\_invoke at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gf.c:2895 [inlined]  
ijl\_apply\_generic at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gf.c:3077  
loss at /home/tighemcasey/StickySquaresProject/Optimize3x3/BestDeltasRun.jl:38  
#1 at /home/tighemcasey/StickySquaresProject/Optimize3x3/BestDeltasRun.jl:80 [inlined]  
#37 at /home/tighemcasey/.julia/packages/OptimizationBase/KIIy3/ext/OptimizationForwardDiffExt.jl:98 [inlined]  
#39 at /home/tighemcasey/.julia/packages/OptimizationBase/KIIy3/ext/OptimizationForwardDiffExt.jl:102 [inlined]  
vector\_mode\_dual\_eval! at /home/tighemcasey/.julia/packages/ForwardDiff/PcZ48/src/apiutils.jl:24 [inlined]  
vector\_mode\_gradient! at /home/tighemcasey/.julia/packages/ForwardDiff/PcZ48/src/gradient.jl:96  
gradient! at /home/tighemcasey/.julia/packages/ForwardDiff/PcZ48/src/gradient.jl:37 [inlined]  
#38 at /home/tighemcasey/.julia/packages/OptimizationBase/KIIy3/ext/OptimizationForwardDiffExt.jl:102  
unknown function (ip: 0x7f1e5ef88369)  
\_jl\_invoke at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gf.c:2895 [inlined]  
ijl\_apply\_generic at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gf.c:3077  
jl\_apply at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/julia.h:1982 [inlined]  
do\_apply at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/builtins.c:768  
#19 at /home/tighemcasey/.julia/packages/OptimizationOptimJL/hDX5k/src/OptimizationOptimJL.jl:297  
unknown function (ip: 0x7f1e5ef8436c)  
\_jl\_invoke at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gf.c:2895 [inlined]  
ijl\_apply\_generic at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gf.c:3077  
value\_gradient!! at /home/tighemcasey/.julia/packages/NLSolversBase/kavn7/src/interface.jl:82  
value\_gradient! at /home/tighemcasey/.julia/packages/NLSolversBase/kavn7/src/interface.jl:69  
value\_gradient! at /home/tighemcasey/.julia/packages/Optim/ZhuZN/src/multivariate/solvers/constrained/fminbox.jl:94  
unknown function (ip: 0x7f1e5eef77f0)  
\_jl\_invoke at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gf.c:2895 [inlined]  
ijl\_apply\_generic at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gf.c:3077  
value\_gradient! at /home/tighemcasey/.julia/packages/Optim/ZhuZN/src/Manifolds.jl:50  
ϕdϕ at /home/tighemcasey/.julia/packages/LineSearches/jgnxK/src/LineSearches.jl:83  
HagerZhang at /home/tighemcasey/.julia/packages/LineSearches/jgnxK/src/hagerzhang.jl:233  
HagerZhang at /home/tighemcasey/.julia/packages/LineSearches/jgnxK/src/hagerzhang.jl:102 [inlined]  
perform\_linesearch! at /home/tighemcasey/.julia/packages/Optim/ZhuZN/src/utilities/perform\_linesearch.jl:58  
update\_state! at /home/tighemcasey/.julia/packages/Optim/ZhuZN/src/multivariate/solvers/first\_order/l\_bfgs.jl:204  
optimize at /home/tighemcasey/.julia/packages/Optim/ZhuZN/src/multivariate/optimize/optimize.jl:54  
optimize at /home/tighemcasey/.julia/packages/Optim/ZhuZN/src/multivariate/solvers/constrained/fminbox.jl:390  
unknown function (ip: 0x7f1e5ef83d6d)  
\_jl\_invoke at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gf.c:2895 [inlined]  
ijl\_apply\_generic at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gf.c:3077  
\_\_solve at /home/tighemcasey/.julia/packages/OptimizationOptimJL/hDX5k/src/OptimizationOptimJL.jl:321  
solve! at /home/tighemcasey/.julia/packages/SciMLBase/HkaWE/src/solve.jl:188  
unknown function (ip: 0x7f1e5eed9a65)  
\_jl\_invoke at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gf.c:2895 [inlined]  
ijl\_apply\_generic at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gf.c:3077  
#solve#647 at /home/tighemcasey/.julia/packages/SciMLBase/HkaWE/src/solve.jl:96  
solve at /home/tighemcasey/.julia/packages/SciMLBase/HkaWE/src/solve.jl:93  
unknown function (ip: 0x7f1e5f06e5f0)  
\_jl\_invoke at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gf.c:2895 [inlined]  
ijl\_apply\_generic at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gf.c:3077  
yield\_vs\_Tmax at /home/tighemcasey/StickySquaresProject/Optimize3x3/BestDeltasRun.jl:83  
\_jl\_invoke at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gf.c:2895 [inlined]  
ijl\_apply\_generic at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gf.c:3077  
jl\_apply at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/julia.h:1982 [inlined]  
do\_call at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/interpreter.c:126  
eval\_value at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/interpreter.c:223  
eval\_stmt\_value at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/interpreter.c:174 [inlined]  
eval\_body at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/interpreter.c:617  
jl\_interpret\_toplevel\_thunk at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/interpreter.c:775  
jl\_toplevel\_eval\_flex at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/toplevel.c:934  
jl\_toplevel\_eval\_flex at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/toplevel.c:877  
ijl\_toplevel\_eval\_in at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/toplevel.c:985  
eval at ./boot.jl:385 [inlined]  
include\_string at ./loading.jl:2076  
\_jl\_invoke at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gf.c:2895 [inlined]  
ijl\_apply\_generic at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gf.c:3077  
\_include at ./loading.jl:2136  
include at ./Base.jl:495  
jfptr\_include\_46447.1 at /snap/julia/104/lib/julia/sys.so (unknown line)  
\_jl\_invoke at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gf.c:2895 [inlined]  
ijl\_apply\_generic at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gf.c:3077  
exec\_options at ./client.jl:318  
\_start at ./client.jl:552  
jfptr\_\_start\_82798.1 at /snap/julia/104/lib/julia/sys.so (unknown line)  
\_jl\_invoke at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gf.c:2895 [inlined]  
ijl\_apply\_generic at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gf.c:3077  
jl\_apply at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/julia.h:1982 [inlined]  
true\_main at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/jlapi.c:582  
jl\_repl\_entrypoint at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/jlapi.c:731  
main at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/cli/loader\_exe.c:58  
unknown function (ip: 0x7f1f8e031d8f)  
\_\_libc\_start\_main at /lib/x86\_64-linux-gnu/libc.so.6 (unknown line)  
unknown function (ip: 0x4010b8)  
unknown function (ip: (nil))  
Allocations: 46790740815 (Pool: 46438082700; Big: 352658115); GC: 59429`

---

<div class="post-metadata">

### Author: ![ufechner7](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/ufechner7/32/51363_2.png) [@ufechner7](https://discourse.julialang.org/u/ufechner7)
#### Post date: [September 1, 2024, 7:37am UTC](https://discourse.julialang.org/t/optimization-program-which-previously-worked-is-no-longer-working/118900/2 "2024-09-01T07:37:02Z")

</div>

> [@tighem](#):
>
> Since then my package environment _has_ changed so perhaps this is the issue?

Always make a backup copy of `Manifest.toml` when your code works to be able to reproduce the working environment if needed.

And it is hard to debug the exception you got without seeing your code and the output of `using Pkg; Pkg.status()`

---

<div class="post-metadata">

### Author: ![ChrisRackauckas](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/chrisrackauckas/32/77_2.png) [@ChrisRackauckas](https://discourse.julialang.org/u/ChrisRackauckas)
#### Post date: [September 1, 2024, 3:00pm UTC](https://discourse.julialang.org/t/optimization-program-which-previously-worked-is-no-longer-working/118900/3 "2024-09-01T15:00:53Z")

</div>

This would need reproducible code.

> [@tighem](#):
>
> `[545371] signal (15): Terminated  
> in expression starting at /home/tighemcasey/StickySquaresProject/Optimize3x3/BestDeltasRun.jl:110  
> jl\_gc\_szclass at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/julia\_internal.h:442 [inlined]  
> jl\_gc\_alloc\_ at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/julia\_internal.h:473 [inlined]  
> jl\_gc\_alloc at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/gc.c:3583  
> ijl\_new\_bits at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/datatype.c:970  
> ijl\_get\_nth\_field at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/datatype.c:1598  
> do\_apply at /cache/build/builder-amdci4-4/julialang/julia-release-1-dot-10/src/builtins.c:707
> 
> - at ./operators.jl:587 [inlined]  
> f1! at /home/tighemcasey/StickySquaresProject/Optimize3x3/Equations.jl:177  
> Void at /home/tighemcasey/.julia/packages/SciMLBase/HkaWE/src/utils.jl:482 [inlined]  
> CallWrapper at /home/tighemcasey/.juli

Are you using Julia v1.10.5? There’s GC fixes that are helpful. If you are on v1.10.5, are you using `@inbounds` anywhere in your code? What’s at that line of code?

---

<div class="post-metadata">

### Author: ![tighem](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/tighem/32/210486_2.png) [@tighem](https://discourse.julialang.org/u/tighem)
#### Post date: [September 1, 2024, 7:53pm UTC](https://discourse.julialang.org/t/optimization-program-which-previously-worked-is-no-longer-working/118900/4 "2024-09-01T19:53:09Z")

</div>

Here is the output of Pkg.installed():

```julia
"CSV" => v"0.10.14"
  "DifferentialEquations" => v"7.13.0"
  "DataFrames" => v"1.6.1"
  "Optimization" => v"3.27.0"
  "SciMLSensitivity" => v"7.65.0"
  "ChangePrecision" => v"1.1.0"
  "OptimizationOptimJL" => v"0.3.2"

```

Here is the optimization code, the file Constants.txt simply contains the numerical constants, f1! assigns the value of dC as a non-linear function of deltas and C.

````using

#=
THIS FILE PRODUCES OPTIMAL BOND ENERGIES AND YIELDS AT THOSE BOND ENERGIES AT A GIVEN SET OF TIMESTEPS,
OUTPUT IS TO FILES 'YieldvsTime.txt' AND 'BindingEnergyvsTime.txt', 
NOTE THAT THE OUTPUT OF BINDING ENERGIES IS IN LEXICOGRAPHICAL ORDER
=#

# ------------------------------------ I/O and Inclusions ------------------------------------ #

# Setup Files
include("Equations.jl")
constantspath = joinpath(@ __DIR__ ,"Constants.txt")

# Output Files
Yieldpath = joinpath(@ __DIR__ ,"YieldvsTime.txt")
Bindingpath = joinpath(@ __DIR__ ,"BindingEnergyvsTime.txt")

#Read Necessary Parameters from 'Constants.txt'
params = readdlm(constantspath,skipstart=2)
concentration = params[4,1]::Float64
numblocks = params[7,1]::Int64
numEdges = params[2,1]::Int64
numVertices = params[3,1]::Int64
C = Vector{Float64}(params[5,1:numblocks])
dC = Vector{Float64}(params[6,1:numblocks])

# Write output to file
function writeoutput(Yieldpath, Bindingpath, concvec, deltaarray)
    #Write Header
    maxtime = tvec[end]::Float64
    deltastring = ""
    for k in 1:numEdges
        deltastring = deltastring*"-delta$k: "
    end
    yieldheader = "density = $concentration ; Initial deltas = $deltas ; \n Initial concentrations = $C ;\
     max time = $maxtime \n final C(Target Structure): maxtime: \n"
    bindingheader = "density = $concentration ; Initial deltas = $deltas ; \n Initial concentrations = $C ;\
     max time = $maxtime \n"*deltastring*" maxtime: \n"
    write(Yieldpath,yieldheader)
    write(Bindingpath, bindingheader)

    #Write Data
    CSV.write(Yieldpath,Tables.table(concvec),append=true)
    CSV.write(Bindingpath,Tables.table(deltaarray),append=true)
end

# ------------------------------------ Input Parameters ------------------------------------ #

tvec = [1000.0,10000.0,100000.0] # Time Steps Optimized At
deltas = [-0.3546379129258437, -1.0306384019336112, -2.765023007323766, -2.5384636391235658, -1.816160723724639, -6.0154802607957905, -4.734737685358027, -5.7959495615945364, -0.5149777996841498, -3.536189921414051, -9.460076287814163, -2.9729010182065396]
tspan = (0.0,100.0) # Time span for initial solve
maxenergy = -15.0 # Boundary condition on bonding energies

#Parameters for ODE solver and optimization
adtype = Optimization.AutoForwardDiff()
maxoptiterations = 200

#define ode problem
problem = ODEProblem(f1!,C,tspan,deltas)
sol = solve(problem,save_on = false)

# ------------------------------------ Optimization ------------------------------------ #

#Loss function
function loss(newp,t)
    newproblem = remake(problem,p = newp, tspan = (0.0,t))
    sol = solve(newproblem, Tsit5(), save_on = false)
    if Sys.free_memory() / 2^30 < 15 #Possibly helps with memory issues
        GC.gc()
    end
    loss = -sol[end,end] #Final concentration of target structure
    return loss, sol
end

function yield_vs_Tmax(tvec,problem,adtype)
    #Initialize storage for data
    deltaarray = Array{Float64}(undef,size(tvec)[1],numEdges+1)
    concvec = Array{Float64}(undef,size(tvec)[1],2)
    i = 1 #Used for indexing
    #Optimization code for each Tmax:
    for t ∈ tvec
        optf = Optimization.OptimizationFunction((x, p) -> loss(x,t), adtype) #Define optimization loss function, in this case maximizes concentration of final structure
        optprob = Optimization.OptimizationProblem(optf, deltas, lb = ones(numEdges).*maxenergy, ub = zeros(numEdges)) #define problem and boundary conditions
        result_ode = Optimization.solve(optprob, LBFGS(), maxiters = maxoptiterations)

        for j ∈ 1:numEdges+1
            if j ≤ numEdges
                deltaarray[i,j] = -result_ode.u[j]
            else
                #deltaarray[i,:] = sort(deltaarray[i,:],rev=true) #In the case of symmetry a sorted output might make more sense
                deltaarray[i,j] = t
            end
        end
        concvec[i,1] = -result_ode.objective*numVertices/concentration
        concvec[i,2] = t
        i += 1
        println("Done optimization for t = $t")
    end
    #Write results to file
    writeoutput(Yieldpath,Bindingpath, concvec, deltaarray)
end

println("Started Parameter Optimization")
yield_vs_Tmax(tvec,problem,adtype)
println("Completed Parameter Optimization")```
````

---

<div class="post-metadata">

### Author: ![tighem](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/tighem/32/210486_2.png) [@tighem](https://discourse.julialang.org/u/tighem)
#### Post date: [September 1, 2024, 7:56pm UTC](https://discourse.julialang.org/t/optimization-program-which-previously-worked-is-no-longer-working/118900/5 "2024-09-01T19:56:15Z")

</div>

I am not using @inline anywhere.

---

<div class="post-metadata">

### Author: ![ChrisRackauckas](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/chrisrackauckas/32/77_2.png) [@ChrisRackauckas](https://discourse.julialang.org/u/ChrisRackauckas)
#### Post date: [September 1, 2024, 9:25pm UTC](https://discourse.julialang.org/t/optimization-program-which-previously-worked-is-no-longer-working/118900/6 "2024-09-01T21:25:03Z")

</div>

And Julia v1.10.5?

---

<div class="post-metadata">

### Author: ![tighem](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/tighem/32/210486_2.png) [@tighem](https://discourse.julialang.org/u/tighem)
#### Post date: [September 1, 2024, 9:39pm UTC](https://discourse.julialang.org/t/optimization-program-which-previously-worked-is-no-longer-working/118900/7 "2024-09-01T21:39:56Z")

</div>

Yes

---

<div class="post-metadata">

### Author: ![ChrisRackauckas](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/chrisrackauckas/32/77_2.png) [@ChrisRackauckas](https://discourse.julialang.org/u/ChrisRackauckas)
#### Post date: [September 1, 2024, 9:46pm UTC](https://discourse.julialang.org/t/optimization-program-which-previously-worked-is-no-longer-working/118900/8 "2024-09-01T21:46:54Z")

</div>

> [@tighem](#):
>
> `include("Equations.jl")`

and?

---

<div class="post-metadata">

### Author: ![Elrod](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/elrod/32/22461_2.png) [@Elrod](https://discourse.julialang.org/u/Elrod)
#### Post date: [September 2, 2024, 12:24am UTC](https://discourse.julialang.org/t/optimization-program-which-previously-worked-is-no-longer-working/118900/9 "2024-09-02T00:24:02Z")

</div>

`@inline` and `@inbounds` do two very different things, but I’m assuming you do neither. You can start Julia with `--check-bounds=yes` to confirm that it isn’t some other package using `@inbounds` causing the problem.

---

<div class="post-metadata">

### Author: ![tighem](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/tighem/32/210486_2.png) [@tighem](https://discourse.julialang.org/u/tighem)
#### Post date: [September 19, 2024, 12:20am UTC](https://discourse.julialang.org/t/optimization-program-which-previously-worked-is-no-longer-working/118900/10 "2024-09-19T00:20:06Z")

</div>

Here is the equations function; it is over the character limit substantially

[Equations.jl](https://discourse.julialang.org/uploads/short-url/fhUKx7gd8By4JrJBzAbjUj2WEJs.jl) (145.9 KB)

---

<div class="post-metadata">

### Author: ![tighem](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/tighem/32/210486_2.png) [@tighem](https://discourse.julialang.org/u/tighem)
#### Post date: [September 19, 2024, 12:30am UTC](https://discourse.julialang.org/t/optimization-program-which-previously-worked-is-no-longer-working/118900/12 "2024-09-19T00:30:39Z")

</div>

yes, I use neither. I will try running the program with this startup option.

---

<div class="post-metadata">

### Author: ![tighem](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/tighem/32/210486_2.png) [@tighem](https://discourse.julialang.org/u/tighem)
#### Post date: [September 22, 2024, 9:15pm UTC](https://discourse.julialang.org/t/optimization-program-which-previously-worked-is-no-longer-working/118900/13 "2024-09-22T21:15:06Z")

</div>

The program has the same error when run with these options.

---

<div class="post-metadata">

### Author: ![ChrisRackauckas](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/chrisrackauckas/32/77_2.png) [@ChrisRackauckas](https://discourse.julialang.org/u/ChrisRackauckas)
#### Post date: [September 22, 2024, 10:16pm UTC](https://discourse.julialang.org/t/optimization-program-which-previously-worked-is-no-longer-working/118900/14 "2024-09-22T22:16:41Z")

</div>

What’s the reproducer and the error?
