David P. Sanders
David P. Sanders
```jl julia> 3.5 ∈ interval(3, 4) false ``` This should either work or give me an error message telling me that this is not supported for intervals and telling me...
Currently the user has no way of knowing which domain a function is defined over. We should define e.g. ```jl domain(sqrt) = 0..Inf ``` or `domain_interval(sqrt)`.
Here's a possible implementation for piecewise functions: ```jl using IntervalArithmetic using IntervalArithmetic.Symbols # to use the constructor `..` # represents a piecewise 1D function struct Piecewise{T x -> x +...
Rationals do not behave well inside `CuArray`s. Not sure there's much to do about this, except implement a simpler Rational type: ```jl julia> v = CuArray([1//3]) 1-element CuArray{Rational{Int64},1}: 1//3 julia>...
[Part of the [review](https://github.com/JuliaCon/proceedings-review/issues/148) for the JuliaCon proceedings.] As far as I can tell the use of a global store for tracking means that the package is not thread-safe. Is...
```jl julia> isapprox(1..2, 3..4) ERROR: ArgumentError: `==` is purposely not supported for intervals. See instead `isequal_interval` Stacktrace: [1] ==(x::Interval{Float64}, y::Interval{Float64}) @ IntervalArithmetic ~/Dropbox/packages/IntervalArithmetic/src/intervals/real_interface.jl:86 [2] #isapprox#809 @ ./floatfuncs.jl:308 [inlined] [3] isapprox(x::Interval{Float64},...
Currently to compute `(1..2)^(1//7)` we are converting the rational to an interval. This loses accuracy. We should instead use the corresponding MPFR function, `mpfr_rootn_si`. Unfortunately Julia does not seem to...
`relative([ C[4] / 4, D, E] )` Should compute octave and duration automatically.