Implement decomposition based solve functions for solving Ax = b

That's one thing, and generally we want to avoid using RCP, unless we have to, as they are slower.