Daan Huybrechs

Okay, I'll try to look into it, if it doesn't become too comprehensive. I'm not very familiar with BlockArrays.jl yet, though for another package the functionality is just what I...

That was doable, I updated the pull request #168.

Big effort, no doubt, but this would be great!

For reference, see also https://github.com/AndrewGibbs/NSDpackage (this is a Matlab package and not related to moment-based quadrature)

Neat. I wasn't aware that `eig` on `BigFloat` tridiagonal matrices would work, since it doesn't for full matrices. Modifying the expansions would indeed be simple and, especially for Laguerre, it...

Perhaps I'm not using the right version then. I receive a `MethodError` related to `eigfact` for a dense `BigFloat`. It does work for tridiagonal matrices. I'll look into it.

Still, the fact that it works for tridiagonal matrices is very useful and is good to know of.

I think the ship has sailed, up to some details. The choice in this package is that an interval is defined by the membership function `a eltype(A) Int64 julia> 2.0...

Perhaps we could have `1..2` turn into `1.0..2.0`, because the notation implies the intention of a continuous set, while leaving `Interval(1,2)` as it is?

> > Perhaps we could have 1..2 turn into 1.0..2.0, because the notation implies the intention of a continuous set, while leaving Interval(1,2) as it is? > > We currently...