Nick

R^n->R interpolation is more difficult and it's been in the back of my head to add it for some time. Do you have an arbitrary pointset in R^n or will...

I think that the Makima spline might be the way to go for this: https://blogs.mathworks.com/cleve/2019/04/29/makima-piecewise-cubic-interpolation/

The main reason this is a good method is that it has a natural generalization to higher dimensions; see Akima's 1974 paper: https://doi.org/10.1145/360767.360779

@tdegeus : I should give you some more context: I have a pybind11-based library which has non-trivial code on both sides of the language boundary and this is distributed as...

@spectre-ns : This will be absolutely amazing!

@spectre-ns : Could not complete the benchmarking: ``` ~/t/build ❯❯❯ ./fftw_vs_xtensor Unable to determine clock rate from sysctl: hw.cpufrequency: No such file or directory 2022-04-30T10:08:54-07:00 Running ./fftw_vs_xtensor Run on (8...

@spectre-ns : Yeah, I added the `-fsanitize=address -fsanitize=undefined` flags. I didn't go searching for these until I hit a segfault though.

``` $ clang --version Apple clang version 13.1.6 (clang-1316.0.21.2.3) Target: arm64-apple-darwin21.4.0 ```

@spectre-ns : Thanks, that fixed it! That check would certainly be nice for foolish noobs like myself. Still have an UBSan error but arising from deep in the library: ```...

Now we have another problem that the complexity is better fit by a cubic in length than NlogN: ``` ~/t/build ❯❯❯ ./fftw_vs_xtensor Unable to determine clock rate from sysctl: hw.cpufrequency:...