Daan Huybrechs

Linear indexing into block arrays

3

This issue relates to #168, where linear indices in a BlockIndex led to a stack overflow. A question that arises is whether linear indexing should be supported and if so,...

sqrt of a BigFloat matrix

1

Any chance of getting `sqrt` to work? It seems like LinearAlgebra makes assumptions about what schur returns in [this function](https://github.com/JuliaLang/julia/blob/8cbe1297f98aa477da5d98ebfaa4027205b35440/stdlib/LinearAlgebra/src/dense.jl#L884). Example: ```julia julia> using GenericLinearAlgebra s julia> A = rand(BigFloat,...

move definition of Domain to DomainSetsCore package

13

This is an up-to-date pull request for the possible move of Domain{T} to [DomainSetsCore.jl](https://github.com/JuliaApproximation/DomainSetsCore.jl), see #169 I did not move these two lines: ```julia Base.IteratorSize(::Type{

I had made a minimal implementation of DomainSetsCore.jl, a lightweight package that defines the Domain type and a domain interface. I've updated the package this week: see https://github.com/JuliaApproximation/DomainSetsCore.jl. As things...

DomainSetsCore.jl package?

17

I believe there was a suggestion at some point to create a new package that would define a Domain type. I can't seem to find the comment. Anyway, I think...

Computing generalized eigenvalues

4

This is the PR mentioned in #126 and based on code provided to me by @thijssteel. It implements `eigval(A, B)` for generic matrices A and B using the QZ algorithm....

Generalized eigenvalue decomposition

2

The generalized eigenvalue decomposition is not implemented in GenericLinearAlgebra: ```julia julia> using GenericLinearAlgebra, DoubleFloats, LinearAlgebra julia> A = rand(Double64, 5, 5); B = rand(Double64, 5, 5); julia> eigvals(A, B) ERROR:...

Roadmap for future development of DomainSets

3

Here are a few areas I'm thinking of to improve DomainSets. I'm linking to existing issues. 1. Factor out the map code. It has no place in this package, and...

package extensions and interoperability with different domain types

1

At the start of DomainSets we identified a host of packages that did similar work (#14). There is an uncountable number of Julia packages defining a Point type! The combination...

As mentioned [here](https://github.com/JuliaMath/IntervalSets.jl/issues/183): ```julia julia> 1..2 * 2 1 .. 4 julia> 1..2 .* 2 1 .. 4 ``` We probably only want the latter and not the former, because...