Skip to content

Commit

Permalink
curtail large output of SymTridiagonal
Browse files Browse the repository at this point in the history
  • Loading branch information
MikaelSlevinsky committed Mar 21, 2024
1 parent 1bfa33b commit 04d8e87
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion examples/halfrange.jl
Original file line number Diff line number Diff line change
Expand Up @@ -43,7 +43,7 @@ qvals = k-> ichebyshevtransform(q(k))

# With the symmetric Jacobi matrix for $P_n^{(-\frac{1}{2}, 0)}(y)$ and the modified plan, we may compute the modified Jacobi matrix and the corresponding roots (as eigenvalues):
XP = SymTridiagonal([-inv((4n-1)*(4n-5)) for n in 1:n+1], [4n*(2n-1)/(4n-1)/sqrt((4n-3)*(4n+1)) for n in 1:n])
XQ = FastTransforms.modified_jacobi_matrix(P, XP)
XQ = FastTransforms.modified_jacobi_matrix(P, XP);

# And we plot:
x = (chebyshevpoints(Float64, n+1, Val(1)) .+ 1 ) ./ 2
Expand Down

0 comments on commit 04d8e87

Please sign in to comment.