Mikael Slevinsky

I mean that dsbev.f is slow compared with dsyev.f (for any bandwidth). Have a look at the summary of #101 for more info. N.B. the name changes as part of...

could have `antiband` for Hankel matrices 😆

Using non-contiguous strides should be part of a wrapper, not the core code: contiguity is utilized to maximize SIMD.

Hi @mauricelanghinrichs, thanks for the bug report! In fact, 3F2 has the special method [here](https://github.com/JuliaMath/HypergeometricFunctions.jl/blob/0da7d8e44b93ccaa7165fe2309bb068d570368b9/src/weniger.jl#L324), which is a rational approximation. I think the explanation is that it's exactly hitting a...

This may sound backwards, but I would probably implement the partial derivative of the incomplete gamma function $\gamma(a,z) = \frac{z^a}{a}M(a,a+1,-z)$ with respect to $a$ with $M$ already implemented [here](https://github.com/JuliaMath/HypergeometricFunctions.jl/blob/master/src/confluent.jl) as...

I think the reason I chose this way is because if you are near the origin you don't need to call `gamma(c)` as the Maclaurin series will do. If we...

I think they can be added for now, and if need be extracted to another package, say AppellFunctions.jl. If it's a wrapper, then perhaps the Julia community-friendly version is to...

Thanks for filing the issue! Evaluating `U` near the branch cut `(-∞,0]` requires a different approach, such as Kummer's transformation #61.

Thank you for the bug report! To be brief, I do know how to fix it: it appears to be a pole/pole situation. I will try to get a fix...

I'm confused by their equation (3.1). Given a diagonal matrix X = diag(x_1,...x_m), isn't exp(tr(X)) a scalar?