IntervalArithmetic.jl
IntervalArithmetic.jl copied to clipboard
JuliaIntervals
Library for validated numerics using interval arithmetic
```jl julia> x = -Inf..(-1e308) [-∞, -9.99999e+307] julia> x^3 ∅ ``` This should *not* be empty.
Currently, ```julia julia> a = (1..2) × (3..4) [1, 2] × [3, 4] julia> b = (1.5..1.9) × (3..4) [1.5, 1.90001] × [3, 4] julia> b ⊂ a false ```...
Does anyone know why this is happening? ```julia using IntervalArithmetic a = interval(Float16(1.0), Float16(3.0)) a + a ``` gives: ```Julia ERROR: MethodError: no method matching rounding_raw(::Type{Float16}) Closest candidates are: rounding_raw(::Type{BigFloat})...
```jl IntervalBox{N,T}(x...) where {N,T} = IntervalBox{N,T}(x) ```
I found these issues while fixing #466 , ``` julia> a = 1//1..2//1 [1//1, 2//1] julia> b = 0//1..2//1 [0//1, 2//1] julia> a/b ERROR: MethodError: no method matching Rational{Int64}(::Float64, ::RoundingMode{:Up})...
## `isapprox` wasn't overloaded for intervals and was using the fallback `isapprox(a::Number, b::Number)`. This PR defines the method `isapprox(a::Interval, b::Interval)` ## before ```julia julia> a = 0.3..0.7 [0.299999, 0.700001] julia>...
in the comparison below, i expected `m ≈ interval(0.295, 0.395)` to hold. ```julia julia> m [0.294999, 0.395001] julia> m.lo ≈ 0.295 true julia> m.hi ≈ 0.395 true julia> m ≈...
I'm giving a seminar in a few weeks using QR and just realised that `LinearAlgebra.reflectorApply!` has issues with if the matrix is not diagonally dominant. Here's a simple demonstration of...
Metadata
Owner
Metadata
Library for validated numerics using interval arithmetic