optimal combination forecasts for hierarchical time series .
Algorithm 583: LSQR: Sparse Linear Equations and Least Squares Problems