Sheehan Olver
Sheehan Olver
What you mentioned is exactly the type of application we have in mind for the package. If you give a MWE I can try to get it working
There’s no way a MWE requires knowing about spin-weighted SH. Usually for SH the operators decouple across Fourier modes so there is probably a simple example of a banded matrix....
There is nothing in the gist that I can obviously copy-and-paste so its not very helpful
You just need low order discrete eigenvalues? If so it would be easy to port over the ApproxFun code: https://github.com/JuliaApproximation/ApproxFunBase.jl/blob/master/src/eigen.jl
That could be doable but assuming you only have discrete eigenvalues Eigenvectors should be an ∞ matrix btw. For continuous spectrum it would have to return a QuasiArray. Actually I...
As long as we don’t claim to guarantee an ordering it isn’t that hard I think: just do simealtanous inverse iteration with Rayleigh shifts. Since eigenvectors have decay we would...
A PR would be greatly appreciated!
I don't have time to look at this, would you be able to submit a PR?
Fixed: ```julia julia> ApplyArray(*, ones(10,10)) (10×10 Matrix{Float64}): 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0...