Add normal equation solver for ordinary least squares with not many features. The approach requires one pass to collect AtA and Atb, then solve the problem on driver. It works well when the problem is not very ill-conditioned and not having many columns. It also provides R-like summary statistics.
We can hide this implementation under LinearRegression. It is triggered when there are no more than, e.g., 4096 features.