Implement decomposition based solve functions for solving Ax = b

Why should it use RCP? Given that matrices are mutable, it doesn't make sense to use RCP.