IntervalRootFinding.jl

IntervalRootFinding.jl copied to clipboard

Allow keywords arguments for roots

5

I new issue found while writing the documentation. The current situation is awkward. I think we could go with the following scheme: **Positional arguments** 1. The function. 2. The search...

Kolaru

Example 5.5 from Smiley and Chun (2001) solved inconsistently

18

The example in `test/smiley_examples.jl` has 41 roots, however, only 20 are found on some setups. On other setups all 41 are identified correctly. This inconsistency was first discovered and discussed...

gwater

Remove old code

3

Remove old code: - [x] bisection.jl - [x] move out definitions of `N` used in `branch_and_prune` - [ ] krawczyk.jl - [ ] newton.jl - [ ] find_roots stuff

dpsanders

No method roots(f, v::Vector{IntervalBox})

2

I'm sure this used to work.

dpsanders

Failure to find solutions for a 10-dimensional problem

13

Problem taken from https://discourse.julialang.org/t/solvers-fail-on-nonlinear-problem-which-has-solutions/20051: ``` function f(X, p) x1,x2,x3,x4,x5,x6,x7,x8,x9,x10 = X return SVector( -p[17] * x1 + -2 * p[1] * (x1 ^ 2 / 2) + 2 * p[2]...

dpsanders

Support function returning IntervalBox and only use them internally

1

Currently, the functions passed to `roots` must return `SVector` which has several drawbacks: - Require one more explicit import from the user. - `IntervalBox` and `SVector` work in a quite...

Kolaru

Add new root status when the presence (but not uniqueness) of a solution is known

6

For Krawczyk and Newton method, the fact that a contracted interval is subset of the original interval is a sufficient condition for the presence of a solution. The proposition here...

Kolaru

Newton fails with abs

6

``` julia> f(p) = abs(1 / ( (1+p)^30 ) * 10_000 - 100) f (generic function with 1 method) julia> roots(f, -1000..1000, Newton) 2-element Array{Root{Interval{Float64}},1}: Root(∅, :unique) Root(∅, :unique) ```

dpsanders

Done in #110.

dpsanders

Add centered forms

2

Add mean value forms for scalar and vector functions, and a 3rd-order Taylor form for scalar functions.

dpsanders