Alternative to lm() or lm.fit(). Fast matrix operations, no bells or whistles.