MikaelSlevinsky / FastTransforms

:bullettrain_front: Fast orthogonal polynomial transforms :surfer:

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

Fast jac2jac with integer steps

dlfivefifty opened this issue · comments

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)