Christian Schilling

I reopen here because this was tagged as a "bug". @mforets, is there some actionable point or should we change the tag?

> About the third item in your list, is it about > > https://github.com/JuliaIntervals/IntervalRootFinding.jl/blob/754c005c3354faf8a46ef803a2845378b7b4c065/test/slopes.jl#L12-L24 > ? If so, things like `interval(T(...), T(...))` should be written `interval(T, ..., ...)`. No, that...

> Mhmm so since promotion between `Number` and interval types are disallowed, the fix may be to define `f(z) = z^2 - interval(1)` (since promotion between `Compl` and `interval` seems...

We also have two functions to compute the remainder: `_exp_remainder` and `_exp_remainder_series`. Currently we always use the first one. This function could be the parameter to compute the remainder, so...

I agree. In the application we only want to increment the power symbolically, i.e., from `A^k` we want to compute `A^{k+1}`. Proposal (not tested): ```julia struct IntervalMatrixPower M::SymbolicMatrix # `SymbolicMatrix`...

> It may be convenient to subtype If we just use the result of M^k computed symbolically that is not going to help too much. I believe that knowing the...

I think the wrapper should combine #83 and #88: Every square number can and should be computed exactly (and rather efficiently). (Note that this is more precise than the symbolic...

... and use the function in the existing code base.

The problem is that `sort!` gives the wrong order due to the `-0.0`. It uses the order relation `

Sorry, should have written that. I propose to define the dimension of a hybrid system globally. Different modes can have different dynamics, constraint sets etc. But the dimension should be...