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)