Improve performance of WhittakerSmooth() and AirPls() by using sparse matrix operations. Both benchmarks were done using a spectral dataset containing 33 spectra with 239 variables each. The results show the average time and standard deviation of 70 runs.