hafen/operator: Loess and STL operator matrix computation for time series

Version 1.0

Computes the operator matrices for local polynomial smoothing methods (loess and STL) applied to time series (or any equally-spaced design setting). Includes capabilities for model selection, prediction, and confidence intervals.