* Use lapack instead of `scipy_linalg.cho_solve`

* Use lapack instead of `scipy_linalg.lu_factor`

* Use lapack instead of `scipy_linalg.solve_triangular`

* Add empty test for lu_factor

* Tidy imports

* remove ndim check