Add more methods to DenseMatrix class

Add new methods for LU, LDL factorizations and LU solve.
Implement wrappers for LU, LDL and LU_solve
Add tests
Remove `matrix_c2py()` function