JuliaNLSolvers / NLsolve.jl

Julia solvers for systems of nonlinear equations and mixed complementarity problems
Other
329 stars 66 forks source link

StackOverflowError #277

Open ccmejia opened 2 years ago

ccmejia commented 2 years ago

Hi there, Today I have installed the 1.7.0 Julia version and I'm experiencing some issues with my code (which works perfectly in 1.6.2 version). The problem appears when the size of the system of equations goes above some threshold... Here is the output:

Stacktrace: [1] dogleg!(p::Vector{Float64}, p_c::Vector{Float64}, p_i::Vector{Float64}, r::Vector{Float64}, d::Vector{Float64}, J::Matrix{Float64}, delta::Float64) @ NLsolve ~/.julia/packages/NLsolve/gJL1I/src/solvers/trust_region.jl:61 [2] trustregion(df::OnceDifferentiable{Vector{Float64}, Matrix{Float64}, Vector{Float64}}, initial_x::Vector{Float64}, xtol::Float64, ftol::Float64, iterations::Int64, store_trace::Bool, show_trace::Bool, extended_trace::Bool, factor::Float64, autoscale::Bool, cache::NLsolve.NewtonTrustRegionCache{Vector{Float64}}) @ NLsolve ~/.julia/packages/NLsolve/gJL1I/src/solvers/trust_region.jl:170 [3] trust_region (repeats 2 times) @ ~/.julia/packages/NLsolve/gJL1I/src/solvers/trust_region.jl:235 [inlined] [4] nlsolve(df::OnceDifferentiable{Vector{Float64}, Matrix{Float64}, Vector{Float64}}, initial_x::Vector{Float64}; method::Symbol, xtol::Float64, ftol::Float64, iterations::Int64, store_trace::Bool, show_trace::Bool, extended_trace::Bool, linesearch::Static, linsolve::NLsolve.var"#27#29", factor::Float64, autoscale::Bool, m::Int64, beta::Int64, aa_start::Int64, droptol::Float64) @ NLsolve ~/.julia/packages/NLsolve/gJL1I/src/nlsolve/nlsolve.jl:26 [5] nlsolve(f::Function, initial_x::Vector{Float64}; method::Symbol, autodiff::Symbol, inplace::Bool, kwargs::Base.Pairs{Symbol, Real, Tuple{Symbol, Symbol, Symbol}, NamedTuple{(:ftol, :iterations, :show_trace), Tuple{Float64, Int64, Bool}}}) @ NLsolve ~/.julia/packages/NLsolve/gJL1I/src/nlsolve/nlsolve.jl:52 [6] modosLρ(Lini::Int64, Lend::Int64, ρ0ini::Float64, ρ0end::Float64, γ::Int64, λ::Int64, C0::Int64, cs_ini::String, cs_lat::String, cs_pot::String) @ Main ./In[8]:10 [7] top-level scope @ ./timing.jl:220 [inlined] [8] top-level scope @ ./In[9]:0 [9] eval @ ./boot.jl:373 [inlined] [10] include_string(mapexpr::typeof(REPL.softscope), mod::Module, code::String, filename::String) @ Base ./loading.jl:1196

Any idea about was going on here?

pkofod commented 2 years ago

I think you forgot the actual error message. Even better, provide the example! :)

ccmejia commented 2 years ago

Oh, I forgot to bring news. 1.7.1 Julia version fixed these kind of overflow problems. Thank you!