Use “Spline” instead of “Hermite” interpolation for P(k), which
matches GSL cspline better.