Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Cholesky instead of Lanczos #131

Open
wants to merge 6 commits into
base: master
Choose a base branch
from
Open

Cholesky instead of Lanczos #131

wants to merge 6 commits into from

Conversation

dlfivefifty
Copy link
Member

Fixes #130.

julia> n = 100; P = SemiclassicalJacobi(2, 1/2, 0, 1/2); @time P[0.1,n]
  0.064808 seconds (18.94 k allocations: 6.870 MiB)
-1.9278461993824534

julia> n = 200; P = SemiclassicalJacobi(2, 1/2, 0, 1/2); @time P[0.1,n]
  0.149148 seconds (38.25 k allocations: 14.203 MiB, 6.32% gc time)
-0.4280837544670504

julia> n = 400; P = SemiclassicalJacobi(2, 1/2, 0, 1/2); @time P[0.1,n]
  0.247942 seconds (76.85 k allocations: 29.288 MiB)
0.6259343524965154

julia> n = 800; P = SemiclassicalJacobi(2, 1/2, 0, 1/2); @time P[0.1,n]
  0.548767 seconds (180.02 k allocations: 62.887 MiB, 2.32% gc time)
0.9986605877062638

Speed is still very slow but at least its O(n).

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

Package should be O(n) not O(n^2)
1 participant