Direct construction of normal equations implemented.  Memory consumptions are significantly lowered for large data set problems.