Implement decomposition based solve functions for solving Ax = b

Removed this since this is not how we are going to do things.