Now it is calculated in two stages: first a bspline grid, then
the real gradient is interpolated using the grid.

einspline library yet doesn't compile on Windows, but speaking truly
the whole branch doesn't compile on Windows :) This issue will be fixed
later by porting einspline into c++.