Implement decomposition based solve functions for solving Ax = b
Overload << for DenseMatrix class