The current cubic interpolators are unnecessarily re-inverting the coefficient matrix for each cell. This can be precomputed. This would also improve the numerical accuracy of (beyond 5-6 decimal place) by avoiding an unnecessary pre/post normalisation that occurs around the inversion.