We provide an efficient algorithm for computing the smoothing spline score estimator of Cox (1985). The algorithm exploits the banded structure of the linear algebra involved. Calls are made to the LAPACK, Level 1, 2, and 3 BLAS subroutine libraries designed to be efficient on a wide range of modern high-performance computers
