Switching to QR-based (dense) least squares solution saves factor of two speed.  (Not a big deal but I like to understand it by QR.)  Add test code for matrix method.