Enhanced implementation of Whittaker–Henderson smoothing
