Added functionality for Matrix for finding rref (solving linear
systems), invert matrix, and compute determinant. Only needs LU
factorization now.

Note: None of this is parallelized. I'm just trying to ensure we have
algorithms (we can use this as our course-grained, lock-based
implementation).