-
-
Notifications
You must be signed in to change notification settings - Fork 56
Open
Labels
bugSomething isn't workingSomething isn't working
Description
Hi,
Following this discussion, here is the description of a bug I seem to have encountered. I'd be happy to help if this is a good first issue!
Edit: the error only appears with v4.12.0, everything works with v4.11.0.
Describe the bug π
Specifying a trace level different than TraceMinimal() (TraceAll() or TraceWithJacobianConditionNumber()) yields an error.
Expected behavior
No error should be raised as these arguments are valid according to the diagnostics API.
Minimal Reproducible Example π
using NonlinearSolve
function nlfunc(resid, u0, p)
resid[1] = u0[1] * u0[1] - p
resid[2] = u0[2] * u0[2] - p
resid[3] = u0[3] * u0[3] - p
nothing
end
prob = NonlinearProblem(nlfunc, [1.0, 3.0, 5.0], 2.0)
solve(prob; show_trace = Val(true), trace_level = TraceAll())Error & Stacktrace
Algorithm: NewtonRaphson(
descent = NewtonDescent(),
autodiff = AutoForwardDiff(),
vjp_autodiff = AutoFiniteDiff(
fdtype = Val{:forward}(),
fdjtype = Val{:forward}(),
fdhtype = Val{:hcentral}(),
dir = true
),
jvp_autodiff = AutoForwardDiff(),
concrete_jac = Val{false}()
)
---- ------------- ----------- -------
Iter f(u) inf-norm Step 2-norm cond(J)
---- ------------- ----------- -------
0 2.30000000e+01 0.00000000e+00 5.00000000e+00
ERROR: MethodError: Cannot `convert` an object of type
Nothing to an object of type
AbstractArray
The function `convert` exists, but no method is defined for this combination of argument types.
Closest candidates are:
convert(::Type{T}, ::T) where T
@ Base Base.jl:126
convert(::Type{T}, ::T) where T<:AbstractArray
@ Base abstractarray.jl:16
convert(::Type{T}, ::LinearAlgebra.Factorization) where T<:AbstractArray
@ LinearAlgebra ~/.julia/juliaup/julia-1.11.4+0.x64.linux.gnu/share/julia/stdlib/v1.11/LinearAlgebra/src/factorization.jl:104
...
Stacktrace:
[1] #update_trace!#56
@ ~/.julia/packages/NonlinearSolveBase/2E600/src/tracing.jl:207 [inlined]
[2] update_trace!
@ ~/.julia/packages/NonlinearSolveBase/2E600/src/tracing.jl:186 [inlined]
[3] #update_trace!#57
@ ~/.julia/packages/NonlinearSolveBase/2E600/src/tracing.jl:227 [inlined]
[4] update_trace!
@ ~/.julia/packages/NonlinearSolveBase/2E600/src/tracing.jl:220 [inlined]
[5] step!(cache::NonlinearSolveFirstOrder.GeneralizedFirstOrderAlgorithmCache{β¦}; recompute_jacobian::Nothing)
@ NonlinearSolveFirstOrder ~/.julia/packages/NonlinearSolveFirstOrder/gwWsK/src/solve.jl:333
[6] step!
@ ~/.julia/packages/NonlinearSolveFirstOrder/gwWsK/src/solve.jl:232 [inlined]
[7] #step!#163
@ ~/.julia/packages/NonlinearSolveBase/2E600/src/solve.jl:500 [inlined]
[8] step!
@ ~/.julia/packages/NonlinearSolveBase/2E600/src/solve.jl:494 [inlined]
[9] solve!(cache::NonlinearSolveFirstOrder.GeneralizedFirstOrderAlgorithmCache{β¦})
@ NonlinearSolveBase ~/.julia/packages/NonlinearSolveBase/2E600/src/solve.jl:240
[10] __solve(::NonlinearProblem{β¦}, ::GeneralizedFirstOrderAlgorithm{β¦}; kwargs::@Kwargs{β¦})
@ NonlinearSolveBase ~/.julia/packages/NonlinearSolveBase/2E600/src/solve.jl:228
[11] __solve
@ ~/.julia/packages/NonlinearSolveBase/2E600/src/solve.jl:223 [inlined]
[12] macro expansion
@ ~/.julia/packages/NonlinearSolveBase/2E600/src/solve.jl:419 [inlined]
[13] __generated_polysolve(::NonlinearProblem{β¦}, ::NonlinearSolvePolyAlgorithm{β¦}; stats::SciMLBase.NLStats, alias_u0::Bool, verbose::Bool, initializealg::NonlinearSolveBase.NonlinearSolveDefaultInit, kwargs::@Kwargs{β¦})
@ NonlinearSolveBase ~/.julia/packages/NonlinearSolveBase/2E600/src/solve.jl:376
[14] __generated_polysolve
@ ~/.julia/packages/NonlinearSolveBase/2E600/src/solve.jl:376 [inlined]
[15] #__solve#152
@ ~/.julia/packages/NonlinearSolveBase/2E600/src/solve.jl:349 [inlined]
[16] __solve
@ ~/.julia/packages/NonlinearSolveBase/2E600/src/solve.jl:346 [inlined]
[17] #__solve#18
@ ~/.julia/packages/NonlinearSolve/JTN4J/src/default.jl:27 [inlined]
[18] __solve
@ ~/.julia/packages/NonlinearSolve/JTN4J/src/default.jl:24 [inlined]
[19] #__solve#153
@ ~/.julia/packages/NonlinearSolveBase/2E600/src/solve.jl:360 [inlined]
[20] __solve
@ ~/.julia/packages/NonlinearSolveBase/2E600/src/solve.jl:352 [inlined]
[21] #solve_call#139
@ ~/.julia/packages/NonlinearSolveBase/2E600/src/solve.jl:154 [inlined]
[22] solve_call
@ ~/.julia/packages/NonlinearSolveBase/2E600/src/solve.jl:120 [inlined]
[23] solve_up(::NonlinearProblem{β¦}, ::Nothing, ::Vector{β¦}, ::Float64; originator::SciMLBase.ChainRulesOriginator, kwargs::@Kwargs{β¦})
@ NonlinearSolveBase ~/.julia/packages/NonlinearSolveBase/2E600/src/solve.jl:108
[24] solve_up
@ ~/.julia/packages/NonlinearSolveBase/2E600/src/solve.jl:101 [inlined]
[25] solve(::NonlinearProblem{β¦}; sensealg::Nothing, u0::Nothing, p::Nothing, wrap::Val{β¦}, kwargs::@Kwargs{β¦})
@ NonlinearSolveBase ~/.julia/packages/NonlinearSolveBase/2E600/src/solve.jl:81
[26] top-level scope
@ ~/Documents/code/julia/tests/bug_nonlinearsolve/bug.jl:12
Some type information was truncated. Use `show(err)` to see complete types.Environment (please complete the following information):
- Output of
using Pkg; Pkg.status()
[8913a72c] NonlinearSolve v4.12.0- Output of
using Pkg; Pkg.status(; mode = PKGMODE_MANIFEST)
[47edcb42] ADTypes v1.18.0
[7d9f7c33] Accessors v0.1.42
[79e6a3ab] Adapt v4.4.0
[4fba245c] ArrayInterface v7.20.0
[4c555306] ArrayLayouts v1.12.0
[70df07ce] BracketingNonlinearSolve v1.5.0
[d360d2e6] ChainRulesCore v1.26.0
[38540f10] CommonSolve v0.2.4
[bbf7d656] CommonSubexpressions v0.3.1
[34da2185] Compat v4.18.1
[a33af91c] CompositionsBase v0.1.2
[2569d6c7] ConcreteStructs v0.2.3
[187b0558] ConstructionBase v1.6.0
[163ba53b] DiffResults v1.1.0
[b552c78f] DiffRules v1.15.1
[a0c0ee7d] DifferentiationInterface v0.7.8
[ffbed154] DocStringExtensions v0.9.5
[4e289a0a] EnumX v1.0.5
[f151be2c] EnzymeCore v0.8.14
[e2ba6199] ExprTools v0.1.10
[55351af7] ExproniconLite v0.10.14
[9aa1b823] FastClosures v0.3.2
[1a297f60] FillArrays v1.14.0
[6a86dc24] FiniteDiff v2.28.1
[f6369f11] ForwardDiff v1.2.1
[069b7b12] FunctionWrappers v1.1.3
[77dc65aa] FunctionWrappersWrappers v0.1.3
[46192b85] GPUArraysCore v0.2.0
[3587e190] InverseFunctions v0.1.17
[92d709cd] IrrationalConstants v0.2.4
[82899510] IteratorInterfaceExtensions v1.0.0
[692b3bcd] JLLWrappers v1.7.1
[ae98c720] Jieko v0.2.1
[ba0b0d4f] Krylov v0.10.2
[5078a376] LazyArrays v2.7.0
[87fe0de2] LineSearch v0.1.4
[7ed4a6bd] LinearSolve v3.41.0
[2ab3a3ac] LogExpFunctions v0.3.29
[1914dd2f] MacroTools v0.5.16
[bb5d69b7] MaybeInplace v0.1.4
[2e0e35c7] Moshi v0.3.7
[77ba4419] NaNMath v1.1.3
[8913a72c] NonlinearSolve v4.12.0
[be0214bd] NonlinearSolveBase v2.0.0
[5959db7a] NonlinearSolveFirstOrder v1.9.0
[9a2c21bd] NonlinearSolveQuasiNewton v1.10.0
[26075421] NonlinearSolveSpectralMethods v1.5.0
[d236fae5] PreallocationTools v0.4.34
β
[aea7be01] PrecompileTools v1.2.1
[21216c6a] Preferences v1.5.0
[3cdcf5f2] RecipesBase v1.3.4
[731186ca] RecursiveArrayTools v3.37.1
[189a3867] Reexport v1.2.2
[ae029012] Requires v1.3.1
[7e49a35a] RuntimeGeneratedFunctions v0.5.15
[0bca4576] SciMLBase v2.120.0
[19f34311] SciMLJacobianOperators v0.1.11
[c0aeaf25] SciMLOperators v1.7.2
[53ae85a6] SciMLStructures v1.7.0
[efcf1570] Setfield v1.1.2
[727e6d20] SimpleNonlinearSolve v2.9.0
[276daf66] SpecialFunctions v2.6.1
[90137ffa] StaticArrays v1.9.15
[1e83bf80] StaticArraysCore v1.4.3
[10745b16] Statistics v1.11.1
[2efcf032] SymbolicIndexingInterface v0.3.45
[a759f4b9] TimerOutputs v0.5.29
[3a884ed6] UnPack v1.0.2
[1d5cc7b8] IntelOpenMP_jll v2025.2.0+0
[856f044c] MKL_jll v2025.2.0+0
[efe28fd5] OpenSpecFun_jll v0.5.6+0
[1317d2d5] oneTBB_jll v2022.0.0+0
[0dad84c5] ArgTools v1.1.2
[56f22d72] Artifacts v1.11.0
[2a0f44e3] Base64 v1.11.0
[ade2ca70] Dates v1.11.0
[8ba89e20] Distributed v1.11.0
[f43a241f] Downloads v1.6.0
[7b1f6079] FileWatching v1.11.0
[9fa8497b] Future v1.11.0
[b77e0a4c] InteractiveUtils v1.11.0
[4af54fe1] LazyArtifacts v1.11.0
[b27032c2] LibCURL v0.6.4
[76f85450] LibGit2 v1.11.0
[8f399da3] Libdl v1.11.0
[37e2e46d] LinearAlgebra v1.11.0
[56ddb016] Logging v1.11.0
[d6f4376e] Markdown v1.11.0
[ca575930] NetworkOptions v1.2.0
[44cfe95a] Pkg v1.11.0
[de0858da] Printf v1.11.0
[9a3f8284] Random v1.11.0
[ea8e919c] SHA v0.7.0
[9e88b42a] Serialization v1.11.0
[6462fe0b] Sockets v1.11.0
[2f01184e] SparseArrays v1.11.0
[fa267f1f] TOML v1.0.3
[a4e569a6] Tar v1.10.0
[cf7118a7] UUIDs v1.11.0
[4ec0a83e] Unicode v1.11.0
[e66e0078] CompilerSupportLibraries_jll v1.1.1+0
[deac9b47] LibCURL_jll v8.6.0+0
[e37daf67] LibGit2_jll v1.7.2+0
[29816b5a] LibSSH2_jll v1.11.0+1
[c8ffd9c3] MbedTLS_jll v2.28.6+0
[14a3606d] MozillaCACerts_jll v2023.12.12
[4536629a] OpenBLAS_jll v0.3.27+1
[05823500] OpenLibm_jll v0.8.1+4
[bea87d4a] SuiteSparse_jll v7.7.0+0
[83775a58] Zlib_jll v1.2.13+1
[8e850b90] libblastrampoline_jll v5.11.0+0
[8e850ede] nghttp2_jll v1.59.0+0
[3f19e933] p7zip_jll v17.4.0+2- Output of
versioninfo()
Julia Version 1.11.4
Commit 8561cc3d68d (2025-03-10 11:36 UTC)
Build Info:
Official https://julialang.org/ release
Platform Info:
OS: Linux (x86_64-linux-gnu)
CPU: 14 Γ Intel(R) Core(TM) Ultra 7 165U
WORD_SIZE: 64
LLVM: libLLVM-16.0.6 (ORCJIT, alderlake)
Threads: 1 default, 0 interactive, 1 GC (on 14 virtual cores)
Environment:
JULIA_EDITOR = code
JULIA_VSCODE_REPL = 1Metadata
Metadata
Assignees
Labels
bugSomething isn't workingSomething isn't working