Implement decomposition based solve functions for solving Ax = b

That's fine.