Function solve_unit_upper_triangular_in_place 

pub fn solve_unit_upper_triangular_in_place<T: ComplexField, LhsT: Conjugate<Canonical = T>, N: Shape, K: Shape>(
    triangular_unit_upper: MatRef<'_, LhsT, N, N, impl Stride, impl Stride>,
    rhs: MatMut<'_, T, N, K, impl Stride, impl Stride>,
    par: Par,
)

solves $U x = b$, replacing the diagonal of $U$ with ones, and implicitly conjugating $U$ if needed, and stores the result in rhs