I have some code that makes constant use of extrapolate, and while probably not a bottleneck in my use case, I was wondering if it would be possible/desirable to modify extrapolate so that we can pass the buffer used to store the diagonal of the neville tableau. Here is what I get right now:
julia> @btime extrapolate(cos, 1)
73.269 ns (3 allocations: 160 bytes)
(1.0, 1.4518386493023172e-12)
The goal would be to allow e.g. extrapolate(cos,1;buffer=...) to avoid allocations provided buffer is an empty vector of sufficiently large capacity.
I would be happy to put together a PR if this seems desirable.
I have some code that makes constant use of
extrapolate, and while probably not a bottleneck in my use case, I was wondering if it would be possible/desirable to modifyextrapolateso that we can pass the buffer used to store the diagonal of thenevilletableau. Here is what I get right now:The goal would be to allow e.g.
extrapolate(cos,1;buffer=...)to avoid allocations providedbufferis an empty vector of sufficiently large capacity.I would be happy to put together a PR if this seems desirable.