SpectralMeasures.jl
SpectralMeasures.jl copied to clipboard
Julia package for finding the spectral measure of structured self adjoint operators
@marcusdavidwebb any objections to this change? I think this package really treats things as (infinite) matrices, so ApproxFun's domain and range space support just gets in the way.
There are issues with the current connection coefficient matrix introducing bad conditioning: ```julia V = SymTriOperator([zeros(5); -ones(6); zeros(10); -ones(7)],zeros(4)) J = -Δ + V Λ, U = eig(J) norm(Matrix(U[1:100,1:100])) #...
This dependency is a bit annoying. I might make a SpectralMeasuresPlots.jl backage
J=Δ+SpectralMeasures.SymTriOperator(0.2*[0.,0.,0.,0.,-1.,-1.,-1.,-1.],[0.,0.])
qlIteration is not updating the shift at the moment.
PertToeplitz should be closed under inv, *, transpose
Essentially L is the conversion operator between new polynomial space and Chebyshev U. So we could plug in to the ApproxFun ODE solving. Will have to think of a good...