A fast Gauss-Seidel-like splitting algorithm for convexly constrained spline smoothing