@dlfivefifty I believe I am getting unexpected behaviour for the evaluation of TwoBandJacobi. Occasionally it's the negative of what I would expect. E.g.
As far as I am aware there's no (-1)^k for some k, factor in the definition. But I cannot figure out where getindex for TwoBandJacobi is defined. The closest I can see is a definition for Interlace but that's not called for evaluating TwoBandJacobi as far as I can see.
@dlfivefifty I believe I am getting unexpected behaviour for the evaluation of
TwoBandJacobi. Occasionally it's the negative of what I would expect. E.g.ρ=0.2; t=inv(1-ρ^2); x=0.4; τ = t*(1-x^2); TwoBandJacobi(ρ,1,1,0)[x,3] ≈ -SemiclassicalJacobi(t,1,1,-1/2)[τ,2]Note that I require a negative sign in front of the RHS. But for a different index
ρ=0.2; t=inv(1-ρ^2); x=0.4; τ = t*(1-x^2); TwoBandJacobi(ρ,1,1,0)[x,5] ≈ SemiclassicalJacobi(t,1,1,-1/2)[τ,3]Note the lack of the negative sign.
As far as I am aware there's no (-1)^k for some k, factor in the definition. But I cannot figure out where
getindexforTwoBandJacobiis defined. The closest I can see is a definition forInterlacebut that's not called for evaluatingTwoBandJacobias far as I can see.