Implement decomposition based solve functions for solving Ax = b

Okay, I'll add it.