Skip to content

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).

@codecov
Copy link

codecov bot commented Jan 27, 2025

Codecov Report

Attention: Patch coverage is 93.10345% with 2 lines in your changes missing coverage. Please review.

Project coverage is 94.98%. Comparing base (b3e4cb3) to head (aa78cc6).
Report is 1 commits behind head on master.

Files with missing lines Patch % Lines
src/SemiclassicalOrthogonalPolynomials.jl 93.10% 2 Missing ⚠️
Additional details and impacted files
@@           Coverage Diff           @@
##           master     #131   +/-   ##
=======================================
  Coverage   94.97%   94.98%           
=======================================
  Files           3        3           
  Lines         657      658    +1     
=======================================
+ Hits          624      625    +1     
  Misses         33       33           

☔ View full report in Codecov by Sentry.
📢 Have feedback on the report? Share it here.

@dlfivefifty dlfivefifty merged commit e06d88a into master Jan 28, 2025
10 of 11 checks passed
@dlfivefifty dlfivefifty deleted the dl/nolanczos branch January 28, 2025 16:41
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)

2 participants