Specifying maxiters and CubaSUAVE in solve causes julia to crash (segmentation fault)
Nobody has claimed this yet.
- Dominant language
- Julia
- Stars
- 245
- Forks
- 36
- Avg merge
- 4h 43m
- Merged PRs (30d)
- 15
Description
Specifying maxiters with CubaSUAVE() causes julia to crash, e.g.
using Integrals, IntegralsCuba
f(x,p) = sum(sin.(x))
prob = IntegralProblem(f,ones(2),3ones(2))
sol = solve(prob,CubaSUAVE(),reltol=1e-3,abstol=1e-3, maxiters=10)
Also I'd like to check - is maxiters meant to be a hard limit for the number of function calls? I would have expected the below to have very different results, but they give the exact same value.
using Integrals
f(x,p) = sum(sin.(x))
prob = IntegralProblem(f,ones(2),3ones(2))
sol1 = solve(prob,HCubatureJL(),reltol=1e-8,abstol=1e-8, maxiters=1).u
sol2 = solve(prob,HCubatureJL(),reltol=1e-8,abstol=1e-8, maxiters=10).u
sol1 == sol2
Contributor guide
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
Research direction
Start by running the reported CubaSUAVE reproduction with maxiters=10 and compare the HCubatureJL maxiters=1 and maxiters=10 examples. Trace the solve paths for CubaSUAVE and HCubatureJL to determine how maxiters is handled. Done means the CubaSUAVE case no longer segfaults and maxiters behavior is made consistent with the documented expectation.
Written by the indexing model from the issue text.
Assessment
- Tech stack
- julia
- Domain
- backend
- Issue type
- Bug
- Difficulty
- 4/5
- Estimated time
- 3-5 days
- Activity status
- Stale
- Clarity
- Mostly clear
- Newbie friendliness
- 35/100