Bsplines are fully non-local

This is easily seen by interpolating a vector with only zeros except for a single one: the basis function is exponentially decaying. If I understand correctly this is not standard Bspline behavior, which are compactly supported. This looks related to the "prefiltering" mentioned in the code (which solves a tridiagonal system and is therefore nonlocal). Is it possible to turn it off?