I am not sure what you mean.
If data is equispaced then the interpolation is, as the name suggests, becomes a convolution operation.
Since you did not use any other package, I was wondering how fast your convolution implementation.
Have you compared your implementation of 1D convolution to other implementation in the eco system?
For instance: