This way the double-to-float conversion is done at a more relevant stage.
uni.size = N;
typename Interpolate::Spline<double, 3, N>::Value value = spline(x);
for(unsigned i=0; i<N; ++i)
uni.size = N;
typename Interpolate::Spline<double, 3, N>::Value value = spline(x);
for(unsigned i=0; i<N; ++i)
- uni.values[i] = Interpolate::SplineValue<float, N>::get(value, i);
+ uni.values[i] = Interpolate::SplineValue<double, N>::get(value, i);