* Use `solve_triangular` instead of in `psd_solve_with_chol`

* Add unittest for `psd_solve_with_chol`

* Specify `mode=FAST_RUN` in test

* Relax `test_psd_solve_with_chol` `atol` and `rtol` for half-precision tests