Matrix multiplication with slicing in gth_solve
Docstrings elaborated