search

MikaelSlevinsky / FastTransforms

:bullettrain_front: Fast orthogonal polynomial transforms :surfer:
MIT License
58 stars 9 forks source link

Fast jac2jac with integer steps #78

Open dlfivefifty opened 11 months ago

dlfivefifty commented 11 months ago 
julia> n = 1000; x = randn(n); @time FastTransforms.lib_jac2jac(x, 0, 0, 1, 1);
  0.021386 seconds (6 allocations: 8.109 KiB)

julia> n = 10_000; x = randn(n); @time FastTransforms.lib_jac2jac(x, 0, 0, 1, 1);
  0.676580 seconds (7 allocations: 78.344 KiB)